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Limitations on Electron Beam Density of Unipotential 
Electron Guns at Low Voltages 


J. Arol Simpson and C. E. Kuyatt 


(Received April 1, 1963) 


The problem of saturating a particular configuration of space with electrons is con- 


sidered. 


It is shown that with unipotential guns of axial symmetry derived from space 


charge limited diodes there exist both a minimum energy and a maximum width-to-length 


ratio of the space that can be saturated. 
flow between concentric spheres. 


These limits are derived for the case of convergent 
Similar limits for partial saturation of the space are 


derived for cases of convergent, parallel, and divergent flow. 


One of the constantly recurring problems in elec- 
tron physics is the generation of low energy electron 
beams in free space. These beams are usually re- 
quired to be circular in cross sections of some specified 
length and diameter and have some energy between 
land 1,000 ev. The highest electron current possible 
within these specifications is usually desired. The 
authors recently gave [1]! a procedure for designing 
such guns using the multistaging principle, i.e., the 
use of an intermediate anode operated far above the 
final beam voltage. Since added complexity and 
practical difficulties, such as the increased probability 
of ionizing the residual gas and increased problems 
of interelectrode leakage, are inherent in such designs, 
it appears desirable to delineate the regions where 
the failure of unipotential guns makes multistaging 
essential. 

It is essential to understand the basic physical 
limitations on achievable electron flow patterns for 
it is the nature of these limitations which determines 
the beam profile and density. These limitations are 
of two classes, one the result of space charge repul- 
sion, the other essentially of thermodynamic nature. 
The resulting restrictions function both individually 
and most importantly in combination. 

Let us consider the usual case figure 1. It is 
desired to irradiate a volume defined by a length 
1 and a diameter 27). Alternately we may consider 
the space defined by / and the angle y. 

As is well known [2], the maximum current, Jax, 
that can be forced through such a space against the 
dispersive forces of space charge is 


Dp 
ce E 
, r 38.51%” tan? y 


| 


38.5 EB 


where current is in microamperes and energies FE. are 
inev. It can be further shown that this results in a 
maximum current density at the point of minimum 
cross section (where the beam diameter is 27)/2.35) 
of 
9 
Tmax => E?? tan? y 


Tlo 


1 Figures in brackets indicate the literature references at the end of this paper. 
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where current density is in microamperes per square 
length unit used in measuring 7. 

To achieve this maximum the beam must be 
launched through a hole the same size as the diameter 
of the space and directed so that in the absence of 
space charge there would be a crossover in the center 
of the volume. We note that while the current 
depends only on the shape but not the size of the 
space this is not true of the density and moreover no 
absolute upper limit on density is predicted. 

One solution to the problem of generating space 
charge limited beams and launching them at the 
proper angle was found by the design techniques 
developed by Pierce [3] and expanded by Samuel [4], 
Field [5] and other more recent workers [6, 7]. 

These solutions consist of a section of a spherical 
diode where the effects of the radial flow beyond the 
beam limits are simulated by proper electrode profiles 
or auxiliary electrodes. When proper account is 
taken of the lens formed by the hole in the inner sphere 
(anode lens) such guns can be designed for currents 
up to the point where the effect of the anode hole on 
the field at the cathode can no longer be compensated. 
For guns without an anode grid this point occurs 
when the hole size approaches the diode spacing or 
at a point where J (u amp) approaches /°” (vy), ie., a 
microperveance of unity. Particularly useful design 
curves for such guns may be found in Spangenberg’s 
work [8]. 

Field [5] seems to have been the first to point out 
that the effects of the thermal energies of the electrons 
at the cathode spread the beam and may be an im- 
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Definition of space to be saturated and of defining 
angle y. 
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portant factor in limiting beam current and density. 
This beam spread is a consequence of a fundamental 
law of optics which in various forms carries the names 
of Helmholtz, Lagrange, Clausius, and Langmuir. 
This law, based upon general arguments of classical 
mechanics and derivable from Liouville’s theorem, 
states, for electrons, that 
Ev? sin 6,d0,dz,= EF}? sin 6,d6.dzz. 

In this case #, and F, are the energies of the 
electrons in planes 1 and 2, @ the convergence angles 
in these planes and dz, and dz, are the corresponding 
elements of length in these planes. 

As written in differential form, this law applies to 
any system whether or not images are formed. It is 
generally not very useful since to determine the 2’s 
and @’s a complete trajectory calculation must be 
performed. If, however, planes 1 and 2 are conjugate 
object and image, or object and crossover then the 
equation may be integrated to give 

EF}? sin 6,=ME?’ sin 6, 
when .M is the Gaussian magnification. 

If we specialize further to the case where plane 1 
is a thermionic cathode and express £, in terms of the 
cathode temperature and .\/f as a ratio of current 
densities we can obtain the Langmuir [9] form: 


J j1+£/(7/11,600)] sin*6, ~J,(11,600F/T) sin 74. 


(3) 


The general consequences of these equations are well 
known and are treated at length in the standard 
texts [10]. 

If we compare eqs (2) and (3) we see that @ priori 
there is no reason why he density required by eq (2) 

cannot be obtained subject to eq (3) since J, is not 

Senlied. In practice of course there are material 
limits on the cathode, but it can be shown that for 
energies below a few hundred volts modern cathodes 
are sufficiently “bright” to satisfy almost any 
demand. 

However, there is a third restriction on unipotential 
guns which in certain voltage regions prohibits the 
densities required by eq (2) to saturate a given space. 
This restriction is the bound on obtainable J; which 
is set by space charge within the gun structure. 
Meltzer [11] has shown that a limitation of the form: 


kE 


J=~— 


(4) 
where S a linear dimension is a general characteristic 
of space charge limited diodes [12]. If we now 
consider a gun based on such a diode rp) can serve 
to determine the scale and play the role of Sin eq (4). 
The value of the multiplier & is set by the overall 
geometry necessary to give the convergence angle y 
in free space beyond the gun. Hence k is a function 
of y. There are, of course, many conceivable guns 
which give an angle y but for each general class k 
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remains a function of y. The result is a limitation 
on cathode current density of the form: 


Ky) E*” 


Yo 


x a (5) 


Substitution of this value of J) for J; in eq (3) 
gives the maximum obtainable current density, Jo, 
in the space to be saturated. If J, is less than Jax 
from eq (2), the space cannot be saturated. Combin- 
ing the equations shows that if 
eG oe tan *y ‘2 
E<— —— = — (6) 
x 11,600 /(y) sin *@, 
cannot be saturated. The connection 
between y, the launching angle, and 6, the angle at 
minimum beam density, is somewhat subtle involving 
as it does the detailed nature of the trajectories 
between the anode and the point of minimum beam 
cross section as well as the angle and energy with 
which a given electron leaves the cathode. However, 
without this detailed knowledge one fact is clear: 0, 
cannot exceed y without some electrons passing 
outside of the space to be saturated. Hence, we can 
set 6, equal to y without violating the inequality. 

There is moreover an upper limit on y since the 
anode hole lens of a diode always divergent. 
Hence electrons make a larger angle with the axis 
just before the lens than they do after the lens. 
Considering a class of guns which give saturation 
and taking 90° as a natural limit of the angle inside 
the lens we see that there is always a maximum 
value of y for which the space can be saturated. 
Actually the limit is reached at smaller angles since 
a large y requires a higher perveance which in most 
structures requires closer anode cathode spacing re- 
sulting in higher fields inside the anode and hence 
stronger divergence at the anode lens. 

For guns based on radial flow between concentric 
spheres the quantity f(y) in eq (6) may be obtained 
from the well-known optics of Pierce guns [8]. 
Pierce [14] has calculated the values of the gun 
parameters where the space charge spread equals 
the thermal spread. In our notation this procedue. 
is equivalent to the assumption that y equals 6. 
Danielson, Rosenfeld, and Saloom [6] have done 
careful numerical calculations on the trajectories in 
a series of Pierce guns and obtained the final size 
of the minimum beam cross section when both 
space charge and thermal spread are important. 
This procedure, of course, accounts fully for the 
detailed relationship between y and 6. Over the 
range of gun parameters studied, there is no doubt 
that the latter is the more precise calculation and 
should be the final check on any proposed design. 
However, both calculations focus their attention on 
the gun and require the calculation of a series of 
guns each of which will saturate the space. It is 
more useful, especially in the early stages of design 
of an experiment, to concentrate on the space so 
that the feasibisity of the experiment may be judged 
prior to detailed gun design. 


the space 


is 





To achieve this end we have evaluated eq (6) 
using an f(y) derived from the results of Pierce. 
We assumed a current efficiency of 0.9 and a corre- 
sponding intensity efficiency of 0.39. The anode 
lens was assumed to be of the Davisson [15] form, 
the ~10 percent correction of Danielson et al., not 
being significant for our purposes. We assumed y 
equal to @. Our calculations are hence much in the 
spirit of those of Pierce but extended to account 
more fully for the effect of the anode lens. They 
are somewhat less optimistic than his, but still must 
be considered outer bounds of possible operation. 

The limits found are independent of the scale of 
the space, depending only on y. If the seale is re- 
duced beyond a critical minimum, it is possible that 
the cathode material limit may be reached. How- 
ever, if a cathode-anode spacing of 1 mm is con- 
sidered a practical lower limit 1 amp/em? cannot be 
exceeded below 265 v. 

The results of the calculation are shown in figure 2. 
On this plot lines of constant perveance, //V*?, ap- 
pear as straight lines sloping upward to the right. 
When the space is saturated, as in the upper portion 
of the graph, eq (1) gives the relation between the 
perveance and corresponding value of y. In this 
portion of the graph the unity perveance line corre- 
sponds to a y=9.2°. The area to the left of this 
line is generally not accessible to unipotential guns 
because of their limited perveance. In the area 
above the almost horizontal line J//J_.,=1, spaces 
characterized by gammas up to 9.2 deg may be 
fully saturated. Below this ///,,..=1 line the spread 
of the beam due to thermal velocities dominates the 
situation and frustrates attempts to operate at lower 
voltages by reducing the current drastically below 
the value predicted by extending the lines of con- 
stant y. It will be noted that there is no unipoten- 
tial gun which will saturate a space below ~220v 
beam energy. 

If one is prepared to relax the current requirement 
so that only 0.5 of the maximum possible current is 
passed, the guns need be less convergent and the 
area between the perveance=1.0 curves now marked 
y=12.8 and J/Jyax=0.5 is available. The entrance 
diameter (anode hole of the gun) instead of being 
equal to the space diameter is now found to equal 
0.52 of this diameter. 

If one still further reduces the current requirement 
to 0.25 I,,.; it is possible to use guns whose injection 
angle into the space is 0°. Under this circumstance 
the anode hole diameter should be 0.425 of the tube 
diameter. The details of the calculation are slightly 
different since now the maximum density lies within 
the gun structure but the results are of the same 
form. The minimum usable voltage is decreased to 
10 ev and spaces up to y=17.85° may be one-quarter 
saturated. 

A sort of natural limit to this process is reached at 
the dotted curves above and to the right of which 
the highly divergent parallel diode guns are usable. 
Although such guns can be used to a minimum of 
almost 1 ev and to space angles of y=74° (we 
assumed a unity spaced diode of microperveance 
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Figure 2. Graph of current versus voltage showing regions 
accessible to guns capable of saturating a given space with I 


max- 


Cases of subsaturation, 0.5 Imax, 0.25 Imsx and 0.05 Imax are shown. The usable 
range lies to the right and above of all curves. On the right margin are given 
the ratio of cathode/anode radius of the optimum Pierce gun. 


2.33 whose anode diameter is very much smaller 
than the tube diameter) the fact that they will only 
put 0.05 Imax through any given space makes them 
not very attractive. 

The curves of figure 2 show clearly that the range 


| of applicability of unipotential guns lies at compar- 


atively high energy and wide angles. To operate 
outside these limits it is necessary to effectively 
decouple the thermal energy and space charge limit 
by adopting the ei c principle and operating 
a diode or triode at comparatively high voltage 
and decelerating to ‘the desired low energy. 
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New Radiofrequency Mass Spectrometer Having High 
Duty Cycle’ 
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A new type of time-of-flight mass spectrometer has been developed which, unlike most 
other such instruments, utilizes a large fraction of the nonpulsed stream of ions. The ions 
are accelerated through a d-c voltage and then pass into a drift tube. A radiofrequency 
voltage is applied at the exit and entrance gaps of the drift tube in such a way that the field 
at the entrance is the negative of that at the exit. The number of ions which pass through 
the drift tube with such a velocity that the energy change from the rf fields is a minimum, is 
measured by placing an electrical barrier in front of the collector and differentiating the 
cutoff curve. Velocity modulation of the beam resulting from the rf field at the first gap 
limits the instrument to moderate resolution. However, the instrument has a 60 percent 
duty cycle, a value significantly greater than that achieved by most other time-of-flight 
mass spectrometers. It also has the advantages of electrical and mechanical simplicity. 
Due to the nature of the instrument’s output, it will be most useful in applications requiring 
the analysis of a simple mass spectrum and where high duty cycle is of real Value. The 
instrument has resolved the isotopes of rubidium at masses 85 and 87, and has identified 


one component ion beams ranging in mass from helium to cesium. 


1. Introduction 


This paper describes a new type of time-of-flight 
spectrometer which has a large duty cycle along 
with moderate resolving power. It was developed 
for use in flame research investigations at the Na- 
tional Bureau of Standards. Two types of experi- 
ments have been planned. The first is that of 
measuring the appearance potential of certain 
halogenated compounds. In experiments of this 
type, one looks for the appearance of a single ion 
peak as the energy of the ionizing electron is in- 
creased. High sensitivity is important if the ap- 
pearance potential is to be determined accurately. 
The second application is that of analyzing ions 
drawn directly from a flame. Other laboratories 
that have used flames as the ion source have ex- 
perienced difficulty in drawing enough ions from the 
flame for analysis. The high duty cycle of the mass 
spectrometer described in this paper should make the 
instrument efficient in using the ions which are 
available in these two experiments. 

Most of the previous efforts to increase the sensi- 
tivity of mass spectrometers have been directed at 
the collector end by using electron multipliers. An 
electron multiplier suffers from a changing gain over 
a period of time and was not used in the instrument 
described here. However, it would be possible to in- 
corporate one if an even greater sensitivity were 
desired. 

In general, 
eters depend on the fact 
heavy ions take a longer 
distance than do light ions. 
ion varies according to: 


time-of-flight (T.O.F 


.) Mass spectrom- 
that, given equal energy, 
time to traverse a given 
The velocity of any 


* A portion of the funds for support of the work reported was provided by the 
Bureau of Ships, Department of the Navy. 
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/2qV 


\ m 


where qg is the charge of the ion, V is the total voltage 
through which the ion has been accelerated, and 

is the mass of the ion. If a pulse of ions of various 
enters a field-free drift tube, after being 
accelerated, the ions will travel in bunches corre- 
sponding to their mass. ‘The distance between bunches 
will increase the further the ions travel down the 
drift tube. The mass spectrum can then be meas- 
ured by observing the ion current collected at the 
exit end of the drift tube as a function of time. 

All instruments of this type have to be pulsed at 
the source and the time between pulses must be 
great enough to clear the drift tube of ions from the 
previous pulse. The pulsing greatly reduces the 
duty cycle of the mass spectrometer which thus 
lowers its time average collector current. The 
instrument of this type which is best known was 
reported by Wiley and McLaren in 1955 [1].’ Ions 
are created by a 0.25 usec electron pulse once every 
100 psec [2], giving a duty cycl le of 0.25 percent. 

Another type of T.O.F. instrument which is 
described in the literature [8, 4| distinguishes masses 
by the energies the ions receive from a series of rf 
fields. A continuous beam of ions passes through 
10 or more electrodes which are perpendicular to 
the beam path. The geometry of these electrodes 
and frequency of the rf fields maintained between 
them are such that ions of one mass receive more 
energy than lighter or heavier ions. An electrical 
barrier is placed in front of the ion collector. Only 
ions which have received extra energy from the rf 
fields are able to overcome the barrier and reach the 
collector. This instrument gives good resolution, 


(1) 


masses 


i Figures in brackets indicate the literature references at the end of this paper 





but has several drawbacks. The spacing between 
the electrodes is critical for proper operation 
which makes construction and modification difficult. 
The numerous electrodes through which the ions 
must travel reduce the transmission of the instru- 
ment. And, finally, this type of instrument collects 
pulses of ions even though a continuous beam of 
ions enters the analyzing region. This effect is 
due to the fact that an ion of correct mass must 
also have the correct initial phase with respect. to 
the rf voltages as it enters the fields in order to gain 
the necessary amount of energy to reach the collector. 
Thus the sensitivity of the instrument is reduced 
in the same way as with the pulsed type of spectrom- 
eter. The duty cycle is given in reference [4] as 
» where N is the number of stages and Ff is the 
Rs v3 
resolution. An instrument of this type having 20 
stages thus has a duty cycle less than 3 percent at a 
resolving power of 50. 


2. Theory of Operation 


The mass spectrometer reported here resembles the 
last instrument described in that ions of different 
masses are distinguished by their energies rather 
than their spatial positions. It differs in that 
initial phase requirements are much less severe. 

Consider a drift tube of length ZL which has rf 
electric fields applied at its ends as shown in figure 1. 
Tungsten screens, not shown in the figure, are 
attached to the ends of the drift tube and define the 
entrance and exit. It is seen that the grids near 
each end of the drift tube are at the same potential. 
Therefore, at a given instant, the rf fields between 
the drift tube and the grids are exactly equal and 
opposite at the entrance and exit ends of the tube. 
These rf fields are represented by the upper and lower 
curves In figure 2. 

The operation of the instrument can be explained 
most easily if the flight time of a given mass through 
the drift tube is assumed to be constant. This is 
not entirely a valid assumption since the effects of 
the variation in flight time (velocity modulation) 
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FIGURI 1. Simplified diagram of the’ mass spectromete Fx 
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Figure 2. Radiofrequency fields at the entrance and exit of 


drift tube. 


A resonant ion entering at time A: and leaving at 
in energy A typical nonresonant ion entering 
time B: experiences a net gain in energy. 


12 experiences no change 
t time B,; and leaving at 


resulting from the rf voltage at the entrance of the 
drift tube are quite noticeable. Nevertheless, 
velocity modulation will be neglected until the 
resolution and duty cycle of the instrument are 
discussed. 

Suppose an ion enters the drift tube at some 
time A,. If the frequency of the rf voltage and the 
velocity of the ion are such that the ion emerges at 
the other end of the drift tube after one period of 
the rf wave, the ion will pass through the second 
rf field at time A». It is seen that no matter what 
part of the cycle the ion enters the drift tube, the 
energy which it obtains at the first field is exactly 
balanced out by the energy it loses at the second. 
Consider now, a heavier ion which takes a longer 
time to traverse the drift tube. For example, if it 
enters the drift tube at time B,, it might pass through 
the exit field at time B,. Any ion which does not 
pass through the drift tube in an integral number 
of cycles will undergo an increase or decrease in 
energy, depending upon when it entered the drift 
tube. A beam of such ions which is homogeneous 
in energy at the entrance of the drift tube will 
emerge with a distribution of energies. However, a 
similar stream of ions which passes through the 
tube in an integral number of rf cycles will emerge 
with all the ions unchanged in energy. These are 
the ions which the detection system measures. 
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3. Detection System 


As stated before, the detection system must 
identify the resonant ions by their energy. This 
measurement is made by the retarding field method 
[5, 6]. An electrical barrier is placed immediately 
in front of the Faraday cage collector. For a given 
barrier potential, “, measured with respect to the 
ion source, Only ions with energies equal to E electron 
volts or greater can pass. Or 


S(f) N(e)de (2) 
JE 


where S(/) is the number of ions per unit of time 
that penetrate a barrier of height E and N(e) is the 
number of ions approaching the barrier per unit of 
time with energy e at the barrier. N(e) 1s obtained 
from (2) by diffe ‘rentiating with respect to L 


dS(E) 


dE —N(E). (3) 


Since the detection system must measure the number 


of ions which have had their energies unchanged by 
the rf fields, and, since these particular ions arrive 
at the barrier with the same energy that they had 
at the source, V(0) is the quantity sought. Figure 
3 is a set of cutoff curves showing collector current 
as a function of barrier height in volts, obtained 
with an z-y recorder. The center curve is the. cut- 
off curve for a potassium beam which is in resonance 
with the rf voltage, taken with n equal to 12. Here, 
n is equal to the number of rf cycles which an ion 
takes to traverse the drift tube. The other curves 
were taken with n equal to 11.5, 11.7, 12.2, and 12.4, 

The differentiation represented in “eq (3) is per- 
formed in the instrument by introducing a variation 
in E about the point H=0 shown as AF in figure 3. 
The resulting a-c component in the collector current 
is proportional to dS(/). And dS(£) is proportional 
to N(e) if dE is infinitely small. A square wave of 
magnitude +3.9 v was used to produce this varia- 
tion and this value will be assumed in the remainder 
of the paper. Using a dE of finite size means that 
an average value of N is actually measured over 
the range H=0+3.9 v. If the ion beam is com- 
posed of more than one mass, the resulting a-c 
signal is the sum of the contributions from the 
individual masses. 

The variation in the ion current was amplified 
by a-c amplifiers and then fed to a synchronous 
detector, AE supplying the reference phase used 
in the detector. An z-y recorder was connected 
to the output of the detector. The masses were 
swept by changing the frequency of the rf voltage. 
A potentiometer was attached to the tuning knob 
of therf generator. The voltage from the potentiom- 
eter was thus related to the rf frequency and was 
connected to the other input of the z-y recorder. 
The potentiometer and tuning knob were driven 
by a 1 rpm motor. 


Figures 4, 6, 7a, 7b, and 8 were made using this 
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Figure 3. Cutoff curves obtained from an x-y recorder for 
n=12.0 (central curve) n=11.5(1), n=11.7(2), n=12.2(3), 
and n=12.4(4) showing the difference in output signals for 
resonant and nonresonant ions. 


A square wave was used for AE instead of the sine wave shown 


The abscissa is not linear in fre- 
quency. No development work has been done to 
make this instrument a fast-sweep instrument, 
capable of displaying its output on an oscilloscope, 
but this seems feasible. 

The advantage of picking out the ions with zero 
energy change should be noted. Instruments of 
the type in reference [4] must maintain a constant 
ratio of barrier height to rf voltage to obtain the 
desired resolution. “However, if the number of ions 
with zero energy change is measured instead of the 
number with maximum energy change, the strong 
dependence on resolution no ‘Tonger exists and less 
care need be taken in controlling the rf amplitude 
as the rf frequency is changed. 


arrangement. 


4. Harmonics 


As has been mentioned before, ions may stay in 
the drift tube any whole number of cycles in order 
to be detected. “Thus, as the frequency is varied, 
the mass peaks will repeat corresponding to various 
values of n as in figure 4. The two rf frequencies 
at which adjacent harmonics of a particular ion 
appear must be measured to determine the ion’s 
mass. The ratio of these two frequencies is equal to 
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FIGURE 4, 


Che eight through twentieth harmonics of K 


the ratic of their harmonic numbers. That is: 


(4) 


This equation can be solved for n, or n can be 
found conveniently on a slide rule. The procedure 
is to look for two consecutive integers which have 
the same ratio as /,/f, Knowing » and the 
corresponding frequency, f equal to 1/T,,), the flight 
time of the ion, ¢t,is found fromt=nT,. If the length 
of the drift tube is LZ, eq (1) can be solved for m, 
using v= L/t. 

The determination of m becomes difficult 


if the 
beam contains ions of many masses. 


The harmonics 
of the mass peaks become too cluttered to pick out 
two adjacent harmonics. In applications requiring 
the study of complex beams, it would be desirable 
to eliminate all but one of the harmonies for each 
mass. Partial elimination of the harmonics has 
been achieved by applying two radiofrequencies, 
referred to as f; and f,, as in figure 5. Frequency 
fi is generated from f, by frequency division so that 
1 


the ratio of f, to f; is a known whole number, say r. 








Figure 5. Drift tube arrangement used to eliminate harmonics. 
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The a-c collector current as a function of rf fre quency as plotted by an X-\ recorder. 


ind K 


re shown. 


With the present equipment, r can be chosen to 
be equal to 2, 4, 8, 16, or 32. In order for the ions 
to escape a change in energy, they must pass through 
two drift tubes emerging from each tube after an 
integral number of rf cycles. Let n, be the number 
of rf cycles an ion takes to travel through the first 
drift tube, and n, the number of rf cycles an ion takes 
to travel through the other tube. Since both tubes 
are of equal length, the transit times for both tubes 
must be equal. That is, n;7;—m7). By defini- 
tion, f;=1/7T; and f,=1/T). 
Then 


If n, is equal to one, nm, must be equal to r. If n;is 
equal to any integer, k, m, must be equal to kr. In 
practice the possibility of k being greater than one 
for a particular mass is eliminated by sweeping over 
an rf frequency range of less than one to two. Thus, 
an ion which would pass through the second drift 
tube with n,#rdoes not pass through the first drift 
tube ina whole number of rf cycles. The voltage at 
frequency f/f; discriminates against the unwanted 
harmonics while the voltage at frequency /, resolves 
the masses as before. As is the case in using a 
single drift tube, each rf voltage can be increased 


until the improvement in resolution is balanced 


| out by the adverse effects of velocity modulation. 


Figure 6 was obtained by using this principle of 
harmonic elimination. It shows the frequency range 
where the 10th through 28th harmonics of potas- 
sium would normally appear. The largest peak at 
the center is the 16th harmonic of K*. The two 
small peaks on either side are the 15th and 17th 
harmonics of K*. Figure 6 was made using the 
maximum low frequency voltage which the equip- 
ment could generate. 
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Ficure 6. 


rn Na, A i ne V 
Awan 
in x-y recorder plot of the a-c collector current as a function of rf frequency using 
the harmonic elimination equipment. 
¢ large peak is the 16th harmonic of K®, The smaller peaks on either sid 


ure the 15th and 17th harmonic of K 


n=30 
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| 


} 
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Figure 7. The 


a-c collector current as a function of rf frequency showing the harmonics of 
Rb® and Rb* (a) and of cesium (b). 


The plots were obtained by using an z-y recorder 
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FIGURE 8. An x-y 


Shown are harmonics of K**, K*!, Cl 


5. Resolution and Duty Cycle 


In order to find the resolution and duty evele 
which can be expected from this instrument (with- 
out harmonic elimination), the effect of the rf 
voltage at the entrance of the drift tube on the flight 
time of the ion must be considered. Neglecting 
initial energy, the energy of an ion after passing 
through both rf gaps can be expressed as: 

qi VactVinF|sin wr—sin w(7+t)| (6) 
where 

Va is the d-e accelerating voltage, 

Vr is the amplitude of the radiofrequency voltage, 

tis the time of flight of the ion through the drift tube, 

ris the time at which the ion enters the drift tube, 

w is the angular frequency of the rf voltage, 

F is the transit-time factor due to the fact-that 

the gap widths are of finite size. 

F can be shown to be very nearly equal to unity by 
a calculation similar to that in reference [7]. Under 
the conditions that the gap width is 2 mm, Z is 
1 m, n=20, and that Vi.<<V,, F is found to be 
0.997. Since the above conditions are valid approxi- 
mations for this instrument, F will be assumed to 
be unity. 


By writing ¢ ; in eq (1 


recorder plot showing a-c collector current as a function of rf frequency. 


and Cl 


where ¢, is the flight time of an ion which enters the 
drift tube when the rf voltage is zero and is equal to 


m 


i= = (8) 
\ 2q\ “ 


Substituting (7) in (6), the net energy in electron 
volts which the ion receives from the rf fields isJseen 


to be: 


; SF “t 
€re=V rt { sin WT sin] wr t wty ( 1+ V.. sin wr ) }} 


(9) 


Thus ¢,, is not zero under the condition that f) is a 
whole number of cycles as was assumed to be the case 
when velocity modulation was neglected. 

If the magnitude of e,; is larger than 3.9 ev, the 
ion will not be represented in the collector output 
signal. The theoretical peak shape was found by 
plotting e,; as a function of wr between zero and 2x 
for various wf). Then for each value of of) the 
fraction of the cycle over which —3.9<e,. <+3.9 
ev was measured, and this fraction was plotted 


. . » lo = 4 7 
as a function of —— eas figure 9 over the range 19.5 
wt, 


20.5. 
») . 


Values of V,, and V,. used were the 
Qn 

measured values when f, equaled 20 times TJ for po- 
tassium (K*): 24 volts and 3000 volts, respectively. 
It is seen in figure 9 that as a result of velocity 
modulation (a) the maximum duty cycle does not 


Wy 0: 
actually occur when =~, is a whole number as was 
¥ T i 





the case when velocity modulation was neglected, 
and (b) the duty cycle at its maximum is equal to 
61 percent. 

The duty cycle actually obtained under the above 
conditions was measured by presetting the rf fre- 
quency for maximum collector output signal. The 
ratio of the a-c collector current with rf voltage on to 
that with rf voltage off was found to be about 60 
percent. 

For a given radiofrequency, 
9 may be converted to mass numbers by 


oh ty 
With * 


corresponds to mass number 41. 

Up to a point, the difference in e,, between two 
ions slightly different in mass increases as the time 
rate of change of the rf voltage increases. The time 
rate of change can be increased by increasing the rf 
frequency or by increasing the magnitude of the rf 
voltage. However, increasing the rf voltage increases 
the amount of velocity modulation, and it can be 
seen in figure 9 that velocity modulation causes a 
broadening in the peak shape. Increasing the fre- 
quency of the rf voltage increases the adverse effect 
of whatever velocity modulation is present. It 
thus reasonable to assume that there exist values of 
V., and w» which give a maximum resolution. 

An analysis of the peak shape was made only for 
the values of V,, and wt) given above. These values 
were chosen because they produced good resolution 
between K* and K* (see fig. 4). It is possible that a 
further study might show that other values would 
result in better resolution. 

It is assumed in figure 9 that 
energetic at the source, (b) they do not lose energy by 
collisons they travel through the instrument, 
(ec) they are in perfect focus at the barrier [7], and 
(d) ions with insufficient energy are cut off with 
perfect sharpness at the barrier.2 The center curve 
in figure 3 indicates a final energy spread in the order 
of 30 electron volts, probably due mainly to the 
failures of assumption (a) and (d). Failure of any 
of the above assumptions results in a lower current 
efficiency, as well as a broadening of the peak shape. 
A careful study of these factors has not been made. 


6. Results 


The following positive ions have been detected 
with this instrument: K*, KK“, Na®*, Het, Nj, Nt*+, 
CE, Cl, Rb®, Rb‘, and Cs", All except the he- 
lium and nitrogen ions were obtained by placing 
salts on a hot tungsten wire filament [9]. Helium 
and nitrogen were ionized by heating a bare tungsten 
filament in the respective atmospheres at reduced 
pressures. 

The d-c accelerating voltage was maintained at 
approximately 3,000 v, and the drift tube was roughly 
1.12 m long (0.5¢ 
monic eliminator). 
sity of measuring V 


the abscissa of figure 
(8). 


wo 
—= 20.5 
2r 


eq 


equal to 20 for mass number 39, 


is 


(a) the ions are mono- 


aus 


the neces- 
-alibration 


In order to avoid 
and LZ accurately, a 

2 For sharp cutoff, the barrier must be such that no field can penetrate through 
it. The collector must be able to collect ions of large and small energy with equal 


efficiency and be able to prevent ions which have experienced elastic reflections 
from escaping the collector. 


| 





} m long when used with the har- | 
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FiGuRE 9. 


Computed peak shape, (a) with velocity modulation, 
and 


(b) when velocity modulation is neglected. 
If w is set so that t.=20T' for K®, at./2x=20.5 corresponds to K“!, 
] ym 

Do Vb 
each of the ions listed above. 
equal to 0.678+0.003 (~0.4%) 
atomic mass units (O” 
megacycles. 

This figure was obtained by using known values of 
mass for the above ions. The frequencies at which 
the various harmonics appeared were measured 
with a timer-counter. The percent fluctuation of 
the calibration factor is an indication of the pre- 
cision of which the instrument is capable. With 
the present source, it is necessary to apply a focusing 
field in the region where the filament is located. 
The proper focusing voltage changes each time the 
filament is removed to replace the salt. It has been 
found that this change in focusing voltage produces 
a slight change in the calibration factor, accounting 
for a portion of the 0.4 percent variation. For a 
fixed value of the focusing voltage, the calibration 
factor is found to vary about 0.1 percent. Figure 
7a shows the 30th through 48th harmonics of Rb* 
and Rb‘; figure 7b shows the 15th through 37th 
harmonics of Cs". Figure 8 is a mass spectrum 
showing K*, K*!, Cl®*, and Cl’. Some variation in 
the peak heights of the harmonics can be observed 
in figures 7b and 8. Duplicating these graphs 
shows that this fluctuation is due to the instability 
of the source. Little care has been taken, up to 
date, to ensure linearity in peak heights. 

In figure 10, the ion source at the right, the tube 
containing the drift tube, the collector shielding box, 
and some of the electronic equipment are visible. 

Because of the cluttering of the harmonic peaks, 
this instrument will be useful in experiments requir- 
ing the analysis of ion beams containing a small 
number of mass components such as those mentioned 
in the introduction. In experiments where the clut- 
tering not a serious limitation, the instrument 
offers important advantages. In addition to the 
high duty cycle, it has light weight, and electrical 
and mechanical simplicity. There are no critical 
alinements or spacings. 


actor d efined as CF = was measured for 


It was found to be 
if m is expressed in 
16.00) and 1/7 in 


is 
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A Power-Series Buildup Factor Formulation. Application 
to Rectangular and Off-Axis Disk Source Problems’ 
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The response (e.g., dose-rate) of an isotropic detector to primary radiation from 
finite plane source may, for points less than a mean-free-path distant, be evaluated as the 
This series is derived by expanding the exponential dependence 
wor) a8 a power series in wor and integrating, numerically or analytically, each term over 
If the medium is effectively infinite and homogeneous, the scattering 


sum over an infinite series, 
exp | 
the source array. 


a 


properties of this medium can be characterized by a point isotropic source buildup factor. 
Buildup factor data are often approximated by formulas having simple analytic dependences 


OD Mor 


and numerical parameters independent of yor. 


Any such set of parameters can be used 


to generate a set of weight functions 6, for an infinite series buildup factor representation 


exp ( 
n 


Which can be used with the above primary radiation series solution to give, instead 
series solution including both primary and scattered radiation. 


buildup only 


tT Mor) i > Dn 2% 


a 


Tables of b,’s 


(air-dose 


derived from coefficients of cubic polynomials fitted to the Goldstein- Wilkins 


data (N YO-3075) are given for 0 <n <13 at primary photon energies of 0.5, 0.7, 1.0, 1.5, 2.0, 


d.09, Dd, 


7.5, and 9.5 Mev in water, Al, Fe, Sn, W, Pb, and U. 
a comparison with dose-rates measured by a detector separated from a Co-60 rectangular 
plaque food irradiator by a layer of steel and a layer of water and (2) 


The method is applied to give: 


dose-rate profiles at 


constant heights across a cleared circular area in a Co-60 infinite plane source in air, 


1. Introduction 


A number of methods for evaluating radiation 
fields from distributed sources have been developed 
as part of a general program at the National Bureau 
of Standards [1, 2]' and other agencies. Some of 
these methods, which provide easy and rapid solu- 
tions to radiation shielding problems, rely on the 
approximate proportionality between dose-rate and 
the solid angle subtended by a uniformly radiating 
surface. This approximation has been shown to be 
quite good if the material interposed between source 
and detector is of the order of a mean free path in 
thickness [3]. 

For situations involving much less than a mean 
free path of interposed material, as in the case of a 


detector separated from a surface deposition of | 


radioactive material by only a thin roof or a few feet 
of air [4, 5], the power-series method developed by 
Sievert [6] for circular disk sources has been extended 
(7, 8] to include rectangular sources and to take into 
account point-source buildup factor data in poly- 
nomial form [9]. 


*Work supported jointly by the Office of Civil Defense (DOD), the Navy 
Bureau of Yards and Docks, and the Defense Atomic Support Agency. 
i Figures in brackets indicate the literature references at the end of this paper. 
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In addition, it has been reported [10] that buildup 
data from a number of other analytic formulations 
[11-15] can also readily be transformed for direct 


use in a Sievert-type power-series solution. The 
present work constitutes a published account of the 
material reported in [10] and includes a quantitative 
comparison of the buildup data formulations con- 
sidered [9, 11-15] for Co-60 gamma rays in aluminum 
for zero to ten mean-free-path penetrations, a table 
of buildup data transformations, a sample table of 
transformed buildup data, and sample calculations 
and results. 

Besides fallout radiation studies [5], the power 
series method is applicable, or has already been ap- 
plied [16] to such diverse distributed source studies 
as: 

(a) Tracer studies involving an elongated rec- 
tangle as the initial distribution of radioisotope- 
tagged sand in a beach-erosion sand-drift experi- 
ment [16], 

(b) food irradiator design studies involving plaque 
sources [3, 17, 18], and 

(c) calculation of the dose distribution in tissue 
adjacent to disk or other shaped radioactive appli- 
cators for treatment of malignancies [6, 19]. Also, 
buildup data in the form of the series coefficients 
presented here are directly applicable to a method 





of barrier shape optimization described by Jain and 
Sharma [20] who use, in effect, the first three terms 
of this series for their analysis. 


2. Buildup Factor: Definition and Discussion 


As in references [9] through [15], the source and 
detector environment is idealized throughout this 
paper as a homogeneous infinite medium, shown 
schematically in figure 1. In such an idealized 
situation the scattered radiation reaching the de- 
tector can be taken into account by means of a 
point isotropic (PTI) source buildup factor [21] 


Berm (Eo, Z, r)=(D°+ D*)/D (1) 


which depends on the effective atomic number Z 
characterizing the medium and the distance r be- 
tween the detector and a monoenergetic point 
source of photons of energy ? &. 

The quantities D® and D* in (1) refer to the un- 
scattered and the scattered components, respec- 
tively, of the radiation received by the detector. 
For a given physical situation, D® and D* (and 
hence Bpy:(Eo, Z, 7)) can assume a variety of nu- 
merical values depending on what kind of detector 
response (e.g., photon number flux [22], energy 
current [22] at an interface, energy dissipation in a 
medium- or air-equivalent detector, etc.) is used to 
describe the radiation field. The buildup factor 
data in this work correspond to a detector which 
measures radiation flux in terms of energy dissipation 
in air. In this case Bpy (2, Z, r), D® and D* in (1) 
are the same as the quantities By, Ja and IJ@ 
defined in [21]. 


2 Here, and in what follows, the source is considered to be either monoenergetic 
or an effectively monoenergetic component of a photon energy distribution over 
which the final results must be summed. 
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Figure 1. Point isotropic source geometry. 


| energy dissipation in water (D 


The total detector response )) to a point isotropic 
source is 


D=DP+D'=D Bp (Ep, Z,1r) 
(k 4m)(exp ( — po?") | Berr( Eo, Z, r) (2) 


in which pol=“(, Z)| is the narrow-beam attenua- 
tion coefficient [23] at photon energy , for the 
medium. The source strength factor k, for the de- 
tector response-type (= J,) specified above is 


k=nX Ey X uen(E, Z) X (dimensional factors) (3) 


where nv is the total number of photons of energy EZ 
emitted per second by the source, p,»(/o, Z) is the 
energy absorption coefficient [24] for air,> and the 
‘dimensional factors’’ provide the units [e.g., (rads/ 
hr)/curie}] desired in the final answer. 


3. Analytic Formulations for Bp7;(Eo, Z, r) 


For applying buildup data to distributed source 
problems, the dependence of Bpy(, Z,r) on r is 
usually approximated by some analytic function. 
This not only provides some economy in data tabu- 
lation, but the function can be so chosen that its 
presence entails little or no additional complication 
in the integration 


. 


D (o 41) lexp( — pio?” ) r|Bem(Eo, Z, r)dS (4) 
JS 


of (2) over a distributed source S with differential 
elements dS, as schematized in figure 2. The con- 
stant factor in (4) 


o=dk/dS (5) 


3 Note, however, that example 1 in section 7 requires an answer in terms of 
Ip) rather than in air (D=TJa), so, that pen (Eo,.H20) 
is used in the computation in 7.1, rather than wen (£o, air) corresponding to the 
data in table 3 in section 5. 
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FicgureE 2. Finite plane isotropic source geometry. 





is the source strength per unit of source area or 


volume, with & defined as in (3). 
rT’ . . 
lhe following functions (a) to (f) have been used 
to represent point isotropic source Bpy;(p, Z,r) data 
in order to obtain analytic solutions to (4). 
(a) Linear approximation [11] (one parameter): 
Bern (Eo, Z, r) >| +- Apo?” (6) 
where 
a (uo— Men) Mens (7) 


(b) Exponential-linear empirical fit [12] * (two 
parameters ): 


Bor (Eo, Z, r= l +o? exp ( —Dyor). (8) 


(c) Exponential-polynomial fit [13] (one or three 
parameters ): 


Bory( Ee. Z r) 


‘ exp (+-por) (br+-1) (9a) 


or, for larger thicknesses: 


Ber (Eo, Z, 7) =exp (+ mor) (ar?+br+c). (9b) 


(d) Low 


parameters) : 


order polynomial fit [9, 25] (three 


Bory ( Eo, Z, 7) =1+-Bywor +B2(uor)? +B3(uor)®. (10) 
(e) Sum-of-exponentials fit [14] (three parameters) : 
Bprs (Eo, Z, 7) =A exp (—ayuor) +(1—A) exp (—aouor) 
(11) 


TABLE 1. 





The point isotropic source buildup factor (“exposure dose”’ 


(f) Exponential-linear moments calculation result 
[15] (four parameters): 


Bpri (Eo, Z, 7) =1+- mor { A,B, exp [(1— By) wor] 


+ A,B, exp [(1—Bz2)uor|}. (12) 

Table 1 contains data calculated at uyr—1, 2, 4, 7, 
and 10 mean free paths according to approximation 
(a) above, empirical fits (b) to (e) and direct-result 
parameters (f), in addition to earlier calculated 
results from NYO-3075, [26] (g), all for “exposure 
dose” (air-equivalent detector) in aluminum and 
interpolated to 1.25 Mev. 

The column headed ¢, beside each column of 
calculated data contains percent deviations from the 
recent Berger-Spencer [15] results (f), 


B—B (Berger-Spencer 


) : 
i 100% 
B (Berger-Spencer) x10 \M% 


€s (13) 


and the column headed e, contains percent deviations 
from the Goldstein-Wilkins results (g) to which (b), 
(d), and (e) were fitted, 


B—B (Goldstein-Wilkins) os 
‘=~ B Goldstein-Wilkins) 1°? 


(14) 

As seen in table 1, the Goldstein [11] linear 
approximation (a), above, involves neither fitting 
nor transport theory calculations. It is very useful 
for penetrations up to one or two mean free paths, 


| and gives results within a factor of two up to ten 


mean free paths in aluminum. 


: atr-equivalent detector), Ber(Eo,Z), for 1.25-Mev gamma rays 


in aluminum evaluated from five different analytic approximations (a to e) at 0, 1, 2, 4, 7, and 10 mfp, compared with the 
corresponding moments-calculated analytic (f) and numerical (g) results 





| (a) Goldstein [11] l-parameter | (b) Chilton et al. [12] 2-pa- 


linear approximation (6) 


Bett Bett 


1.00 
2. 06 
3.01 
5.04 
8. 04 
11.0 36 | —37 


rameter exponential-linear fit 


(ec) Leshchinskii [13] 1-param- 
eter exponential-linear fit 
(Ya) 


Bett 


1.00 
1. 80 
2. 41 
—18.9 
— 1490 
— 52100 








(e) Taylor 

(d) Capo [9] 3-parameter poly- | 

nomial (cubic) fit (10) Lakey [28] data 
Bett 


Beri 


0. 997 

1. 93 

3.05 

5. 81 

11.1 

10 17.6 


[14] 3-parameter | (f) Berger-Spencer [15] 4-pa- 
re ye fit (11); 


(g) Goldstein-Wilkins [26] 
moments calculation re- 
sults; Al data in tabulated 
form only 


rameter exponential-linear 
moments calculation re- | 
sults (12) 








* e-=([Beri— Bpri(Berger-Spencer) }/ Bpr1( Berger-Spencer) } X 100%. 
** ¢-=([(Breri— Brri(Goldstein- Wilkins) ]/ Ber1(Goldstein-W ilkins) } X100%. 
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The two-parameter exponential-linear empiric ot 
fit (b) of Chilton, Holoviak, and Donovan [12] 1 
surprisingly good for penetre ations up to ten mean 
free paths. This formulation, similar in form to 
the calculated result (f) of Berger and Spencer 
[15], can thus provide economical storage and easy 
application of buildup data over a wider range of 
penetration thicknesses. Reference [12] contains 
parametric data for aluminum only, which was the 
deciding factor in selecting the reference 
substance for table 1 

Leshchinskii [13] supplied data, for the reference 
situation of Co-60 gammas in aluminum, only for 
the one-parameter exponential-linear fit (9a), in- 
tended for penetrations of 0 to 10 em (0 to 1.49 
mean free paths), although for Co-60 gammas in 
water he did use the three-parameter fit (9b) for 
penetrations of 2 to 70 em (0.127 to 4.44 mean 
free paths . In each case the fit is within about 5 
percent up to the specified upper limit and then 
departs wildly. However, this formulation does 
have the advantage that the exponential dependence 
exp Mol is canceled out of the distributed source 
integral 4), and it will be seen that the power-series 
solution of (4) [S, 10] is, in effect, a generalization 
of the Leshchinskii result. 

The slight departure of the Capo data (d) from 
the condition Bpy(, Z, 7 l at por—0 results 
from additional fitting as a polynomial in 
energy or reciprocal energy, in which it is difficult 
to preserve this normalization. This discrepancy 
is not large enough to detract from the usefulness 
of the data, however, and in this specific comparison 
Mey photons in Al) the departures ey and Ey 
calculated data (f) and (g) are seen to be 
over the given range of yor values, than cor- 
responding departures of (a), (b), (c), and (e). 


this as 


cross 


from 


less. 


In addition to the empirically fitted Capo tabula- 
tion, 8; polynomial coefficients for Co-60 and 
Cs-137 gamma rays in concrete and water wer 
directly calculated by Spencer [1] using the moments 
method and were punched on IBM cards for use 
in further calculations. In [1] the <A, coefficients 
in (Al), page SO, are identical with coefficients 
if (a) A; refers to the zeroth Legendre harmonic 
(J=0), and (b) A; refers to integral dose data rather 
than to the differential energy data indicated by 
the dependence on /; in (Al 

The sum-of-exponentials three-parameter fit (4) 
of Taylor [14] is perhaps the most well-known and 
widely used buildup factor formulation in current 
use |27 
aluminum seem to be those of Lakey [28], and these 
were used for the comparison in table 1 

The exponential-linear four-parameter formula (f) 
of Berger and Spencer [15] represents direct calcu- 
lated results by the moment method [21]. Parametric 
data in [15] are tabulated for 0.0341- to 10.22-Mev 
gamma-rays in aluminum and concrete, and addi- 
tional formulas are supplied for applying these build- 
up parameters to such geometries as the isotropic 
disk source (on-axis) and the 


isotropic spherical 


The only published parametric data for 


surface source. The buildup data in table 1, obtained 
by evaluating (12) using parameters interpolated to 

1.25 Mev, agree within 4 percent with the earlier 
moments calculation [21] results of Goldstein and 
Wilkins [26], (g), to which the formulas (b), (d), and 
(e) were fitted. 


4. Power-Series Solution for Unscattered 
Radiation From a Distributed Source 


It has been shown [6, 7] that for finite plane sources 
of simple shape, a solution of (4) for )” (i.e., omitting 
the factor Bpr, (4, Z, r) within the integrand) may be 
obtained in the form of an infinite series 


IP=(e 4r) >) ¢ (veom) - (upa) (15) 


n=0 


which converges quickly when the detector is less 
than a mean free path from the source. The dimen- 
sionless g,(geom) coefficients in (15) depend only on 
the shape of the detector-source geometry, and are 
derived from a power expansion of the 
exponential in (4) 


series 


exp ( 
such that 


r/x)"dS/(r? -n!'). (17) 


q; (geom ) 


The quantity yo, also dimensionless, is the mean- 
free-path thickness of material between source and 
detector. 

Formulas and tables of ¢,,(geom) coefficients (tables 
include values for 0<n<9) have been generated for 
a rectangular source |7] and for a circular disk source 
with detector off-axis [8]. As a check, it was noted 
[8] that for a detector on-axis over a circular disk 
source the series (15) reduces to the familiar exponen- 
tial-integral solution, in which the exponential 
integrals /,(#) are replaced by power series in z. 


5. Adaptation of Analytic B).(E,, Z, r) For- 
mulations to the Distributed Source Power- 
Series Solution 


In the series solution (15) for the unscattered 
component J” in section 4, the material between 
source and detector can either be distributed through- 
out the intervening space, or can be arranged as one 
or more uniform slabs or sheets lying parallel to the 
source-plane. In this section, on the other hand, 
the source and detector are assumed embedded, 
discussed section 2, in an infinite homogeneous 
medium. 
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Replacing the polynomial in the Leshchinskii for- | 
mula (9b) by & power series we obtain the buildup 
factor formulation 


>>, @ 


n= 


Bor (Eo, Z, r) exp (+ Uo? (Ey, Z) - ( Mor)" n! 


(18) 


reported in [10], from which the coefficients * 5, (2, Z) 
can be used with the series solution (15). The only 
change in (15) is the multiplication of each series 
term q,(geom)-(uor)" by the additional factor 
b, (2, Z) to give the total dose-rate solution 


D=ITP , Z) + dn(geom) - (uor)” 


(19) 


(o/4r) >) 5, (Ey 


n= 


Equations (18 and (19) are, of course, little 
practical value if an infinite number of b,(4, Z) co- 
efficients must be fitted and tabulated. However, by 
expanding as power-series all expone — appearing 
in the analytic formulations (a) to (f) i 
then equating coefficients of like mail of por, a 
transformation for generating 6,(/, Z) coefficients 
from the parameters of any of these formulations may 
be obtained. Table 2 presents these “parasite’’ trans- 
formations, and includes (18) for comparison. 


section 3, 


‘The simpler notatior 
8’, (Ko, Z 


work, is equivalent to 


used it 


TABLE 2. 


in which i 
nerating the b, 
the 


The power-series 
a tof) Is use 
was used ti 


formulation, 
ad for ge 


» generate 


sample 6,(F A 


Reference Number of 


parameters 


Beti( Eo, Z, 1 


Goldstein [11] 1 | 1+ay,r 


{Chilton, Holoviak, ° : 
\ Donovan [12] 2 | I+paner exp 


exp (+,,r) (6r+1) 
Leshchinskii [13] 


Summary of analytic formulations (a to f) for the point isotropic source buildup factor Bpri(Fy, Z) 


” for number of parameters (column 3 
Zz coefficients, has been added for comparison 


, point isotropic source 


— bor) 


The surprisingty simple transformation (e) of 
Taylor [14] and Lakey [28] parametric data results 
because inclusion of a Bp rif), Z,7) of this form into 
(4) produces no change in the analytic form of the 
integrand. It can also be seen in table 2 that the 

Goldstein {11] linear approximation (a) is a special 
case of the polynomial fit (d) used by Capo [9], and 
that the four-parameter calculated result (f) of 
Berger and Spencer [15] is identical in analytic form 
with the two-parameter empirical fit (b) of Chilton 
et al. [12]. The similarities are reflected in the 
respective transformations in the last column. 

To illustrate the use of table 2, transformation (d) 
was used to generate the 6,,(/), Z) values for 0<n<13 
in table 3 from “dose” (air-equivalent detector) 
buildup data for water, Al, Fe, Sn, Pb, W, and U in the 
form of 6,(/£y, Z) polynomial coefficients. The input 
B,(, Z)’s were taken from reference [9], tables 1B, 
4B, 7B, 10C, 13B, 16C, and 18B which in turn were 
fitted by least squares to data calculated by Gold- 
stein and Wilkins [26] using the moments method. 
Of the 15 gamma-ray energies from 0.4 to 9.5 Mev 
presented in [9| for reactor “operating” and “shut- 
down”’ shielding calculations, the 10 energies L,=0.5, 
0.7, 1.0, 1.5, 2.0, 2.5, 3.5, 5.5, 7.5, and 9.5 Mev were 
arbitrarily selected. 

The en ince in table 3 of b)(4), Z Bo( Ey, Z) 

B(Ey, Z,7) at por=0} values which are not unity 
follows foes the discussion of the Capo data in 
section 3. 


(column 4 


may be replac ed by 1, 2, 3, 4, or N, 


The last column contains the 


m which of the formulations 


depending « 
”’ transformations, of which (d 


** parasite 


in table 3 from the Capo [9] tables of 8; polynomial coefficients 


Eq No. 


exp (+ ,7r) (ar?+-br+e) 


N 

Capo [9] : N z Bi( Eo, 
i=0 

{ Taylor, [14] 

\Lakey [28] 


Berger, Spencer [15] 


Power This work 


series 


is 8] 
exp (+,7) z= 
n=0 


Z) (uor)* 


A exp (—dipr)+(1—A) exp (—aayor) 


1+por{ Ai B,; exp [—( Bi 
+ AB: exp [- 


(n—i)'\Bi( Eo, Z 


—A)(1+<a3)* 


1) por] 
Br—1)yor)} 


n 
1-> i(%) (4.B(Bi-1 i-t 
i=l 


+A:2B2(B.—1 


b»(Eo, Z 
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TABLE 3. *Coefficients b,(Eo, Z) of the power-series buildup factor formulation (18) derived from the Capo [9] 8; polynomial coe fii- 
a n ’ A } J J J y 
cients according to transformation (d), table 2. 
n (Eo, Z) coefficients may be used directly in the distributed source series solution (19) with geometry coefficeints g, (geom) (17) such as are given in refe: 
7] and {8} for rectangular or off-axis disk sources. 
bo( Eo, HO) 


These 


rences 


1. 0008 (0) 9. 9678 (- 9.9480 (—1) 

5. 0643 9. 7950 2 6752 (—2) 
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0506 (0 
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2 
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—5 
a 
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5215 (—1) 9. 0059 
2030 (0) 5. 6633 
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2242 (0) 0904 
3060 (0) 4128 
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0729 (—2) 5912 
3421 (0) 2502 
2. 8321 (0 
5870 (0) 
5691 (0 


9640 
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0190 
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R28 
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We te ett fet eet et et 
1 } 
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9.9849 (—1) 0008 (0) 
2042 (—1) 5.6586 (—1) 
3354 (—1) 3872 (—1) 3. 73 
6191 (—1) 8031 (—1) 1564 
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6359 (0) 0928 (0) 5622 

2. 0785 (0) 4856 (0) 6859 
4896 (0 8690 (0) 2741 
8676 (0) 2427 (0) 5724 
2110 (0) 2. 6069 (0) 8635 
1 2.9598 (0) 1471 
7882 (0) 3. 3026 (0) 2. 4230 
0189 (0) 3. 6343 (0) 6910 
2092 (0) 3.9548 (0) 9510 (0) 


0032 (0) 
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TABLE 3. *Coefficients b,(E, Z) of the power-series buildup factor formulation (18) derived from the Capo [9] 8; polynomial coe ffi- 
cients according to transformation (d), table 2.—Continued 
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*See footnote at end of table. 
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TABLE 3. *Coefficients b,(Eo, Z) of the power-series buildup factor formulation (18) derived from the Capo [9] 8; polynomial coe ffi- 
cients according to transformation (d), table 2.—Continued 
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TABLE 3. *Coefficients b,(Eo, Z) of the power-series buildup factor formulation (18) derived from the Capo [9] 8; polynomial coe ffi- 


cients according to transformation (d), table 2. 
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0.064656 


6. Utility and Limitations of the Power- 
Series Method 


With respect to use of the infinite series (19), 
gamma-ray distributed source problems fall into 
three classes: 

(a) In the region 0< x2 0.01, as is often the case 
in an air medium [29], all terms in (19) may be 
neglected except for the zeroth term (¢ 4m) qo(geom). 
This term is the complete solution for pure inverse 
square law attenuation as assumed by Meredith 
[29], and depends only on relative rather than abso- 
lute dimensions of the detector-source geometry. 

(b) In the region 0.01 SyorS1.0 the series (19) 
has its greatest usefulness, either with or without 
inclusion of scattered radiation by 6,(2, Z) coef- 
ficients. This will be demonstrated in the numerical 
examples in section 7. 
as is the case in 
fallout shelter calculations, the series (19) is usually 
not useful, and must other 
available methods [1, 2]. 


(c) In the region 1.0 uor< @, 


recourse be made to 


It is interesting to note that the speed of con- 
vergence of (19) is usually enhanced over that for 
(15 BR 


rather than the converse. A possible explana- 
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which the adjacent entry is to be multiplied; 


tion for this enhancement (with at a 
proof) is offered as follows: 

In an absorbing and scattering medium the 
point-source buildup factor Bpr:(, Z, 7) is bounded, 


roughly, by the limits 


no attempt 


1< Bers (Eo, Z, r)> exp(+ mor). (20) 
The lower bound, Bpr: (2, Z, r)=1, corresponding to 
pure absorption (D'=0 in eq (1)), would lead to 
buildup coefficients all equal to unity, ie., 
b.(. 2)=1, n=0, 1, 2, (21) 
in which (19) would be identical with (15). 
The upper approximate bound, Bpr(,Z, 1) 
exp(+-yor), corresponding to a pure scatterer (pho- 
tons may be scattered, but never absorbed or degrad- 
ed in energy) is the same as simple inverse-square-law 
attenuation, such that the point-source eq (2) would 
reduce to 


cause 


99) 


(22 


D=(k/4x)/r’. 
This buildup-factor behavior can be applied to (19) 
by setting all except the zeroth buildup coefficient 
equal to zero, 1.e., 





b( & . B) 
b,, (Eo, Z 0, 


Hence, in the intermediate region indicated in (20) 
we might expect the observed enhancement. 


7. Numerical Examples 


7.1. Example 1. Rectangular Co-60 Source, Strati- 
fied Barrier of Steel and Water. Comparison 
With Experiment 


There exist little experimental data in the moder- 
ate-penetration region 0.01 >yor>1.0, discussed in 
(b) of the preceding section, for simple geometries 
amenable to direct analysis using (19). However, 
in connection with food irradiator design studies [17] 
data have been published for the more complicated 
situation of a two-layer barrier. Here, Donovan 
[30] measured the dose-rate received by a detector 
separated from a rectangular Co—60 source by a 3- 
in. laver of water and varying thicknesses of steel 
cladding as shown in figure 3. These results are 
over the range 0.589<yor(water+steel)< 1.225, 
hence they provide an excellent check on (19) up to 
the upper limit of its usefulness as discussed in (b) 
of the preceding section 

The agreement between theory (infinite homo- 
geneous medium) and experiment (layered inhomo- 
geneous medium) in the following example must be 
considered somewhat fortuitous, although the posi- 
tion of the detector within the water phantom tends 
to minimize boundary effects. Since interest has 
been expressed in adapting eq (19) for use in nuclear 
engineering curricula, this example is presented in 
detail, particularly in the steps involving 
dimensional analysis. 


creat 


a. Situation 


A 6-in. thick water-equivalent slab phantom is 
separated by an airgap of 2.5 in. from a 50-in. wide 
by 57-in. long rectangular plaque Co-60 source. 
This source, with a uniform activity of 1 curie per 
em’ of area, is encased in steel cladding of uniform 
thickness 

b. Problem 


Find the dose-rate (rads [31] in H.O) in the center 
of the phantom as the thickness of steel cladding is 
varied from 0.097 in. to 0.679 in. Use the simplifying 
approximations: 

1) To take into account the finite thickness of 
the Co, Z=27), include half the source 
thickness in the thickness of steel (Fe, Z=26) repre- 
senting ‘‘cladding thickness” in the calculation, as 
was done in an alternative analytical study of this 
probli hn ; 

2) Assume that the Co-60 nuclides emit two 
1.25-Mev photons per disintegration, rather than a 

17 Mev and a 1.33 Mev photon as is actually the 


case. 


source 
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FiagurE 3. Donovan [30] food irradiator rectangular source 
geometry treated in example 1. 


(3) Kalos ({21], page 796) has shown that for a 
composite barrier consisting of a layer of high-Z 
material followed by a layer of low-Z material, the 
resulting buildup factor shows a smooth rise from the 
value for pure high-Z material toward that for pure 
low-Z material, the rate of approach depending on 
the source energy. Indeed, substitution of the linear 
approximation (6) for Bpy(,Z,r) in’ Kalos’ 
empirical eq [21], (26.4) gives the linear connecting 
form 


Beri(Eo, Z;, Z2, 11 +12) 


wo( Zr 
bo( Z1) 11+ Mo( Zo) 2 


Uo(Z2) i 
o(Z PM! L )I's 


Bor, (Eo, Z1, 71) 


Bpr(E%o, Z2, 172) (24) 
in which Z,, 7, refer to one layer of the barrier and 
Z>, r. refer to the other layer of the barrier. Since 
the buildup for iron and for water are not greatly 
different [32] over the penetration range in this 
problem, the approximate combinatorial form (24) 
will be used. 


c. Plan of Attack 


Evaluate the dose rate (rads in H,O/hr) using (19) 
assuming the medium (a) to be all water: 


D(H,Q) = (0/42) 546, (1.25 Mev, H,O) - ¢,(geom) - (uoxr)” 
| (19a) 


and (b) to be all steel: 


dD Fe) 


(o 4x) >“, (1.25 Mev, Fe) 


* Jn(geom) - (ugar )” 


(19b) 


in which (19a) and (19b) differ only in the values of 
the buildup factor coefficients 6,(/, Z). Then use 
(24) to obtain a dose rate between these upper and 
lower bounds by weighting the results from (19a) 
and (19b) by the fractions of mean-free-path distance 
traversed by the primary photons in each material: 


D (Mo) 11,0/ (Mo) total D(H.0) 


(492) re! (Mo®) totar} D( Fe). (25) 





d. Source Strength Constant (o/47) 


For an answer from (19) in units of rads/hr per 
curie/em’ of source strength, the factor (¢/47) is 


(0/4) = 1.7005 « 10° « N(photons of energy /,/disin- 


tegration) « &, (Mev/photon) Men( 4, Z)/p} (em?/g) 


(26) 


in which the numerical constant 1.700510 is 
obtained from the dimensional factors 


1.7005 « 10° = (1/47) (steradians~!) 3.71 * 10” 
(disintegrations/sec)/curie} * 1.6 1078 {rads 
(27) 


(Mev/g)} *3.6 10° (sec/hr) 


For the problem stated above the numerical values 
of the other constants in (26) are 


N=2 photons/disintegration 
ky =1.25 Mev/photon, and 
Men(1.25 Mev, H,O)/p=0.0296 cm?/¢5 


Substitution of these values in (26) gives a source 
strength constant (¢/47) of 


1.258 « 10' 


(o/47) 


for values of D(H.O) and D(Fe) from (19a) and (19b) 


in units of {rads (H.O)/hr 


curies/em? 
e. Buildup Factor Coefficients },, (Ey, Z) 


Since reference [9] and hence table 3 of this work 
do not contain data at 1.25 Mev, some interpolation 
is required to obtain the b,,(£, Z) series coefficients 
needed in (19a) and (19b). This interpolation is 
conveniently performed on the “untransformed” 
8(E, Z) polynomial coefficients, from which any 
number of 6,(/, Z) series coefficients may be gener- 
uted using (d), table 2. Table 4 contains poly- 
nomial coefficients: 

(a) 811.25 Mev, H.O) 
table 2, and 

(b) 8,(1.25 Mev, Fe) interpolated from [9], table 
7B 


interpolated from [25],° 


5 Interpolated from [24], table 
The coefficients in [25] as in 


9] were fitted to the calculated results in NYO 
3075[26] 


TABLE 4. Polynomial buildup factor coe flic tents BE 
used as input data in examples 1 and 2 


Z) 


0 
’ 


8;(1.25 Mev, H20)* 8;(1.25 Mev, Fe)** 


1.00 (0 1.013 (0 

9.00 l 7. 455 1 
1.40 1 6. O55 2 
3. 90 3 5.320 (—4) 


*Interpolated from reference [25], table 2. 
**Interpolated from reference [9], table 7B. 


Corresponding values of 6,(1.25 Mev, H,O) and 
6,(1.25 Mev, Fe), computed using the numbers in 
table 4 as input values for formula (d), table 2, are 
given in table 5 for 0<n<30 and are ready for use 
in (19a) and (19b). 


TaBLeE 5. Buildup factor and geometry data used in example 


1, serves solutions (19a) and (19b). 


The b,( Eo, Z) coefficients, of which 6,(1.25 Mev, H2O) were also used in example 
2, eq (19¢), were derived from table 4 using transformation (d), table 2. The 
qn(geom) geometry coefficients were obtained froin [7], table 3 (and unpublished 
higher coefficients from formulas (57) and (58), as 


qn( Donovan geom) =4q,(a,b 


using a=(w/2)/h=5, b=(l/2)/h=5 and the factor of four to apply the detector- 
over-corner data to Donovan’s centered-detector geometry For h=5.5 in 
these a and > values imply plaque dimensions of 55 by 55 in., 
Donovan's 50 by 57 in. 
by less than 1 percent 


’ rather than 
However, the solid angles subtended in each case differ 


bn(1.25 Mev, HO) | 6,(1.25 Mev, Fe) 


a(geom 


f. Geometry Coefficients ¢,(geom 


The parameters describing Donovan’s 


experi- 
mental geometry, as shown in figure 3 are: 


w=plaque width=50 in. 


plaque length =57 in. 
distance (height) of detector above plaque center 
2.5 in. airgap+3 in. water 


5.5 in. 


In the notation of [7] the geometry coefficients 
for this centered-detector situation could be cal- 
culated exactly as 


qn(Donovan geometry) =4q,(4, 5) 





where dn(a, 6) are corner-position rectangular-source 
geometry coefficients and 
(w/2 h 4.5454 


(50/2 


L/2)/h 57/2)/5.5=5. (29) 
However, g,(a. } 
available l \7 
form 


values for a=5, 

for O0O<n<9 and in unpublished 
We can use these tabulated 
coefficients and avoid a two-way interpolation, or an 
exact calculation of Gn\2eom coefficients 
formulas in [7], by the following argument: 

It has been noted [1-3] that veometry effects scale 
approximately as the solid angle, 2,7 subtended by 
the source from the detector for barrier thicknesses 
of the order of one mean-free-path. 


5 are already 


ior O<n<35. 


from 


Hence, since 
the solid angle corresponding to the already tabu- 
lated coefficients differs from the solid angle in 
Donovan’s geometry by less than a percent, i.e., 


1.5454. b 5.1S18) 


sZ( a 


19(4.5454, 5.1818 


5.170 steradians/5.132 steradians 


1.007 (30) 
the available coefficients for a 5, 6=5 
according to (19) and are given 
of table 5 for 0<n<30. 


are used 
as the last column 


g. Penetration Thickness yur 


The last factor. 


uot)”, in each series term of (19a 


and (19b) is for this problem simply the n’th power 
of the barrier thickness in mean 


hick free paths (abbrevi- 
mfp If this thickness is given in 
hickness in mfp is 


ated 


inches, 
the t 


por(mfp)=p( Ey, Z)(em?/e) x p(g/em 


2.54 (em/in. 
where mili. 2, is 
coefficient [23 


the medium. 


the narrow-beam at tenuation 


as used in (2) and p is the density of 


For 1.25 Mev gammas in water we see that 


uox(mfp)/x(in.) =0.0634 (em2/e 


1.00 


0.161 mfp/(in. 


o¢/em’) 2.54 (em/in.) 


of water) (32) 


and similarly for iron that 


wor (mfp)/x(in.) =0.0533 (em?2/e 


7.83 (g/em*) & 2.54 (em/in.) 


1.06 mfp/(in. of steel). 33) 


ntical with the tabulated 


quantity 7 


P,f 


| Thus, for example, 


3 in. water=3 in. <0.161 mfp/in.=0.483 mfp, 


and 


0.097 in. steel=0.097 in. 1.06 mfp/in.=0.103 mfp. 


h. Numerical Results, D® and D 


Donovan’s measured values of D for 3 in. of water 
preceded by various thicknesses of steel [80] are 
shown in figure 4, with his estimated +3 percent 
uncertainty indicated [33]. 

Before evaluating the total dose D according to 
(19a), (19b), and (25), it is of some interest to first 
evaluate the unscattered component 1° according 
to (15). We ean then note the effect of buildup as 
characterized by the 6,(/, Z) coefficients. In- 
serting the data from the preceding sections (d), 
(f), and (g) in (15), sums for a variety of yor values 
yield D® (curve (a) in fig. 4) which is independent 
of the kind of medium and the photon energy 
except for the source strength factor ¢/4%. Indeed, 
except for the factor o/4m the data for curve (a) are 
already tabulated in [7], table 4. 

Next we can evaluate (HO) and D(Fe) from 
(19a) and (19b) using the same set of summations 
as for D® (curve (a), above) except that each series 
term multiplied by the appropriate buildup 
factor coefficient 6,(/2),Z) from table 5. These 
sums D(H.O) and D(Fe) (for Donovan’s geometry 
embedded in infinite media of water and _ steel, 
respectively ) are displayed as curves (b) and (c) 
in figure 4. The circled numbers are a qualitative 
indication of the convergence rate of (19) and are, 
for the indicated calculated points, the number of 
the series term (n+1, including the zeroth term) 
beyond which the partial sum remains within 0.1 
percent of the sum up to n=30. 

Weighting )(H.O) and D(Fe) (curves (b) and (c)) 
according to (25) we obtain J), the hybrid curve 
(d), our desired result. 

The calculated values of )) in curve (d) can now 
be compared absolutely with Donovan’s measured 
values, since the method is free of any arbitrary 
normalization. Agreement seen to be within 
Donovan’s 3 percent estimated uncertainty. 

The results of a more approximate but simpler 
analysis by Moote [3] are shown as curve (e). In 
this analysis, which gives an excellent fit to the data 
in the neighborhood of wrv=—1, the dose rate D is 
given as the simple product of: 

(a) The source strength constant o/47, 

(b) the solid angle Q, indicated in figure 3, sub- 

tended by the source from the detector, 

(c) the exponential integral /;(uor) for the total 

barrier thickness pox, and 

(d) the plane isotropic source buildup factor, 

Bozy(,H20,uor), for the dose rate at a 
distance yor from an infinite plane source 
in water. 


is 


is 





EXPERIMENT [30] 


Dp? 


DOSE -RATE, (RADS/ hr)/(CURIES /cm2) 


[ (N)- NUMBER OF 
TERMS TO SUM 
FOR CONVERGENCE 
TO €01% FROM 
SUM TO n= 30 


0.5 


DiFe) 


oo 
= (0/4) X Ingeom) “(0x)” (7) 
n= 


( 
D(H,0)? (9/42 b,,(Eo,H.0) 4,(9e0m)-(wox)” 


_ 30 n 
= (a0 /4m) 2 Paleo, Fe): q,(geom):(yzox) 
n= 


~ ‘He ™)HoO D(H20) + (po x)FeDiFe) 
(Ho¥)H0 + (woX)Fe 





* (0/4) X(geom): E(wox) Bo, (EoH,0,u0x) [3] 


4 


C080 


(2 125-Mev PHOTONS ASSUMED) 


10 


Mok = (Ho HO + (Ho X)Fe 


Figure 4. Results for the 


Donovan's experimental points 


The dose-rat« 


4), 33] are compared with the theoretic 
component D® due to primaries (eq 15 


Donovan food irradiator rectangular source geometry (example 1). 


il results 


b) The total dose-rate D( HO) assuming an all-water medium (eq 19a). 


The total dose-rate D(F¢ 
1) A weighted average, 
e) The total dose-rate, 


issurming 


D, of curves (b) and 


7.2. Example 2. Dose-Rate Profiles Across a Cleared 
ee Area in a Co-60 Infinite Plane Source 
in ir 


Along the aris of a circular disk source, analytic 
solutions for the flux in terms of exponential integrals 
and the buildup formulations [9, 11-15] are well 
known. However, off-aris solutions including build- 
up data have not, to the author’s knowledge, yet ap- 
peared in the literature. The power-series solution 
(19) does apply to this situation, as will be illus- 
trated in this example. 


a. Situation 


A circular area 600 ft in diameter has been cleared | 
in an infinite plane isotropic source of 1 curie/ft? of | 
Co-60. The surrounding medium is air (no ground- | 
air interface is taken into account). 


an all-iron medium 


eq 19b). 


swccording to (2 
D, according to the solid-angle exponential-integral approximation used by Moote [3] 


}- 


available in parametric form (9, 11-15], but water 
buildup data can be used to good approximation 
(1. 2) 
ji, 4]. 

c. Plan of Attack 


Evaluate the dose-rates (rads/hr in air) at h 
30, and 300 ft above the infinite plane source [25, 
34] using 8,;(Z, Z) coefficients from the polynomial 
buildup formula (10) according to 


9 
», 


Doz lh )=(0/4er ) 27 E, (ugh ) 


rexp (— Loh) By t Bo T 283 


+ (B2+28s) woh t+ Bs(uoh)*}] (34) 
in which E,(yh) is the ordinary exponential integral 
[35] and values 6,(1.25 Mev, H.2O) are already given 
in table 4. Then, using the derived values’6,(1.25 
Mev, H,O) from table 5, evaluate dose-rate profiles 


| across a 600-ft diam disk source with the same source 


b. Problem 


Find dose-rate profiles, relative to the undisturbed 
infinite plane source, at heights 3 ft (0.00626 mfp), 
30 ft (0.0626 mfp) and 300 ft (0.626 mfp) above the | 
plane of the source. Buildup data for air are not | 
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| Dorsx(e, h) 


| strength constant (¢/47) according to (19) 


(o/4m)>5 b, (1.25 Mev, H,0) 


n=0 


 Qn(P, h) + (uph)" (19e) 





where g,(p,h) geometry coefficients are available in 
[8] for 0<p<10 disk radii and 0.1<A<10 disk radii, 
and in unpublished form for 0.01<A<0.1. In this 
problem, A assumes the values 0.01, 0.1, and 1.0 disk 
radit while the penetration thickness [oh assumes the 
corresponding values 0.00626, 0.0626, and 0.626 mfp 
given in the problem. 


d. Numerical Results 


Evaluation of (34) using the source strength con- 
stant o/44—12.24 (rads in air/hr)/(curies/ft?)* gives 
infinite-pl: ine dose-rates Dora(h t) of 423.7, 247.6, and 
78.83 (rads in air/hr)/(curies/ft?) for the he ights 3, 
30, and 300 ft, respectively, as indicated by the hori- 
zontal lines in figure Summation of (19¢) for these 
three heights and a number of p-values gives the 
dose-rate profiles Dprsk(p, h) also shown in figure 5. 
The circled numbers, as in example 1, indicate the 
required number of series terms for convergence to 
within 0.1 percent of the stable value. 

Subtraction of Dorsx(e, h) from Dpzr(h) gives the 
desired profiles shown in figure 6, normalized to 100 
percent for each infinite plane source result Dpz1(h) 
as specified in the problem. 


§ This value of o/4x is equivalent to the irke and Buchanan ((34], table}1) 
value of 14.0 (roentgen i 


CGo60 
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AIR MEDIUM, H,O BUILDUP DATA 
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Figure 6. Relative dose-rate profiles at heights 3, 30, and 300 ft 
above a 300-ft radius cleared circular area in an infinite plane 
Co-60 source in an infinite air medium (example 2). 


These profiles were obtained by subtracting the disk-source profiles Dprsx(p,h) 
in figure 5 from the corresponding infinite plane source dose-rates Dpz(h), 
dividing by Dre1i(h) and multiplying by 100 percent 


Caution must be taken in applying these results to 
the study of radiation fields from fallout deposition 
on the ground, since the ground-air interface has 
been shown by theoretical [36] and experimental [37, 
38] studies to affect the scattered component by as 
much as a factor of three. Qualitatively, this would 
have the greatest effect on the h=3 ft profile in 
figure 5, “‘squaring the shoulders” of the curve and 
depressing the tails on either side. In figure 6 the 
3-ft profile would more closely resemble a ‘square 
well,” but the 30-ft and 300-ft profiles would not be 
appreciably changed. 


8. Discussion 


The power-series formulation for buildup data is 
seen to have two main limitations: 

(a) Convergence of the primary-flux series (15), 
and hence the total-flux series (19), is poor for 
combinations of thick barriers (ur = 1) and detector- 
source separations “e h less than the source dimen- 
sions (i.e., a, b or h< <1). In such cases the 
geometry bth ly sana the infinite plane source 
geometry for which the series (15) for D° becomes, 
in the limit, the series representation ((5), [35]) for 
the exponential integral /\(uor). The latter series 
converges absolutely for O<uoar<o, but for 1< 


’ 
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uor<_@ the convergence rate is rather slow for prac- 
tical use. In these situations it is more economical 
to revert to numerical integration, or to use an 
infinite plane source analytical formulation, such as 
(34) in example 2, as a close approximation. 

(b) Boundary effects such as the density-interface 
effect [36-38] are neglected, but this limitation is 
inherent in any buildup factor formulation. Prob- 
lems in which such effects are important are best 
treated by Monte Carlo techniques at present [36]. 

On the other hand, the series formulation provides 
some advantages: 

(a) The three-way separation variables [(1) unique 
properties of the medium—®#,(/),Z), (2) geometry 
dn (geom) and (3) barrier thickness—por| gives both 
economy of data tabulation and flexibility of ap- 
plication. A given set of buildup data can be applied 
to a variety of distributed source geometries, or 
vice versa, and at the same time the barrier thickness 
fox can be varied in a trivial fashion. 

(b) In some simple idealized geometries the series 
(19) can provide exact answers, especially in the 
region of small barrier thickness (ur<1). These 
answers, in addition to available experimental data, 
can then be used to check the validity of simpler 
engineering-type approximate methods (e.g., [3]) in 
this region. 


The author thanks R. L. Bach for computing 


table 3, R. J. Herbold for computing the additional 


qn(geom) coefficients needed for example 2, and 
M. J. Berger for his interest and helpful criticisms 
regarding this effort. 
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An Alinement Interferometer 
J. B. Saunders 


(July 


Or 


“0, 


1963) 


This paper describes a diverging-beam type of alinement interferometer that permits 
the use of large apertures, provides ruggedness, high sensitivity, and is relatively compact. 


Since this interferometer is fully compensated, white light can be used. 
The narrowness of the zero order fringe, relative to separation of fringes, in such a well 


compensated system permits settings to better than one twentieth of a fringe width. 


This 


corresponds to a lateral displacement of the light source of less than 1.5 mm at a distance 
of one mile for a one-inch aperture instrument and to less than 0.15 mm for a 10-inch aperture 


instrument. 


The increase in dmands for higher accuracy in metrology includes that of pointing. 
The alinement interferometer described here is a precision device that is rugged and relatively 


easy to apply. 


Its performance,because of its high sensitivity, is severally affected by the 
homogeneity of the atmosphere or medium through which the light passes. 


However, 


this effect is also greatly reduced by the proximity of any two interfering rays of light. 


There appears to be a growing need for an accurate 
method of alining objects over distances of several 
miles. The diffraction methods, described by Van 
Heel! and Dy son-, appear to be inadequate for such 
distances. Other methods such as the axicon® and 
the telescope also have limited accuracy. The aline- 
ment interferometer described here is believed to be 
practical and more sensitive than previously de- 
scribed methods. 

An arrangement of optical elements that form a 
precision alinement interferometer is shown in figure 
1. It consists of two double-image prisms‘, P; and 
P,, a beam divider, B, a lens, LZ, a source, S, and two 
pairs of quarter-wave plates, C; and C,. Figure 2A 
shows the course of the light from B to the observer 
located either at /’, or F’’, for measurements in the 
plane of figure 1, and figure 2B shows the same for 
measurement in the plane perpendicular to that of 
figure 1. 

The alinement problem may consist of determining 
the position of the source, S, located a distance, d, 
from the optic axis of L, the precise adjustment of S 
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ion of a Késters double-image prism; J. Res. NBS 


DIVIDING 
PLANE 


FIGuRE 1. 


The dividing plane of the double-image prism is a smi-reflecting optical film, 


on the axis of Z, or the alinement of several sources 
along the axis of L. Obviously, if the problem is to 
determine the position of S relative to one plane 
only (the plane of fig. 1, for instance), then P; and 
B may be eliminated. 

The optical sensitivity of this interferometer is 
proportional to the aperture. If the aperture of the 
prism is sufficient, the lens may be eliminated and the 
entrance face of the prism may be plane or spherical. 
A plane face prism, used without the aperture-in- 
creasing lens L, requires a collector lens as shown in 
figure 3. The 2d prism and beam divider are not 
shown in figure 3. They are, however, necessary 
for alinement in the plane perpendicular to this 
figure. 

This prism is difficult to make in large sizes. Small 
prisms are relatively easy to construct and may be 
used with a lens to increase the aperture of the inter- 
ferometer. Aperture sizes are limited only by the 
difficulty of producing large objectives (lenses or 
mirrors). A mirror objective is as easily used as a 
refractor and is free from chromatic aberration. 

The equation of sensitivity is obtained from figure 
4,in which the aperture is represented by the aperture 
of alens, L. Consider two rays of light that emerge 

| from the lens at points A, and A, that are equally 


Plane view of optical elements. 


It serves as the beam combiner in this 


interferometer, 


307 





FIGURE 2, ire 2A represents the path of the light in the plane of figure 1. Figure 2B repre- 
sents the same in a plane perpendicular to the plane of figure 1. 


tersection of these two planes represent the straight line of interest 


distant (y and —y) from Ao, which lie in the extension 
of the dividing plane of prism P,;. The distance Y 
(approximating the distance from lens to source), 
shown in figure 4, is measured from A». The plane 
of figure 4 is assumed to include the optic axis of L. 
Let p, and p, equal the optical paths from the source 
to the observer (at Fj or F%, fig. 2) through A; and 
Ay, respectively. The optical path from <A, to the 
observer is identical to that from A», because of sym- 
metry in the prism. The observer sees images of A, 
and A, superimposed, as shown in figure 5. The 
observed order of interference, N, at this pair of 
points, multiplied by the wavelength of the light, X, 
is equal to (p;—p.). Since XY is assumed quite large 
relative to y and d, we can ignore the refractive index 
of the lens and obtain, 





Uplics of an alinement interferometer. 
(1) 


p3 +(y+d)’. (2) 


On subtracting eq (1) from (2), and expanding, we 
readily obtain 


P2— Pi=(P2— Pi) (p2+ p1) =4 yd. 


On putting (p.+,)=2.Y (which it approximates for 
large values of Y) and solving for d we obtain 


2X- Nav 
4y 


d 


If A equals the aperture of the interferometer and 
the reference point is chosen at or near the margin 
then y='5A and 

4. Rays of light. p, and p2, are combined in the "AT 
Consequently, points A, and A» appear to coincide d X Nx 


(3) 
eference povnts, indicated in figures 5 and 6. A 





An A 
AN 


AW Ali, Zi 


Figure 5. The appearance of 
point source is: A, above the 
below the line of interest, 


interference fringes when the 
line; B, on the line; and C, 


If the problem is to adjust S on the axis, and the 
prism is adjusted to place the zero order of inter- 
ference on the reference point where d=0, then the 
adjustment consists in moving S perpendicular to the 
axis until N equal zero. If white light is used with a 
reflecting objective (a parabolic mirror), the zero 
order of interference appears to have excellent con- 
trast and an experienced observer can adjust the 
fringes to an accuracy of less than 0.05 fringe spacing. 
The corresponding error in N is 0.05 and for \=5.0 

10~° em (the approximate effective value for white 
light sources), the corresponding error in d is X/A 

<2.5X107° em. This corresponds to an error in 


 & 
d/X (the angular position of S) of= 10-° (A ex- 


A 
pressed in cm). If A equals 6 in. and X equals 1 
mile, the corresponding error in d is approximately 
0.3 mm. 

Two photomultiplier tubes placed in a differential 
arrangement can detect changes in the position of S 
to much less than this. Obviously, this principle 
can be used to guide an astronomical telescope when 
mounted on a satellite (above the atmosphere of the 
earth) to very high precision. 

Visual observations are most easily made with 
a prism that has a “‘built-in-wedge”’ of convenient 
magnitude to produce interference fringes whose 
widths are most appropriate for precision of readings. 
Prisms may easily be adjusted to provide fringes 
that are neither too broad nor too narrow. When 
photoelectric recording is to be used for null settings 
the prism should be adjusted to have no built-in 
wedge. This permits the use of the entire beam 
for detection. 

Figures 5A, 5B, and 5C 
appearance of the fringes in a built-in wedge inter- 
ferometer of 4 in. aperture when S is one second 
of are above, on, and one second below, the dividing 
plane of the prism through which observation is 
made. 

Because of differential polarization of light in 
the prism, unpolarized light cannot be used without 
properly oriented quarter-wave plates or a polarizer. 
Use of polarized light results in a loss of one-half 
the available light. The polarizer is used by placing 
it either between the prism and observer or near the 
source. If a polarizer is adjusted to transmit 
light vibrating parallel to the dividing plane of 
the prism, and if the observed fringe pattern is 
represented by figure 6A, then the pattern will 
change to that of figure 6B when the polarizer is 
rotated 90° (light vibrating perpendicular to the 
dividing plane of the prism). Therefore, unpolarized | 


respectively show the 


COINCIDENT POINTS, A, AND A, 
’ 
(ll 
A 
Figure 6. The appearance of interference fringes when: A, 
the light is polarized parallel to the dividing plane of figure 1; 


B, the light is polarized perpendicular to the plane of figure 1; 
and C, when quarter-wave plates are used instead of a polarizer. 





nm 








DOUBLE IMAGE PRISM 

















FIGURE 7. Arrangement of quarter-wave good 


plate s lo uve 
contrast with unpolarized light. 


light produces two sets of fringes with the maximum 
intensity areas of one set coinciding approximately, 
with the minimum intensity areas of the other set, 
thus eliminating contrast so that no fringes can 
be observed. 

Dr. Charles J. Koester > suggested and demon- 
strated that by cutting a quarter-wave sheet ° of 
mica at 45° to its optic axis and mounting it in 
front of the prism, in the manner shown in figure 7, 
the relative phase changes in the two polarized 
beams are practically equalized. The two sets of 
fringes produced by the two component beams 
(vibrating parallel and perpendicular to the dividing 
plane of the prism) will then coincide in phase. 
Good contrast is thereby obtained with little loss 
of light. The resultant fringe pattern (fig. 6C) 
indicates that the two sets of fringes, shown in 
figures 6A and 6B, are shifted in opposite directions 
by an amount that places them in coincidence. 

Tests made in the laboratory showed that the 
sensitivity of the instrument, obtained from eq (3), 
agreed within the limits of observations. The 
theoretical accuracy claimed above requires that 
the medium through which the light travels be 
quite homogeneous. An evacuated light tube would 
probably be required to obtain the accuracy indicated 
by eq (3) when measurements are made over long 
distances. 


5 Working at NBS under a National Research Council Postdoctoral Research 
Associateship and now with American Optical Co., Southbridge, Mass. 
6 Subsequent tests showed that the relative retardation is not critical for obtain- 
ing good contrast. 
7 


(Paper 67C4-141) 
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Determination of Optical Path Difference for 
Photographic Objective 


Francis E. Washer and Walter R. Darling 


(June 6, 1963) 


\ method of measuring longitudinal spherical aberration and optical path differences 


bv a direct visual means is described. 


and a movable-slit system. 


The method employs a nodal slide optical bench 
The underlying theory of the method is presented together 
with a brief description of the apparatus used. 


The results of measurement on a typical 


lens and a procedure for checking the consistency of measurements are included. 


1. Introduction 


One of the criteria used in evaluating the potential 


image forming qualities of a photographic objective | 
is the magnitude of the longitudinal spherical aber- | 
ration which measures the variation in axial focus 


with radius of annular zone of the lens. It is some- 
times desirable however to use a related quantity 
more meaningful in terms of physical optics, namely 
the optical path difference existing between the par- 
axial region of the lens and the various annular zones 
located between the axial and marginal regions of the 
The optical path diff rence expressed in units 
of wavelength yields information on the manner in 


lens. 


which light converging from a given annular zone of | 


the lens to a specified focal point may be expected to 
contribute to image formation. 
performance of a high quality telescope objective, a 
common requirement is that the optical path differ- 
ence over the entire wave front emergent from a lens 
shall not depart by more than 0.25 wavelength from 
a true sphere: 

A method has therefore been developed for the 


measurement of the departure of the emergent wave | 


front from a true spherical surface in units of wave- 
length. The method is essentially a refinement of 
the method reported by Chalmers! and similar to 
the two slit method used by Viisili.? The present 
method is well adapted for use with lenses having a 
wide range of focal lengths. 

The accuracy of the final results compares favor- 
ably with that of results obtained by more elaborate 
interferometric methods. Moreover the results are 
in such a form as to simplify the computation of 
changes in optical path difference with shift of the 
focal plane. 


2. Theory of the Method 


2.1. Determination of Optical Path Difference, P 


When a lens is affected by longitudinal spherical 


S. 1D. Chalmers: Proceedings of the Optical Convention 2, 156 (1912 
Y. Vilisiili, Ann. Univ enn. Abo. 1 3 (1922 


696—713—63 


In specifying the | 
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aberration, the position of best axial focus changes 
from zone to zone of the lens. The magnitude of 
this focal shift can be readily measured by a variety 
of methods, one of which is reported in an earlier 
paper.’ The effect of longitudinal spherical aberra- 
tion may a distortion of the 
wave front from a true spherical form by action of 
the lens. This distortion of the wave front 
variation in the optical path separating 
different portions of the wave front emerging from 
the lens and the focal point. It is possible to meas- 
ure this variation using a technique similar to that 
reported by Chalmers footnote 1). In this 
method, two parallel slits, separated by distance h 
are placed in front of a lens, the image of an infinitely 
distant bright line object formed by the lens then 
appears as a series of parallel interference fringes in 
the focal plane of the lens. The central fringe is 
the brightest of the group, the others growing fainter 
with increasing order of interference. 

The separation of these fringes is given by the 
relation 


also be described as 


pro- 
duces 


see 


rv] 


S 
h 


where S is the measured separation of the central 
fringe and an adjacent fringe, No is the slit separation, 
f the equivalent focal length of the lens, and \ the 
wavelength of the incident light. 

If the slits are so placed that one is on the axis of 
the lens, the second slit will be at distance hy from the 
axis. A difference in focus between the paraxial 
zone and the annular zone of radius fo, will shitt 
the image of the central fringe laterally by amount 
6, and the change in optical path, AP, is given by 
the relation 


F. E. Washer 
longitudinal spherical aber 





a radius of the 
a manner that the first slit occupies the 
position formerly occupied by the second slit, a new 
value 6, will be found and the change in optical 
path difference AP, is given by the relation AP 


This 


On moving the slits outward along 
lens in such 


process is continued until the last setting 
with the second slit near the margin of the 
and the corresponding value of AP Is AP. 

The total optical path difference P between 
and : Y ‘hh 


1, 2,3, 


centel hone designated by zone 


is given by the expression 
P:=ZAP 


expressed in 
\ elving the 
. which is 


units of wavelength by 
il path difference in 


opt 


waveleneth 
(4 


the aid which 
can be determined 
mt with respect to any 
ther point of the wave front located along the same 


of eq 4), the variation in 7 
Leh be interpreted as phase shift 
any point of the wave f1 


l 


liameter and for the particular image plane in which 


2.2. Effect of Focal Shift 


The relation, shown in eq j 
single position of the focal plane 
separated yy distance Af fro 
a different set of values of 


measured values of 


is valid for any 
For a focal plane 
the initial position, 
will be found. The 
by amount Ad which 
related to Af as shown in the following 


change 
changes are 
equation 
Aé6) 6A 
h 
s the zone height The change An 


the relation 
A? 


For the case of one slit located initially on the lens 
axis, h=ho(n and An, is given by the relation 


h 


specified value of the zone number m, 
be written 


h 
A? -e Ay: 


\ 


For the case of the lens axis located initially midway 
of the two slits, the zone height / is given by the 
relation 

mh. (Y) 


The value of An, is given by 


h 
A ~AfSm 
A} 


which for a specified value of becomes 
h2Afon 
\/ 
2.3. Determination of the Longitudinal 
Aberration 


Spherical 


which the 
distortion of the wave front in units of wavelength 
are determined derive the 
values of loneitudinal spherical aberration Af’, 
expressed as a focal shift. The values of Af’, 6, 
and A the zone height are related as follows: 


The same initial observations 6 from 


may also be used to 


(12) 


1 


For the ease of one slit on the lens axis at the start of 


observations, eq (Ss becomes 


(13) 


and for the case of the zero point of observations 
occurring midway of two slits, eq (S) becomes 
6} 


Ay’ 
ml 


3. Method of Measurement 


The optical bench, shown in figure 1, is used in 
making the measurements required for the deter- 
mination of the distortion of the wave front, n,, 
and the longitudinal spherical aberration Af’. 
Initially the lens under test is mounted in its holder 
on the nodal slide assembly and its optical axis 
alined with the optical axis of the collimator. The 
reticle, located in the focal plane of the collimator, 
consists of a luminous cross composed of intersecting 
vertical and horizontal lines. The collimated beam 
is incident upon the lens under test and forms an 
image of the luminous reticle in the focal plane of 
the lens. The viewing microscope is adjusted using 
the transverse and longitudinal motions provided 
until the microscope is focused on the image formed 
by the lens and the vertical line of the image coincides 
with the intersection of the crosshairs in the ocular 
of the viewing microscope. 
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A diaphragm containing two parallel slits is placed | 
in front of the lens under test. This diaphragm is 
so mounted that the double slit may be moved 
transversely across the front of the lens by meas- 
ured amounts. The slits are alined parallel to the 
vertical line of the target reticle and so positioned 
that the center of the slit system will lie in the 
plane determined by the horizontal front diameter 
of the lens and the center of the target reticle for 
all positions to be occupied by the slit system 
throughout its range of transverse movement. 
At the start of a run, the slit system is brought to a 
position as near the edge of the lens as feasible. 
Light from the collimator on passing through the 
double slit and lens forms a set of interference 
fringes in the focal plane of the lens. A photograph 
of the fringe system is shown in figure 2. The 
crosshair of the viewing microscope is brought into 
coincidence with the image of the central fringe 
using the transverse micrometer movement of the 
microscope to make the setting. This setting of the 
transverse micrometer recorded. The slits are 
now moved transversely by an amount equal to 
the slit separation, ho, and the new location of the 
central fringe is recorded. This process is repeated 
until the entire diameter of the lens is spanned. 
The microscope setting made when one of the slits 
lies on the lens axis is taken as the zero point of 
the displacement, 6, and the corresponding setting 
of the double slit is taken as the zero point of 
the zone number, From the series of observa- 
tions made at the measured intervals across the lens 
diameter, the values of 6 as a function of m are 
obtained for the particular focal setting. In ad- 
dition, it customary to determine the average 
fringe separation S) in the vicinity of the lens axis 
as it can be used in eq (1) for the determination of f 
when \ is known or for the evaluation of 26 in units 
of wavelength in eq (4). 

It is recognized that the use of a collimator with 
the illuminated ‘reticle located in its focal plane 
instead of an object at infinite distance may intro- 
duce a systematic error in the determination of the 
longitudinal spherical aberration and the optical 
path difference in units of wavelength. However, 
such systematic errors may be minimized by using 
a collimator having negligible longitudinal sphe rical 
aberration. The measurements, reported in the 
following section, were made using an off-axis 
paraboloidal mirror to collimate the light from the 
target reticle. For the area of the emergent beam 
used in the measurements the error arising from this 
source does not exceed +0.01 wavelength in the 
determination of m, and hence may be regarded as 
negligible. 


is 


mm. 


is 


4. Results of Measurement 


A series of measurements was made on a photo- 
graphic objective of a type commonly used in aerial | 


photography. The lens has a measured focal length 
of 210.4 mm and a maximum nominal f-number of 
{/6.8. The slit separation was 1.453 mm and the 


| 


Figure 1. Photograph showing optical bench 


position for test. 
The lens is mounted in the holder of the nodal slide. louble slit syst 
mounted in front of the lens on a transverse micrometer slide. 
incident on tre slit system passes through the lens and forms a fringe system 
in the fecal plane of the lens under test. The fringe system is viewed with 
the microscope having a transverse micrometer movement. 


with lens 


un 


The d m is 


Paralle) light 





Figure 2. Photograph of fringe system formed in the 


plane of the lens under test. 

In the photograph, the central fringe is indicated by the vertical arrows. 
zero, first, and second order fringes are clearly visible; higher orders are pres- 
ent and may be seen through the microscope. The bright horizontal line 
is the image of t! 1e horizontal crossline in the collimator target. 


focal 


The 


effective wavelength of the light incident * on the 
lens was approximately 589.3 mu. Observations of 
5 were made in the position of best focus and for two 
positions on either side of best focus. Values of 6 
for these three runs are listed in tables 1 to 3. Ten 
observations were made on either side of the optical 
axis and are designs ited by zone number m ranging 
from —10 to +10. Actual values of the zone height 
h are shown ‘under the column headed (m—4)hy. 

The observed values of the displacement 6 of the 
central fringe with transverse displacement of the 
slit systems are shown graphically in the left hand 
frame of figure 3 which shows the variation of 6 as a 
function of the zone height / for three positions of the 
image plane. 


‘The wavelength of the incident light is controlled by a narrow band inter- 
ference filter whose peak transmittance is at A=589.3. 
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asa function of 2 t_h (m 16) hy 

S were obtained fe focal plane 0.4mm farther from the 

gthof 210.4mm,. The double 

ength of the incident light was 
erration Af’ and the distortion 
ilso given 


lens than 


4.1. Evaluation of Longitudinal Spherical 
Aberration A/’ 


It is usually convenient to evaluate first 
longitudinal spherical aberration Af’ as it can be 
done readily using the observed values of 6 and 
eq (13). The values of Af’ are listed in tables 1 to 
3 and are shown graphically in the right-hand frame 


the 
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of figure 3. It is clear that the form of the curves is 
essentially the same for the three focal positions but 
the curves are displaced with respect to each other by 
amounts closely approximating the nominal separa- 
tions of planes in which the measurements were made. 


TABLE 2. Displacement 5 in the image plane of central fringe 


as a function of zone height, h=(m bo) hy 


These values were obtained for the plane of best focus ofa lens having a focal 
length of 210.4 mm rhe double slits are separated by ho=1.453 mm.  Wave- 
length of the incident light was 589.3mm,. Valuesofthe longitudinal spherical 
iberration Af’ and the distortion of the wave front, na, in units of w ivelength 
re also given 


TaBLE 3. Displacement 6 in the image plane of central fringe 


16) hy 


These values were obtained for a focal plane 0.8 mm nearer to the lens than the 
plane of best focus of a lens having a focal length of 210.4mm. The double slits 
ire separated by ho=1.453mm. Wavelength of the incident light was 589 3mm 
Values of the longitudinal spherical aberration Af’ and the distortion of the wave 
front, , In units of wavelength are also given 


as a function of zone height, h (m 


4.2. Check on Consistency of Observed Data 


In work of this nature it is desirable to have and 
use check methods to ensure that errors are held to 
a minimum and to verify the selection of the form of 
equation used in the determination. This can be 
done by analysis of the variation 6, the average of 6, 





for paired values of +m and —™m, for two focal posi- 
tions. Consideration of eq (12) indicates that the 
difference Ad obtained from the differences in 6 for a 
given value of m at two different focal positions is a 
linear function of m. In table 4, values of Aé obtained 


by taking the differences between the corresponding | 


values of 6 for focal positions +0.4 and —0.8 are 
listed under the caption, observed values of Aé. 
When these values of Aé are plotted against zone 
number m as shown in curve 1 of figure 4, it is clear 
that a linear relation exists between Aé and m—'. 
This indicates that it was proper to use eq (13) in the 
evaluation of Af’. The slope a of the best fitting line 
through the points on the graph is given by the 
expression 


DAs 
a== 
=(m 


t) 
which in this instance gives a value 


a=383/50=7.66 pu. 


The equation of the best fitting line is 


Ab. =a(m—}). 














ZONE NUMBER, m 


Figure 4. 
fringe displacement Ad between two specified focal planes 
with zone number, m. 


Variation of the change in average transverse 


Curve 1 shows the variation in Aé for focal planes designated Af=0.4 and Af=—0.8; 
curve 2 shows the variation in Aé for the focal planes designated Af=0.0 and 
0.8; and curve 3 shows the results for focal planes designated Af=0.4 and Af=0.0 
mm. The equations of the best fitting straight line are shown adjacent to 
each curve. 


| Values of Aé, obtained with this equation for 
a=7.66 uw are listed in table 4 together with the de- 
partures from the observed values of Aé. It is clear 
that the agreement is excellent. Values of Aé for the 
other two combinations of focal position (+ 0.4 and 
0.0, curve 3; 0.0 and —O.8, curve 2) are also shown 
in the graph and indicate that all three sets of data 
are reasonably consistent. 


TABLE 4. Comparison of values of 46 obtained for focal posi- 
tions +0.4 and —0.8 as a function of m with values computed 
from relation Aé 7.66 (m—}) 


4é in microns 


Observed Comp 


From the values of the slope a, for the three com- 
binations of Aé versus m, it is possible to determine 
the average displacement of the focal plane for the 
three focal positions usine the relation 


at 


ho 


A(Af) (17) 


Values of A(A/) derived from the analysis of Aé for 
the three combinations are given in table 5. The 
focal positions were known initially within + 0.1. mm 
so it is clear that good agreement exists between the 
nominal values of A(Af) and those derived from 
analysis of Aé as a function of m. 


TABLE 5. Comparison of nominal values of the relative se para- 
tion A(Af) of the three focal planes with that derived from the 
slopes of Ad versus mM lo and with that derived from the slope 
of Any versus m? 


Value of A(Af 
Focal position Nominal 
From 
Am VS 


From slope of 
Aé VS m 


0.0 and 0 
0.0 and 0 
0.4and —0. 


4.3. Evaluation of Optical Path Difference, n, 


The values of the optical path difference, n,, are 
obtained with the aid of eq (4). Here the value of 
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the term h,/A fis equal to 0.01172 wavelengths/micron, 
Values of 26 and n) for the three focal positions are 
shown in tables 1 to 3. The values of n, are also 
shown graphically in figure 5. 


4.4. Check on Consistency of Values of ny 


While it is possible to evaluate n, from data ob- 
tained for a single focal position, confidence in the 
results is increased by comparison with results for 
a different focal position. Hence it is desirable to 
check the inner consistency of results obtained for 
several focal positions to ensure that errors have been 
held to a minimum and to verify the form of equation 
to be used in determing the change in n, with shift of 
focus. This can be done by analysis of n, the 
average value of n, for paired values of +m and 

m for two focal positions. Consideration 
of eq (8) indicates that the differences An, obtained 

from the differences in n) for a given value of m at 

— two different focal positions is a linear function of m?. 

n,, WAVELENGTHS In table 6, values of An, obtained by differencing the 
corresponding values of mn, for focal positions +-0.4 
and —0O.8 are listed under the caption, observed 
values of Any. When these values of An, are plotted 
against m? as shown in curve 1 of figure 6, it is clear 








Figure 5. Variation of the optical path difference expressed 


in number of wavelengths, n, as a function of zone height h. 


of nm, are shown for three focal positions. Curve 1 shows the results for . ° ° 7 > 
4, curve 2 for Af=0.0, and curve 3 for Af=—O.8. These values of mare | that a linear relation exists between An, and mi’. 


ere ao This also indicates that eq (8) was the proper choice 
for use in the interpretation of the relation between 
An, and m. The slope 6 of the best fitting line 
through the points on the graph is given by the 
expression 


DAnr 


y 2 
wih 


b 


which in this instance gives a value of 


TH 


b 17.53/385=0.04553 


The equation of the best fitting line is 


Sn), WAVELENG 


Any bm (19) 


Values of An, obtained with this equation for 
b=0.04553 are listed in table 6 under the caption, 
computed values of An), together with the departures 
from the observed values under the caption, 0—C. 
It is clear that the agreement is excellent. Values of 
An, for the other two combinations of focal position 
(0.4 and 0.0, curve 3; 0.0 and —0.8, curve 2) are also 
shown in the graph and good consistency is indicated 
. for the three sets of data although the values of 
Ficure 6. Variation of the average change An, between two | Q—C are larger at m=9 and m=10. 
specified focal planes with m2, where m is the zone number. From the values of the slope } for the three 
Carve 1 shows the variation in Am for f cal planes designated Af=0.4 and af combinations it is possible to determine the average 
5 m irve 2 shows the variation for the focal planes designated Af=0.0 : 
3 shows the variation for focal planes designated | displacement of the focal plane for the three focal 
-quations of the best fitting straight lines are shown positions using the relation 
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2rf? 


~«b, 20) 
h,? 


A(Af) 


Values of A(A/) derived from analysis of An for the 
three combinations are also listed in table 5. It is 
clear that good agreement also exists between the 
nominal values of A(Af) and those obtained in the 
foregoing manner. 


TABLE 6. 
position 


Comparison of values of An, obtained for 
+0.4 and 0.8 as a function of m with 


computed from the relation An, =0.04553m? 


focal 
those 


Values of An) 


(‘or iputed 


4.5. Effect of Focal Shift Af on the Values of ny 


One of the manifest advantages of this method of 
measuring optical path difference ny, is that it per- 
mits the evaluation of n, as a function of zone height 
h for any given tocal position separated by Af from 
the focal plane in which measurements were made. 
Hence it is not necessary to make the measurements 
of 6 at the position of best axial focus. It is only 
necessary that the focal plane in which measurements 
are made is that of a reasonably good focus as the 
form of the curve can be computed for a selected 
focal plane if one has the measured values of n as a 
function of m tor a nearby plane. The suitability 
of such a procedure is illustrated by comparison of 
the 3 sets of values of n,, the adjusted values of 
n, Obtained from n., the observed values of n, for 
the three focal positions for which results are given 
in tables 1, 2, and 3. 

Values of no, the observed values of ny, as a function 
of zone number are given for the three focal positions 
in table 7. It is clear that the values differ by 
appreciable amounts. It is worthwhile to adjust 
these values by a focal shift of such amount that n,, 
the adjusted values of ny), shall be zero at m=7. 
This may be done using the relation 


Any n, 
where 


, 


An) (45)m: (21) 


and n)’ is the observed value of n, for zone number 
m=7. Values of n. obtained in this manner are 


also listed in table 7 together with n,, the average for 
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the three focal planes. It is clear that there is good 
agreement among the three values of n, although 
the determinations are based upon measurements 
made at markedly different focal positions. Con- 
sideration of the departures from the average n, 
shows that the average error in n, does not exceed 
+ 0.02 for the range from m=1 to m=8 and does 
not exceed +0.05 for values of m=9 and m=10. 
TABLE 7. Values of ne as a function of zone number 
three focal planes 
Values of n, obtained from no, the observed values of na, for three focal 


Values of Ana are so chosen that n.=0 at m=7. The value of n-., the 
n- for the three focal planes, for each value of m is also given 


Af=0.4 mm Af=0.0 mm 0.8 mm 


An 


0.01 0.03 0.03 
O4 09 09 
Ov 1s 18 
16 2 22 
25 2 19 
36 
49 
tA 


SI 


The average value n, versus zone height h is 
shown in figure 7. It is clear that for this particular 
focal position the optical path difference does not 
exceed % wavelength for the range from m=1 to 
m=8. From this it may be inferred that the optical 











1.0 


Nic, WAVELENGTHS 


average adjusted value, n with 
m is the zone number and 


Figure 7. Variation of the 
zone height h=(m—})ho where 
ho the slit separation. 

The value of 7; is the average value of n) for the three focal planes designated 
Af=0.4, 0.0, and —0.8 following adjustment such that n-=0 at m=7 in all three 
cases. The scale on the right-hand ordinate of the graph shows the value of 
the f/number corresponding to the indicated value of h, the zone height. 





path difference will not exceed 4 wavelength for one- 
half the effective area of the lens at maximum aper- 
ture. Moreover, it can be inferred that this lens 
should yield excellent imagery at f/11, good imagery 
at f/10, and that the quality of imagery can be ex- 
pected to lessen appreciably at f/8 and smaller 
f/values. 


5. Discussion 


In this study, a visual method of determining 
longitudinal spherical aberration and departure of 
the wave front erriving at the image point from a 
true sphere is described. Results are reported for 
measurements made on a single lens for three posi- 


tions of the focal plane. Tests for consistency of 


observations are described which permit an estimate 
of the reliability of the measurements. 


The authors express their appreciation to other 
members of the staff for assistance during this work 
and particularly to Edgar C. Watts who designed the 
slit system and prepared the illustrations; and to 
George Cassiday and Peter Costello who performed 
many of the measurements. 
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have 
but because of the 


Filters 
colorimeters, 
the other, 


the CIE tristimulus functions. In 


been designed for use with photocells and phototubes in 
variation in spectral response of these detectors from one to 
it has not been worthwhile to attempt to achieve the best possible duplication of 
the case 


1963) 


tristimulus 


of a thermoelectric detector, however, the 


response is so nonselective that the same filter designs can be used with any thermopile. 
The new filters permit CIE tristimulus funetions and chromaticity coordinates to be 


obtained directly with reasonable accuracy. 


Results obtained on a thermoelectric 


colorim- 


eter for five standard filters are compared with tristimulus values obtained by computation 


from spectral transmittance data. 


1. Introduction 


Experience has shown that the results obtained in 
visual color measurements by a single observer are 
likely to differ from those of some other observer. 
For this reason, a number of observers must be used 
in order to obtain representative results. Further- 
more, an Observer using the equality-of-brightness 
method of photometry will change his criterion for 
photometric balance over a period of several days; 
some Observers exhibit variation during a single day. 
If the properties of the normal average eye (CIE 
standard observer) can be reproduced in a completely 
physical instrument, these problems can be avoided. 
However, it is to be emphasized that physical color 
measurements must be made with high precision if 
they are to compare favorably with visual judgments 
of color, especially of small color difference. A 
colorimeter comparable to the eye of the standard 
observer in discrimination sensitivity must be capable 
of measuring tristimulus values, Y, Y, Z, to at least 
one part in athousand. For many years, the Bureau 
has been examining glasses and chemicals for possible 
use in physical colorimeters. It is more desirable to 
design tristimulus filters for use with a nonselective 
detector rather than a selective detector because such 
filters are not dependent on the useful life of an indi- 
vidual detector nor limited to a particular one. 

No physical detector has been constructed which 
has the same relative response curve as the average 
human eye for radiant energy of different wave- 
lengths. It is possible, however, to make a non- 
selective receiver, which indicates the amount of 
received energy per unit time independent of its 
wavelength throughout the visible spectrum. If we 
interpose between the source and such a receiver 
filters which duplicate the CIE tristimulus functions, 
we can then measure the tristimulus values and by 


*Guest Worker 
Egypt. 


it the NBS—On leave from Ein-Shems University, Cairo, 


computation obtain the chromaticity coordinates 
without reference to human observers. Such filters 
would permit CIE tristimulus functions and chro- 
maticity coordinates to be obtained directly. The 
object of the present work was to determine if such 
filters could be realized closely enough to determine 
these quantities with res asonable accuracy. 

Numerous instruments have been proposed and 
many constructed for direct color measurement over 
the vears. Ives (1915) [1]! deseribed a thermopile 
tristimulus colorimeter in which he made use of a 
spectrum selectively projected through a series of 
carefully cut templates in place of the filters employed 
in the more recent devices. Twyman and Perry 
(1930) [2] suggested that three filters duplicating the 
then proposed standard observer could be used in a 
photoelectric instrument. Guild (1934) [3] described 
an experimental model of a photoelectric tristimulus 
colorimeter. A few years later Winch and Palmer 
(1937) [4], Dresler and Frithling (1938) [5], and 
Barnes (1939) [6] described ohiyg for the same 
purpose. Gibson (1936) [7] and Van den Akker 
(1937) [8] have discussed the problem of obtaining 
source-filter-photocell combinations giving the closest 
equivalence to the standard observer. Hunter, 
(1942) [9] described a three-filter method for photo- 
electric tristimulus colorimetry which was success- 
fully applied to the Hunter multipurpose reflectom- 
eter [10] to measure surface color. Several new 
instruments have been designed in recent years.’ 
However, the photosensitive elements in these 
colorimeters have been barrier-laver photocells as in 
the instruments of Barnes (1939) [6] and Hunter 
(1942) [9], phototubes in the devices described by 
Glasser and Troy (1952) [11] and Hunter (1958) [12], 
and multiplier phototubes in the instruments de- 
signed by Sziklai (1951) [13], Bentley (1951) [14], 
and Nimeroff and Wilson (1954) [15]. Only in the 
instrument described by Ives (1915) [1] was a ther- 
mopile used. 
ature references at the end of this paper. 


| Figures in brackets indicate the liter 
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2. Filter Design 


The only parts of a colorimeter over which the 
designer can exercise much control are the filters 
used to modify the detector response to obtain the 
%, j, and 2 functions of the CIE standard observer. 
The @ function is bimodal and therefore two filters 
are used; one designated Z#-short represents the 
function from 380 to 500 nm and the other, Z-long, 
the function from 510 to 760 nm. The components 
of these filters are usually absorbing solutions and 
Information on characteristics of various 
solutions and glass cells tested for photometric use is 
given — papers published by Ives and Kingsbury 
{16}, | [18]. In this investigation, spectrophoto- 
me frie curves of various solutions and glasses were 
made ona General Electric recording spectrophotom- 
eter. The data obtained on these filter components 
were used to compute the concentrations and thick- 
nesses required to obtain a good fit for the spectral 
tristimulus values, Z%, ,, 3. Each of the filters 
designed contains as its components a solution, a cell, 
and two or more glasses. 


ols nsSses. 


2.1. Solutions 


Chemicals used for preparation of solutions were 
labeled “analytical reag quality” or “reagent 
quality.” To prevent contamination from cork 
stoppers, glass-stoppered bottles were used as con- 
tainers for all solutions. 

The compositions of solutions used are as follows: 


rent 


a. Solution I (x-short and z) 


Cobalt ammonium sulphate 
CoSO,-(NH,).SO,-6H:0 13.0 § 
Nickelous sulphate (NiSO,-6H,O 19.5 ¢g 
Glacial acetic acid 60 - 
Distilled water to make one liter of solution. 


« 
< 


b. Solution II (x-long) 


Copper sulph: ite (CuSO,- 
Ferric sulphate (Fe,(SO,);-9H,0) ll.ig 
Glacial acetic acid 60 ml 
Distilled water to make one liter of solution. 


5H.O) 11.0 ¢ 


c. Solution III (y) 


Potassium dichromate (K,Cr,O- 
Copper sulphate (CuSO,-5H,O) 
Cobalt ammonium sulphate 
(CoSO,-(NH,).SO,-6H.O) 210 ¢ 
Sulphuric acid (1.835 sp gr) 10 ml 
Distilled water to make one liter of solution. 


0.140 
l7.0¢ 


2.2. Cells 


The filter accuracy obtained with solutions in cells 
depends in part on the spectral transmittance and 
uniformity of the cell. Cells having nonselective 
windows, closely parallel, and with close tolerances 


on window spacing are commercially available. 
Three identical cells with 1-mm thick windows and 
50-mm spacing were selected. Spectrophotometric 
measurements indicated that the windows were 
nonselective in transmittance from 380 to 780 mm 
except for variations than 1.5 percent. The 
solutions I, II, and III prepared as described above 
were placed in cells I, Il, and III, respectively. 


less 


2.3. Glasses 


evlasses which 


The specifications of the 
selected are shown in table 


were 


TABLE 1. Glass components of the four tristimulus filters 


M fer ‘olor 
No Vame spec Chickness 
N 


mm 
Heat resist yrojec bli 5-50 4.4 
Heat. resistant lantern blu 5-0 0.84 
4-97 615 
3-80 67 


Heat absorbing K-48 *5 


about 5 mm as furnished b manufacturer, This glass 
nay tony han » obtainable from Kopp Gla ‘o., but it is similar to Corning 
filter No 4000, ce] ir specification No, 1-69 spectral transmittance of these 
¢ absorbing glasses from 560 to 780 i ir to that of copper sul- 
pr rate solution rt , adjustment ration could be used te com- 
pensate for individual differences in spectral transmittance of the glass; however, 
l l ined ym three manufacturers were found to be sufliciently 


solution concentration was required 


2.4. Tristimulus 


Solution-cells and the glasses listed in table 
combined to provide tristimulus filters as follows: 


ure 


z-short Filter 


Solution-cell I and glasses numbered 1, 2, and 5 


b. v-long Filter 


Solution-cell II and glasses numbered 3 and 5 


c. » Filter 


Solution-cell IT] 
This filter design is a slight 
luminosity filter design by Gibson, Teele, and 
Keegan (1939) [19]. The use of heat-absorbing 
glass necessitated a reduction in copper sulphate and 
potassium dichromate and an increase in cobalt 
ammonium sulphate. 


and 5. 
modification of the 


and glasses numbered 


d. = Filter 


Solution-cell I and glasses numbered 1 and 5 


3. Filter Characteristics 


Spectral transmittance data for each of the four 
filters were obtained on a recording spectrophotom- 
eter [20], [21], [22]. Routine corrections for 100 
percent, zero, and wave length scale were made [23] 
but additional corrections for inertia, slit width, and 
back reflectance were not made [24] since their 
total for these filters is less than the uncertainty of 
a single run on the spectrophotometer. The results 
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obtained are summarized in table 2. For each filter, | Where the subscript for I corresponds with the 
the first column gives the spectral transmittance, | particular colorimeter filter being used. 
the second column gives the relative transmittance 
adjusted to have the same area as the CIE curve, 
and the third column gives the CIE data. For ease 
of comparison the data from columns two and three 
for each filter are shown in figures 1, 2, 3, and 4.) , cera 
The transmittances less than 0.001 for the 2 color- ae 
imeter filter from 590 to 830 nm were found by ce 

Meg: Be: . OBSERVER 
remeasuring the three individual components of the 
filter on a wider range spectrophotometer after O EXPERIMENTAL 
unexplained discrepancies turned up in the check of 
the performance of the colorimeter to be described 
presently. 

For each colorimeter filter the sum of column 
three divided by the sum of column one gives the 
factor used to obtain the values in column two. 
These factors are also used to obtain the relative 
tristimulus values of any light source directly from 
the instrument readings, I. 








J i i 
20 40 60 80 
WAVELENGTH, NANOMETERS 





2.560 


x 
Y =2.033 
Z 


FiGuRe 1. Adjusted transmittance of x-short filter compared 


12.166 . to CIE tristimulus values X). 


TABLE 2. Veasured , T, and adjuste d spectral transmittances for each of the four filters con pare dto CIE 
spectral tristimulus values 


T-long 
Wave 
length 
5.5963 T IE Dy 2 j E 4 12.166 T 


0.0014 0. 0000 0. OOOO OO00 0. OO06 
0042 0007 0014 o0oo1 0034 


0143 0007 0014 0004 0074 

0435 0007 0014 ool 0199 2421 
1344 Ooll 0022 0040 0608 7397 
2839 0030 0061 O116 1242 5111 
3483 0092 O87 0230 1623 9746 


3362 O87 0380 0380 1414 7203 
2908 0297 0604 0600 1094 3310 
1954 0523 1063 0910 O72 0609 


0956 0. 0000 0. 0000 0731 1486 1390 0829 7652 
0320 0016 0090 1020 2073 2080 0394 4793 


0049 0046 0257 1530 3110 3230 0209 2543 

0431 ¢ 2437 4054 5030 0106 1289 

0677 d . 3545 7206 7100 0053 0645 

530) 7 212 1186 5 4250 8639 8620 0026 0316 
540 0005 2317 2¢ 4671 9495 9540 0021 0255 


550 0010 4007 3: 1925 ool! 9950 0036 0438 
560 OO1S 6044 504! 4038 0037 9950 OO5S 0706 
570 OO13 8355 76 46592 9537 9520 0049 0596 
580 oool 7 OO17 Ot 4263 S665 S700 00113 0138 
590 0000 1. 0599 : 3716 7554 7570 00021 0026 


600 1. 0297 522 3077 6255 6310 00013 0016 
610 0. 9267 26 2443 4966 5030 00010 0012 
620 ‘ 7689 f 1832 3724 3810 00007 0008 
630 6033 42 1297 2636 2650 00003 0004 
640 4354 78 O853 1734 1750 00001 0001 


650 2972 t 0544 1106 1070 00001 0001 
660 1892 ( 0316 0642 0610 00001 0001 
670 1186 7 0171 0348 0320 00002 0003 
680 2 0677 j 0090 0183 .0170 00005 0006 
690 0375 227 0043 0087 0082 00012 0015 


700 0196 0013 0026 0041 . 00019 0023 
710 0106 5 0006 002 0021 00016 0020 
720 0056 2 0003 0006 0010 00013 0016 
730 0006 0001 0002 0005 00010 0012 
740 0000 7 0001 0002 0003 00008 0010 


750 0003 0000 0000 0001 00007 0009 
760 0002 0000 0000 0001 00006 0007 
770-830 0001 00020 0025 


Sums 0. 6971 78 SO86 8. 8O87 5. 2569 10. 6855 10. 6857 0.8751 10. 6770 10. 6770 
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STANDARD OBSERVER 


RIMENTAL 








FIGURE 2 of x-lonq filter compared to 


€S X). 


CIE STANDARD OBSERVER -— 
EXPERIMENTAL 





filter compared to 


y 
vn. 


4. Performance 


To check the performance of these four color- 
imeter filters, use was made of the instrument de- 
scribed by Teele |24| which consists essentially of a 
thermopile, potentiometer, and recorder. 

The first check was a measurement of the chro- 
maticity coordinates of an incandescent lamp at a 
color temperature of 2,854 °K (CIE source A). 
The tristimulus values obtained by reducing these 
data by the factors given above were divided by 
their sum in the usual way to obtain values of chro- 
maticity coordinates as follows: s=0.457, y=0.403, 
z=0.140. These chromaticity coordinates differ 
from those obtained by cemputation from the 
spectral distribution of a Planckian radiator at a 
temperature of 2,854 °K by +0.009, —0.004, and 
—0.005, respectively. These differences correspond 
to a combination of errors from possible failure of 
the voltage assigned to the lamp to operate it at a 
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40 60 80 





WAVELENGTH, NANOMETER 


Adjusted transmittance 


CIE tristimulus 


FiGuRE 4. compared lo 


color temperature of 2,854 °K, from failure of the 
lamp to yield spectral distributions characteristic of 
the Planckian radiator, from failure of the color- 
imeter filters to duplicate the CIE standard observer, 
and from the experimental errors of reading the 
thermopile. 

A more detailed check of the performance of the 
colorimeter filters was obtained by measuring five 
NBS standard spectrophotometer-integrator filters, 
NBS 2101 through 2105, master set number 1. 
These are glass filters described as selenium reddish 
orange, carbon yellow, sextant green, cobalt blue, 
and selective neutral. 

They were originally calibrated to check the per- 
formance of spectrophotometers equipped with tri- 
stimulus integrators. [25] The five filters were 
measured on two spectrephotometers (Cary 14 and 





GE); appropriate corrections (such as wavelength, 
zero, 100%, slit width, inertia, and back reflectance) 
were applied; spectral transmittances were adopted 
as weighted means; and uncertainties were esti- 
mated. The spectral transmittances of these filters 
are plotted in figure 5. 

To provide a more detailed test of the performance 
of the colorimeter filters use was made of a device 
already devised by Teele, which consists essentially 
of a thermopile (nonselective receiver for radiant 
energy), potentiometer, and recorder, the incan- 
descent lamp was operated at a color temperature of 
2,854 °K as before. To eliminate back reflection 
effects, a concave-convex lens was positioned before 
and another after the filter to be measured as shown 
in figure 6. The colorimeter filter, say z-short, was 
then introduced between the light source and thermo- 
pile, and the reading was recorded. The filter to be 
measured was placed in the beam between the con- 
cave-convex lenses, and a second reading was 
recorded. This procedure was repeated with each 
of the four colorimeter filters. The ratios, source 
through filter to source alone, of such readings 
through the four colorimeter filters, Z-short, Z-long, 
7, and 2, gives the values: Re snort, Re-1ong, Ry and RZ 
respectively, which may in themselves serve as a 
colorimetric specification of the filter to be measured 
relative to the reference source. The CIE tristimu- 
lus values of the test filter with source A, adjusted to 
make Y, equal to the luminous transmittance of the 
filter in percent, are computed as follows: 

X .=5.37R3-woret 


104.4525. tong, 


Y,=100.00R;, 


Z,.=35.55R;, 


where the coefficients 5.37, 104.45, 100.00, and 35.55 
are related to the source-A tristimulus values. The 
coefficients of R- and Rt, are the Y and Z tristimulus 
values for source A, respectively. The coefficients 
Of Re-snorr ANd LZ_iong are the summations of the spec- 
tral tristimulus values for the Z function and source 
A from 380 to 500 nm and from 510 to 770 nm, 
respectively. 

The chromaticity coordinates, x, y, 2, of the test 
filters are calculated from the tristimulus values, LY, 
Y, Z, in the usual manner: 
2=X/((X+Y+42Z), y=V U(X +Z),2=Z/((X+Y+Z) 

Table 3 compares the measured values of the five 
standard filters with their assigned values. Each 
measured value was computed from 10 independent 
determinations of the ratio R. The uncertainties of 
the assigned values shown in the table are those esti- 
mated by Keegan, Schleter, and Judd as the ranges 
of the values found by three different spectropho- 
tometers |25|; see their table 11. The uncertainties 
shown for the measured values are based upon com- 
putations of huge error (4.9 times the probable error) 
of the mean of the ten settings of the ratios R. These 
computations showed that the huge error (uncer- 
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Figure 5. Spectral transmittances of the five test filters 
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Fiaure 6. Diagram of the optical arrangement. 


tainty at the one-tenth of one percent level) of values 
of the ratio 2 was on the average equal to + 0.0025, 
and these estimates of uncertainty are given directly 
for the values of Y in table 3. Propagation of these 
errors to chromaticity coordinates, z, y, 2, leads to the 
estimates of uncertainties entered in table 3. 

It will be noted from table 3 that the sum of the 
uncertainties listed fails in 11 out of 20 to 
exceed the difference between the measured and as- 
signed values. Of course, since the colorimeter filters 
do not succeed precisely in duplicating the distribu- 
tion functions of the 1931 CIE standard observer 
(see table 2), some discrepancy between measured 
and assigned values is to be expected. Table 4 gives 
in column 2 the differences ascribable to failure of 
the spectral transmittances of the colorimeter filters 
to have precisely the wavelength variation striven for. 
In column 3 are given the observed differences copied 
from the last column of table 3, and in column 4 
are given the sums of the uncertainties listed in col- 
umns 2 and 3 of table 3. If the expected difference 
departs from the observed difference by less than the 
combined uncertainty, we may say that a reasonable 
basis for explaining the discrepancy has been foynd. 
Column 5 gives the amounts of the unexplained dis- 
crepancies. It is computed as the excess of the abso- 


cases 


- 
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lute value of the difference between columns 2 and 
3 over the sum of the uncertainties in column 4. If 
there is no excess, no entry is inserted. 


TABLE 3 
chro? aticit 


Y. for 


Con parison of measured and assiqned values of 
nates, x, v, z, and luminous transmittance, 


A for the VBS standard filters 


/ coord 


CIE source 


+0. 0060 
0059 
0064 


es between the differences 


expected fron the actual 


jained discrep- 
If the absolute 
the difference 
columns 2 
ind 3 exceeds the sum 
of the uncertainties in 
column 4, the amount 
of this excess is recorded 
below, 


value of 


between 


0. 0027 0. 0074 
0027 0056 


0008 0010 


0. 0054 
0026 
-, OO81 


+0). O17 
014 
003 

-. 51 


+-() 0121 
0021 
0100 


1 


+0). OO48 0. 0033 

OO11 0034 
—, 0038 0021 0010 
—.6 45 7 


Table 4 


accounted 


shows that nearly all 
for by the imperfection of the fit of the 


discrepancies are 
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four filters to the respective CIE curves and that 
improvement in performance depends upon finding 
better constitutents for the solutions and different 
classes (or different thicknesses of the same glasses 
that are now used). 

It will be noted from table 4 that for four of the 
five filters the unexplained discrepancies (column 5) 
are few and small. For filter 2101, however, the 
discrepancies in chromaticity coordinates are large. 
Perhaps it is significant that this filter (selenium 
reddish orange) is the only one of the five which has 
a large temperature coefficient of spectral trans- 
mittance. 

By comparing the uncertainties in the measured 
values (column 3, table 3) ascribable to the uncer- 
tainty of 2.5 on a scale of a thousand in the measured 
ratios with the differences (column 2, table 4) to be 
expected because the colorimeter filters do not have 
precisely the ideal spectral character striven for, it 
may be noted that those two sources of error are 
about equally important. The performance of the 
colorimeter could be significantly improved only by 
developing better # and 2 colorimeter filters and at 
the same time reducing the noise from the thermo- 
pile and amplifier. Significant improvement of the 
colorimeter filters might be accomplished by com- 
bining the components in parallel as well as in series 
as Was originally suggested by Dresler [26], and used 
among others by Nimeroff and Wilson [15], Friihling 
and Krempel [27], Geutler [28], and Davies and 
Wyszecki [29]. 


The authors express their appreciation to the 
members of the staff who ran the numerous spectro- 
photometric curves required to develop the filters. 
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In an attempt to eliminate unwanted resonances and transverse motions in an electro- 
dynamic exciter used as a vibration standard, the mechanical suspension-guide for the moving 


element was replaced with air bearings. The 


transverse motions of this exciter showed a 


distinct improvement over the transverse motions of the exciters with flexure plate and ten- 


sioned wire supports. The 


ment between three methods of calibration, 
was very £00 1. Ther 
are presented, 


1. Introduction 


The widespread use of vibration pickups to meas- 
ure the amplitudes and frequencies of vibrations 
present in aircraft, rocket engines, ordnance systems, 
machinery, ships, ete., has led to the need for im- 
proved vibration standards. The { feasibility of using 
certain kinds of electrodynamic vibration exciters as 
standards for the calibration of pickups was demon- 
strated by Levy and Bouche [{1, 2]'. A limitation in 
the use of electrodynamic vibration exciters as 
calibrators has been the presence of undesirable trans- 
verse motions superimposed on the required sinus- 
oidal axial motion of the moving element. 

The present report the modifications 
made in a commercial electrodynamic exciter in 
which the mechanical suspension-guide for the mov- 
ing element is replaced with air bearings. The 
methods used and the results of the calibration of 
this exciter are given and compared with exciters 
using flexure plate and tensioned wire supports. The 
results of calibrations of two types of vibration pick- 
ups using this exciter are also presented. 


describes 


2. Description of the Electrodynamic Exciter 


A rectilinear vibration exciter of the electro- 
dynamic type was equipped with air-bearing support 
guides. The moving element includes a_ velocity 
sensing coil which, when calibrated, permits the 
exciter to be used as an electrodynamic vibration 
standard. The exciter is shown schematically in 
figure 1. The moving element consists of a mount- 
ing table, 7’, a velocity sensing coil, C,, a driving 
coil, C,, and a connecting shaft. The table is set in 
motion by applying alternating current to the driv- 
ing coil which is in a steady magnetic field induced 
by the main magnet, M, when direct current is 
passed through the field coil, C,,.. The motion of the 
table, 7, is measured by the voltage output of the 
velocity sensing coil, C,, in conjunction with the 


1 Figures in brackets indicate the literature references at the end of this paper 


improved stability of the moving element of the air-bearing 
exciter made it possible to calibrate this exciter easily by an interferometer method. 
reciprocity, 
results of calibrations of two types of pickups on this modified standard 


Agree- 


interferometer, and optical target, 


permanent magnet, W, since this voltage is 
portional to the absolute velocity of the 
motion. The pickup to be calibrated is 
mounted on the table, 7. 

The ideal exciter should produce pure axial motion 
free of transverse motions. The earlier use of 
flexural plates or suspension wires caused unwanted 
resonances and transverse motions. 

More detailed discussions of the principles of 
operation and limitations of the electromechanical 
system are given in references [{1, 2, and 3]. 


pro- 
table 
rigidly 


2.1. Air-Bearing Standard 


The vibration exciter modified with an air-bearing 
suspension guidance system is shown in figure 1. 
The moving element is supported on two air bearings 
of the external pressure or hydrostatic type [4-7]. 

The only mechanical connections between the 
moving element and the stationary parts of the ex- 
citer are the flexible input wire to the driver coil and 
the thin output wire from the velocity sensing coil. 
The air bearings are located in the approximate posi- 
tion of the replaced flexure plates. Two different 
sizes of bearings are used. The front bearing is 2.187 
in. in diameter and the back bearing 1.250 in. The 
lengths of the bearings are 1.000 in. The journal of 
the front bearing is also the vibration table of the 
exciter and was designed with a larger diameter to 
make it possible to mount large or small vibration 
pickups on it. The front bearing, shown in figure 2, 
has twelve orifices; the back bearing has six. The 
bearing sleeves have air distribution annuli 0.125 in. 
deep and 0.255 in. wide in their outer surface. Ori- 
fices of 0.02 in. diam. are drilled through these sleeves. 
The bearings have a radial clearance of 0.001 in. 
The bearings are made of brass and have aluminum 
journals. 

The front bearing fixture, which is bolted down on 
three studs, is a 1 in. thick circular aluminum plate 
14 in. in diameter. The back bearing is fixed in the 
same manner on the magnetic shield which is a cir- 
cular steel plate in. thick and 8 in. in diameter. 
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FIGURE 2. 





bearing used to provide one degree of freedom 


in the electrodynamic exciter. 


The velocity coil was modified by cutting off the 


reduced diameter front part. It then was screwed 
and cemented rigidly to the table and the shaft. 
The compressed air used is the building mainte- 
nance airsupply. The air is filtered and each bearing 
has a separate air regulator (see fig. 1). The air is 
transmitted to the bearings by means of flexible tub- 
ing with an inside diameter of 0.125 in. During the 
calibration of the exciter and the pickups, the air 
pressure was maintained at 50 psig for both bearings 
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Figure 3. Setup for calibration of a piezoelectric accelerometer. 
and the air temperature was about 65 °F. When 
varying the air pressure between 20 and 80 psig, no 
significant change is observed in the measured 
transverse motions. 

When the exciter is used in a vertical position, the 
moving element is supported by three rubber bands 
120° apart attached to the moving shaft by screws 
and to a circular plate positioned directly above the 
moving element as shown in figure 3. The circular 
plate can be rotated and the position of the moving 
element thus adjusted with respect to the stationary 
exciter. The equilibrium positions of the coils are 
held constant with respect to the magnets. With 
the exciter in a horizontal position, there is no need 
for any additional support. The leads from the 
driver coil provide sufficient support to keep the 
moving element from rotating and coming into 
physical contact with the stationary part of the 
exciter. An ohmmeter is connected between the 
moving element and the sleeve to indicate if metal 
to metal contact is made between the moving and 
nonmoving parts of the exciter. By careful aline- 
ment of the two bearings, this possibility is eliminated 
because the surfaces are always separated by a thin 
film of air. Having the exciter in a horizontal posi- 
tion with a 5 lb weight mounted on the table, and 
air pressure at 20 psig, no mechanical connection is 
observed. A 5 psig air pressure is adequate to sup- 
port the moving element with no load attached. The 
performance data of the air bearing standard are 
given in table 1. 





TABLE 1. Performance data for air-bearing exciter 


Weight of moving element 2.20 Ib. 

Exciter table load Sib. 

Force rating {0 to 50 Ib to 500 c/s. 
\0 to 25 lb to 1000 ¢/s 

Pounds per ampere 10 

Maximum table excursion 0.5 in. 

Frequency of major resonance Z 2200 to 2300 c/s. 

Frequency of iviagnet structure resonance 10 c/s approx 

1-c fleld resistance 50 to 55 ohms approx. 

Driver coil resistance 2.0 ohms approx. 

Stray Magnetic fleld at table level 5 gauss approx 


2.2. Transverse Motion Measurements 


The transverse motions of the mounting table were 
measured using two piezoelectric accelerometers 
with their principal sensitive axes mutually per- 
pendicular to each other as shown in figure 4 [2]. 
The one accelerometer, which measured vertical 
transverse motion, has its sensitive axis perpendic- 
ular to the axial motion of the table and measures 
motion up and down in the plane of the page. The 
second accelerometer, with a sensitive axis per- 
pendicular to the axial motion and perpendicular 
to the plane of the page measured the transverse 
mot on in the horizontal direction. A third accelerom- 
eter was mounted with its sensitive axis coaxial 
to the axial motion of the standard to measure the 
axial motion or axial acceleration. The three 
accelerometers were mounted on a small cube of 
stainless steel with two dummy accelerometers to 
balance the load of the vertical and _ horizontal 
accelerometers. 

The transverse motions indicated with these 
devices should not be attributed entirely to the 
standard since the accelerometers themselves have 
a transverse response. The transverse motions of 
the air-bearing exciter depend to a great extent on 
the bearing alinement and bearing clearance. 
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FicurE 4. Setup for measuring transverse motion with three 
piezoelectric accelerometers with their principle sensitive axes 
mutually perpendicular. 


Figure 5 shows the transverse motions of an air- 


bearing exciter using two experimental bearings 


with equal diameters but with closer radial clearances 
(less than 0.001 in.) than the bearings subsequently 


installed. Significant transverse motions in the 
vertical direction are observed in the regions of 
750 and 2,190 e/s, and in the horizontal direction 
in the regions of 400 and 2,190 ¢/s. The transverse 
motions near 400 and 750 c/s are attributed to the 
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Ficure 5. 


Ratios of transverse to axial acceleration using experimental air bearings. 
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Figure 6. Meas 


fundamental and probably the second bending mode 
resonances of the moving element. In figure 6 the 
measured transverse motions of exciter No. 676 
using air bearings (fig. 4) are compared with the 
transverse motions of exciter No. 619 using a flexure 
plate suspension and a tensioned wire suspension. 
Significant transverse motions of the air-bearing 
exciter in this case occur in the vertical direction at 
650 and 1,100 c/s and in the horizontal direction at 
730 and 1,050 c/s. The differences between the 
transverse motions plotted here (fig. 6) and the 
transverse motions in the experimental set of bear- 
fig. 5) are thought to be due to differences in 
bearing clearance and alinement. 


ines 


3. Calibration of the Air-Bearing Exciter 


The air-bearing exciter was calibrated by the 
reciprocity method. The results obtained by the 
use of the interferometer fringe disappearance and 
optical target methods are compared with the re- 
sults of the reciprocity method. 

The exciter rests on a 
from floor-borne 
springs. 


concrete 
vibrations by 


block isolated 
compression coil 


3.1. Reciprocity Calibration 
The reciprocity theorem for mechanical systems 
states that if a force is impressed at point A and the 
velocity is measured at point B, the same velocity 
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will be measured at point A if the force is impressed 
at point B. The reciprocity theorem as applied to 
electrodynamic vibration exciters is given in refer- 
ence [1]. In order for the theorem to hold, the 
system must be linear, i.e., the displacement, velocity 
and current must increase linearly with force and 
voltage. The application of reciprocity to this sys- 
tem is limited to frequencies above and below the 
resonances of the moving element. At frequencies 
near each resonance there is a loss in precision due 
to changes in the electromechanical properties of the 
calibrator in that the relationships are no longer 
linear. 

The reciprocity calibration consists of two experi- 
ments: (1) admittance 
between the driving coil and the velocity-sensing 
coil; (2) measurement of the voltage ratio of the open 
circuited velocity coil and driving coil when the 
exciter to be calibrated is driven by an external 
exciter. 


¢ 


measurement of transfer 


Transfer admittance measurements, Y,, were made 
with a series of weights attached one at a time to the 
mounting table of the exciter. Also a zero load 
transfer admittance measurement was made before 
and after each measurement with a weight attached. 
The zero load transfer admittance is denoted by ee 
The weights were increment weights ranging from 
0.1 Ib to 2.0 Ib. 





Using the measured values of Y, and Yo, graphs 


(1,2 «++) 
Y.- Y 20 
for each frequency, where 

increments of mass attached to the table. 
zero intercepts, J and iJ, of the resulting 
straight lines and their slopes, Q and 7Q, are computed 
by a weighted least square method. 

Having determined the voltage ratio P?, J, and Q, 
the sensitivity of the exciter is given by 


of the real and imaginary values of the ratio 


are plotted versus W 1.2 


Wa.2 


WQ ({RT7)"? v/in./see 


S—0.04288-+-0.04286 


where 


f is the frequency in cycles per second and 
jis the unit imaginary vector. 


3.2. Interferometer Calibration 


The air-bearing exciter was calibrated by the 
interferometric fringe disappearance method as given 
in reference [8]. Figure 7 shows the exciter and 
interferometer in position for calibration. The 
sensitivity of the exciter is given by 
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where £F is the rms voltage output of the velocity 
sensing coil corresponding to the first order fringe 
disappearance, s is the amplitude of displacement of 
the moving element (4.11 * 10~° in. for the first order 
disappearance using mercury light and green filter), 


f is the frequency, and S is given in peak to peak 


v/in./sec. The voltage was read by an a-c transfer 
standard using a VT'VM asa null. The data were 
taken with the exciter in a horizontal position. 
Readings for zero load and 1 |b load were taken over 
the range of 350 to 2,500 c/s. 

Each value of voltage, /, is obtained from the 
average of two observations, one for the first fringe 
disappearance during slowly increasing displacement 
amplitude and the second for the fringe disappear- 
ance during slowly decreasing amplitude. For the 
zero load, the mean deviation from the mean of the 
voltage values ranged from 0.1 to 0.7 percent. The 
deviations were computed from sets of six or more 
values at each test frequency over the range of 350 
to 2,500 c/s. The average deviation was 0.3 per- 
cent. The voltage values obtained by two observers 
differed from 0 to 0.8 percent for sets of three or 
more values per observer at each test frequency. 
The average repeatability was 0.3 percent. 

For the 1 lb load readings the reflecting glass was 
attached to the center of the weight. The mean 
deviation for the 1 lb load readings ranged from 0.1 
to 1.6 percent. The average deviation was 0.5 
percent. 
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FiGureE 7. 


Schematic arrangenent for calibration of exciter with interferometer. 
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The use of the interferometer fringe disappearance 
method is greatly facilitated on the air-bearing ex- 
citer as compared to the tensioned wire exciter 
because of the increased stability of the interference 
fringes. 


RECIPROCITY So» ® INTERFEROMETER Sow T 
Sv2e Sie 
Sie @ OPTICAL TARGET 
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3.3. Calibration With an Optical Target 
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FACTOR 


The air-bearing exciter was calibrated optically by 
small wire targets mounted on the exciter table. 
The targets were coated with particles from camphor 
smoke and the diameters measured to ten—thou- 
sandths of a centimeter. Two diameters of wires 
were used. A 0.2042 cm wire was used to calibrate 
at 15, 17, 20, 30, and 45 e/s. A 0.1048 ecm wire was 
used to calibrate at 55 and 70 c/s. The exciter was 
driven with increasing amplitude until the wire ap- a 
peared to vanish. The output of the velocity coil - ST... AE 
was used as the output corresponding to the ampli- se ey ia 
tude of vibration given by the diameter of the wire. 
Knowing the amplitude, frequency, and output of 
the velocity coil, the sensitivity of the exciter was 
determined. The average sensitivity for 13 observa- 
tions at 7 frequencies was 0.1859 v/in./see with a 
mean deviation of 0.0004. This compares with an 
average of 0.1857 v/in./see for the reciprocity calibra- 
tion from 10 to 400 c/s and 0.1857 volts/in./see for 
the interferometer calibration from 350 to 750 ¢/s. 

Figure 8 shows the calibration factors obtained 
for the air-bearing exciter by the three methods 
described above. 


CALIBRATION 


DEGREES 


PHASE LAG, 


FIGURE 8 ‘alibration factors of the exciter velocity se 


3.4. Calibration of Pickups 


After calibrations by the reciprocity method and 
the interference fringe disappearance method, vibra- 
tion standard 676 with air-bearing supports was 
used to calibrate a piezoelectric acceleration pickup 
and an electrodynamic velocity pickup. 

The piezoelectric pickup was calibrated by com- 
paring its voltage output with the voltage generated 
by the velocity coil of the exciter. The sensitivity 
of the pickup is given by 


Figure 9. Calibration factors of a piezoelectric acceleration 
pickup. 


voltage output of pickup 


exciter 


voltage output of velocity coil” 


where the S.cier IS corrected to read v/g instead of 
v/in./see (g=386 in./sec’). This ratio is determined 
by the following transfer admittance ratios: 


oa/V. 


pickup 


The electrodynamic velocity pickup was calibrated ee eee ee ee 
similarly. The results of these calibrations are 
given in figures 9 and 10. 


Figure 10. Calibration factors of an electrodynamic velocity 
| pickup. 
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4. Discussion 


Although the air-bearing suspension resulted in 
eliminating much of the transverse motion for ex- 
citer 676 (see fig. 6), there are still frequencies near 
700 and 1,000 c/s which should be avoided during 
calibrations of vibration pickups. 

It is extremely sta that five of the six de- 
grees of freedom of the mounting table can be 
eliminated due to the geometrical configuration of 
the moving element as shown schematically in figure 
1. Also, there are several materials with different 
moduli of elasticity incorporated into the moving 
element. Due to these factors, axial and flexural 
resonances may be excited in the moving element 
that interfere with its use as a pickup calibrator. 

It is attractive to consider the possibility of con- 
structing moving elements of uniform cross section 
and maximum axial and transverse  stiffnesses. 
Such elements should be much better for use in 
vibration standards. 

Certain advantages accrue from the substitution 
of the air-bearing suspension guidance system for 
flexure plates and tensioned wires. The air bearings 
allow larger axial motions than were possible with 
the mechanical suspension systems. The flexure 
plates and wires themselves were excited in resonance 
at certain frequencies which interfered with the 
utility of the exciter. With air bearings, the moving 
element is free to move axially but is prevented from 
moving laterally; this minimizes transverse motions 
at the mounting table. 


Richard Harwell machined the bearings and sup- 
porting assembly. Ruth M. Woolley performed 
most of the computations. James A. Miller and 
L.. R. Sweetman rendered valuable assistance and 
constructive criticism. 
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Leak-Resistant Rotation Seal for Vacuum Applications 


Frank L. Howard 


(June 6, 1963) 


\ leak-resistant rotation seal for application with vacuum systems is described. 


The 


packing comprises a system of Teflon cones with deformable ridges that act as bearing 


surfaces. 


The seal described herein was designed to convert 
a rotary-tube stirrer used for measurement of freezing 
point at 1 atm to one that would allow measurements 
at reduced pressure, approaching triple point con- 
ditions. 
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Figure 1. Sectional view of seal. 


A. Thermometer well used in this application. 
through a standard packing gland. 

B. Housing for the seal, made of stainless Steel, 0.d, =2 in. 

C, Connection to vacuum system. 

D. Teflon spacer. 

E, F. Packing shown in detail in figure 2. 

G, Compression ring, tightened by six $42 machine screws passing through H, 

I, Connection for guard vacuum, 

J. Rotating tube with polished surface on bearing area. 


The thermometer was attached 
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For this application it was desired that the packing 
material be inert, provide low friction and effectively 
prevent leakage. Polytetrafluoroethylene (Teflon 
brand) was selected for this purpose, and the packing 
sections were designed to partially obviate the cold 
flow properties of this material and to reduce the 


DETAIL OF PACKING 


0.D. 
































NOTE: ALL SLOTS 
0.025 DEEP 
0.040 WIDE 
0.020 BETWEEN SLOTS 


Figure 2. Enlarged view of packing. 


E;. Inner cones of Teflon, the internal diameter is 0.010 in. smaller than outside 
diameter of J (fig 

E2. Outer cones of teflon, similarly larger than inside diameter cf housing B 
(fig. 1). 

E3. Brass cone, a taper of 14° was used. 

F. Guard vacuum collecting ring. 





friction inherent in large bearing surfaces. The parts 
of the seal are illustrated in figures 1 and 2. When 
the plastic packing cones are compressed vertically 
the small brass cones separating them cause them 
to tighten laterally against the rotating shaft and 
the outer housing. Some deformation of the pro- 
tuberances between the grooves enables a tight fit. 

The seal was tested with a helium mass spectrom- 
eter leak detector, while the shaft was rotated at 
20 rpm. It was found that at about 5x<107~ mm Hg 
the leakage rate was less than 210 
gas per minute. At this pressure, the guard vacuum 


was disconnected and no increase in pressure or in 
rate of leakage was noted. 
the temperature of the housins 


During 2 hr of operation 
rincreased about 10 °C. 


- em? (STP) of 


The assistance of Herbert B. Lowey and Tully J. 
Pound, in making this apparatus is gratefully 
acknowledged. James L. Cross suggested use of 
the guard vacuum. 
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Various lots of tungsten and rhenium wire were obtained from leading American manu- 
facturers. Eleven tungsten-rhenium thermocouples were made up from these wire samples 
and tested at approximately 100 deg C intervals up to 2000 °C. A fourth order equation 
and a third order equation were selected to represent the emf of the tungsten-rhenium 
thermocouple in the range from 0 to 1000 °C and from 1000 to 2000 °C respectively. Using 
these two equations, reference tables were established and are presented in 5 deg C intervals 
from 0 to 2000 °C and in 10 deg F intervals from 32 to 3640°F. Inverse tables giving temper- 
ature in °C and °F at 20 microvolt intervals are also included. Information is furnished on 
variations due to wire size and among different wire lots and different wire manufacturers. 
Spectrochemical analyses of the tungsten and rhenium elements that were used to represent 


the tables of temperature versus emf are listed. 


A graphic comparison is made between the 


NBS emf values and values from other investigators. 


1. Introduction 


In 1931 Goedicke [1] ' investigated the possibility 
of using tungsten and rhenium as high temperature 
thermocouple elements. However, this early work 
was quite limited in scope and did not, for example, 
result in a precise relationship between temperature 
and emf for these thermocouple elements. 

Between 1952 and 1956 the Battelle Memorial 
Institute conducted the first detailed investigation 


of various properties of rhenium [2, 3] and in 1957-58 
their studies included the thermoelectric properties of 


rhenium, tungsten, and other high temperature 
thermocouple elements [4, 5]. In this latter work a 
power series Was used to express the temperature-emf 
relationship of the tungsten-rhenium thermocouple 
between room temperature and 2200 °C. The 
workers at Battelle concluded that rhenium had 
several distinct advantages as a high temperature 
thermocouple element e.g., it was found to have 
excellent strength from room temperature to 2000 °C. 
was highly resistant to the water cycle and retained 
it ductility after thermal cycling above its recrystalli- 
zacion * temperature. Unfortunately, tungsten is 
affected by the water cycle and in addition becomes 
quite brittle upon reaching its recrystallization 
temperature. The recrystallization temperature of 
tungsten is not well defined. Smithells [6] reports 
that recrystallization can occur at any temperature 
from 1000 to 2000 °C depending upon the metal- 
lurgical properties of the tungsten, and Hall and 
Sikora [7] report partial recrystallization of commer- 
cially pure sintered tungsten at 1715 °C and complete 
recrystallization at 2090 °C. The brittleness of 
tungsten due to recrystallization can be somewhat 
inhibited by a process called “doping.” This process 
involves the addition to the tungsten of small 


1 Figures in brackets indicate the literature references at the end of this paper. 
2 Recrystallization is defined as the formation of a new, strain-free grain struc- 
ture from that existing in the cold worked metal, usually accomplished by heating. 
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quantities of materials such as compounds of potas- 
sium, silicon, or aluminum. These “doping’’ com- 
pounds reduce grain growth and essentially raise the 
recrystallization temperature of the tungsten. In 
the case of pure tungsten, however, very little can 
be done to inhibit the brittleness resulting from 
recrystallization. 

The water cycle effect which occurs at elevated 
temperatures alters the chemical and metallurgical 
properties of tungsten. In the case of tungsten 
thermocouple elements, this effect will eventually 
erode the metal to a point where the thermocouple 
fails. Langmuir [8] first demonstrated the water 
cycle effect with a tungsten filament lamp. A small 
amount of water vapor was released in a vacuum 
lamp and the filament was heated to incandescence. 
The water vapor coming into contact with the fila- 
ment was decomposed, the oxygen combining with 
the tungsten and the hydrogen being evolved. The 
oxide distilled to the glass envelope where it was 
reduced to metallic tungsten by the atomic hydrogen 
given off by the filament. Reduction of the oxide 
resulted in the reconstitution of water vapor and the 
‘water cycle” was able to repeat itself indefinitely. 
Thus, in applications where tungsten is used in a 
vacuum as a thermocouple element, care must be 
taken to keep the area surrounding the thermocouple 
free from water vapor at elevated temperatures. 

Additional investigations concerning the tungsten- 
rhenium thermocouple were undertaken by Lachman 
and Kuether [9-11] and included useful information 
such as chemical analysis of the thermoelements, emf 
reproducibility of the thermocouple after prolonged 
temperature cycling and a table of temperature 
versus emf from 50 to 4000 °F. 

The above workers concluded that the tungsten- 
rhenium thermocouple has the following advantages: 

1. High thermoelectric potential. 

2. High thermoelectric power. 

3. Very high melting point of both elements of the 

thermocouple. 





4. Chemical stability in vacuum, inert or reducing 
atmospheres up to 2200 °C. 
5. Calibration accuracy of 0.1 mv up to 2200 °C, 
6. Low cost of the tungsten element. 
7. Good reproducibility after thermal cycling at 
2200 °C 
Some of the disadvantages of the thermocouple are 
1. Brittleness of the tungsten element after re- 
crystallization. 
Susceptibility of the tungsten element to de- 
structive erosion due to the water cvcle effect. 
3. High cost of the rhenium element. 
4. It cannot be used in an oxidizing atmosphere. 
In view of the favorable advantages of the tung- 
sten-rhenium thermocouple as a high temperature 
measuring instrument, a program was initiated at the 
National Bureau of Standards to carry out studies on 
this thermocouple. Although various studies were 
conducted on the tungsten-rhenium thermocouple by 
the authors mentioned, it was felt that an effort 
should be made at NBS to further evaluate this 
thermocouple using thermoelements that represent a 
variety of wire lots from various manufacturers. 
The results of this study could then be compared to 
the earlier reported results in order to ascertain the 
thermoelectric reproducibility of these materials. 


2. Thermocouple Materials 


2.1. Tungsten Elements 


The tungsten elements that were used in this 
study were obtained from three major American 
tungsten wire manufacturers and represent a total 
of eight lots of wire. Although this tungsten wire 
was manufactured primarily for electronic applica- 
tions, the chemical purity of the wire as determined 
by spectrochemical analysis see table 1) was suffi- 
ciently high as to render it acceptable for this par- 
ticular work. All of the tungsten elements consisted 


250 


of two types of tungsten. One type of tungsten re- 
ferred to by the manufacturers as “black as drawn” 
tungsten comprises five of the eight lots i.e., lots 
L-2, M-2, M-3, M-4, and N-1. The letters BL 
have been placed in parenthesis after each of these 
lot designations to denote the “‘black as drawn’’ wire 
type. The second type of tungsten is referred to as 
“chemically cleaned” tungsten and comprises three 
lots i.e., lots L-1, M-1, and N-2. The letters CL 
are used to denote this type of wire. According to 
the manufacturers, the ‘black as drawn’ wire re- 
ceives no additional treatment or processing after it 
has been drawn to the desired size. On the other 
hand, the chemically cleaned wire after having been 
drawn to the desired size is immersed into a hot 
caustic solution (usually KOH or NaOH), is rinsed 
and then dried. All of the tungsten elements con- 
tained ‘doping’? compounds as reported by the 
manufacturers. 

The prefixes L, M, and N designate the three 
manufacturers from whom the tungsten lots were 
obtained. One element was selected from each of 
the eight tungsten lots and labeled the same as the 
lot number followed by the letter “A.” Five of these 
eight tungsten elements were used only once to 
represent the positive (+-) leg of a particular tung- 
sten-rhenium thermocouple and the other three 
elements viz, L-1(CL)A, M-4(BL)A and N-2(CL)A 
were used twice as a positive leg. For example, the 
element L—-1(CL)A was paired with the two rhenium 
elements 219A and 228A, the combinations being 
designated thermocouple No. 1 and thermocouple 
No. 9 (see table 3). This arrangement resulted in a 
total of 11 thermocouples. This method of thermo- 
couple selection was chosen rather than a method 
whereby all possible thermocouple combinations 
(each tungsten element rhenium ele- 
ment) could be represented. All thermocouple 
combinations could be equally represented by 
averaging the data in figures 1 and 2, and then 
combining the averages via measurements of a 
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FIGURE 1, 


Emf differences between eight tungsten elements. 


[With N-2(CL)A as the reference element] 
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FIGURE 2. 


Emf differences between eight rhenium elements. 


[With 228B as the reference element] 


single tungsten-rhenium combination (see also in 
Experimental Procedure). However, in order to 
use this method effectively, the following must be 
realized. 

(a) The tungsten and rhenium elements that are 
used as reference elements must be reproducible. 
Since all of the elements could not be tested in one 
furnace run, it was necessary to heat the reference 
elements several times. The tungsten and rhenium 
elements that were tested gave indications of not 
being reproducible from one furnace run to the next 
(Experimental Procedure). 

(b) A considerable amount of weight would be 
placed on the single tungsten-rhenium combination 
measurement. If large errors were encountered in 
this measurement (poor optical pyrometer readings, 
for example), these errors would be grossly reflected 
in the final temperature—emf relationship. 

Each of the eight tungsten elements was spectro- 
chemically analyzed. A listing of the magnitude of 
the impurities is given in table 1. 


2.2. Rhenium Elements 


All of the rhenium elements associated with this 
study were obtained from the Chase Brass and 
Copper Company which is the sole American 
producer of rhenium wire in commercial quantities. 
Eight rhenium elements were paired with various 
tungsten elements listed in table 3. These 
rhenium elements were taken from four lots of wire, 
numerically designated 219, 221, 226, and 228, with 
each element labeled A, B, or C of a particular 
lot *. All eight elements were 0.020 in. in diameter. 
As in the case of the tungsten elements, five of the 
eight rhenium elements were used only once to rep- 
resent the negative leg of a particular tungsten- 
rhenium thermocouple and the three elements, 219A, 
226C, and 228B, were used twice. 


as 


Elements A, B, and C ofa particular lot were remote from each other as far as 
the position of each element on a spool of wire is concerned. 
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TABLE 1. Spectrochemical analysis of tungsten elements 


Tungsten element 
Chemical 
element 
M-2 M-3 
BL)A | (BL)A 


0.1% 

0.01% 
0.01°%-0.001% 
0.001°%-0.0001°%. 

E = Less than 0.0001°% 
Not detected. 


A=1.0% 
B=0.1% 


The eight rhenium elements were spectrochemically 
analyzed by the same laboratory‘ that analyzed the 
tungsten elements and the results of the analysis are 
listed in table 2. 


4 Spectrochemistry Section of the National Bureau of Standards. 


"a 





TABLI & Spectro hen ical analysis of rhenium ele ments 


Rhenium elen 


3. Apparatus 


All of the thermocouple and thermoelement tests 
between 0 and 2000°C were performed in one 
furnace. This furnace consists of a tantalum tube 
that serves as a heating element. The tantalum 
tube furnace is tully described elsewhere. [12] The 
thermocouples and thermoelements were vertically 
suspended inside the tantalum tube with the 
measuring junctions located in a molybdenum black- 
body which was also suspended in the central region 
of the tantalum tube. All vests were conducted in a 
purified helium atmosphere, the pressure of which 
~ page from 360 mm Hg at 20 °C to about 660 mm 
Hg at 2000 °C. The helium was purified by passing 
it Sats 9 a titanium-zirconium purifying apparatus 
[12] before being released into the cir ace chamber. 

In the range trom 100 to 1000 °C the instrument 
used for temperature reread was a Pt 1 
percent Rh versus Pt 30 percent Rh thermocouple. 
A modified commercial optical ——— was used 
in the 1000 to 2000 °C range. The percent-30 
percent” thermocouple was calibr: ceil by direct 
comparison to a standard platinum versus platinum 
10 percent rhodium thermocouple. The latter 
thermocouple received a fixed point calibration as 
described in the literature [13]. An estimate of the 
maximum uncertainty in the comparison calibration 


| ohm standard resistor in series with the lamp. 


! 


is +0.5 °C. The optical ;pyrometer was calibrated 
and used on a basis of temperature versus optical 
pyrometer lamp current. The calibration was 
performed by sighting the pyrometer on a tungsten 
ribbon filament lamp, the temperature of which was 
determined by the Fairchild optical pyrometer ° 
at NBS. An estimatecf the w2amum uncertainties 
in the optical pyrometer calibration is +3.0° at 
1000 °C and + 6.0° at 2000°C. The lamp current was 
determined by measuring the voltage drop across a 1- 
The 
thermocouple emf measurements were made with a 
L & N K-3 type potentiometer having a limit of 
error of about 0.7 uv at 1000 vv and 10 pv at 50,000 pv. 


4. Experimental Procedure 


In the tests conducted in the 0 the 1000 °C range, 
the Pt 1 percent Rh versus Pt 30 percent Rh thermo- 
couple and the tungsten and rhenium elements 
were placed in high purity aluminum oxide insulat- 
ing tubes with the measuring junctions exposed at 
one end. The measuring junction of the 1 percent- 
30 percent thermocouple was mechanically and 
thermally bound to the measuring junction of the 
tungsten and rhenium elements by wrapping with 
platinum wire. Generally, a total six thermo- 
elements (two tungsten, two rhenium, and the two 
platinum-rhodium elements) were placed in the 
furnace at one time. The emfs of thermocouples 
No. 1 through No. 11 (table 3) were measured at 
approximately 100 deg intervals from 100 to 1000 °C 
with the reference junctions maintained at 0 °C. 
In addition to these measured values, an assumed emf 
of 0.000 my at 0 for each thermocouple gave a 
total of 121 data points for the 0 to 1000 °C range. 
At each of the calibration points, the temperature 
in the hot zone of the furnace as determined by the 
1 percent-30 percent thermocouple was allowed to 
stabilize before the emfs of the tungsten-rhenium 
thermocouples were measured. 
~$ The NBS Fairchil d optical pyrometer, also refer 
s “The NBS Optical Pyrometer” is a hich precision 


C, O. Fairchild of NBS and is used at NBS for testir 
und pyrometer lamps submitted for calibratior 


red to in the literature [14] 
instrument designed by 
other optical pyrometers 


TABLE 3. Thermoco ple elements 


All elements were 0.020 in. in diat 


neter.) 


Rhenium 
elements 


Tungsten elements 


Thermocouple 
number 
Element 
number 


Manufac- 
turer 





The eleven tungsten-rhenium thermocouples were 
also tested at approximately 100 deg intervals 
between 1000 and 2000 °C with the reference junc- 
tions at 0 °C. At each of the 100 deg intervals the 
emf of the thermocouple was read three times and 
an optical pyrometer observation was made three 
times but in an alternating sequence. This resulted 
in a total of approximately 360 temperature-emf 
determinations. The data that were recorded for 
each thermocouple at each of the 100 deg intervals 
were averaged. This resulted in 101 averaged 
temperature-emf data points in the range from 1000 
to 2000 °C. In most instances, and particularly 
at temperatures above 1700 °C, it was necessary to 
anneal the thermocouples at a slighty higher tem- 
perature than the temperature at which the data 
were taken. For example, in order to render the 
thermocouples stable at 1800 °C, it was necessary 
to raise the temperature of the furnace to about 
1825 °C for a 5 to 10-min period and then lower the 
temperature to the initial 1800 °C. During the 
time that the thermocouple emf and optical pyrom- 
eter readings were recorded, the temperature in 
the hot zone of the furnace was not fluctuating by 
more than 0.05 percent per minute and in most cases 
the fluctuations were considerably less than this 
amount. 

Since the tantalum tube furnace design allows 
the thermocouple to be freely suspended vertically, 
it was possible to avoid the use of ceramic supports, 
thereby eliminating two possible effects: (a) chemical 
contamination, and (b) electrical errors due to 
conduction in the ceramic material 

A d-c electric arc operating in a helium atmosphere 
was used to fuse the measuring junction of the 
tungsten and rhenium elements. 

All of the emf measurements related to the tung- 
sten and rhenium elements listed in table 3 were 
conducted in four separate test series. These four 
test series, in chronological order were as follows: 
(1) The emf’s of the eleven thermocouples were 
measured between 1000 and 2000 °C, (2) the emf’s 
of the eleven thermocouples were measured between 
0 and 1000 °C, (3) the emf differences between the 
tungsten elements and between the rhenium ele- 
ments were measured between 0 and 1000 °C, and 
(4) the emf differences between the tungsten ele- 
ments and between the rhenium elements were 
measured between 1000 and 2000 °C. All readings 
were made with increasing temperature. After the 
first test series was completed, it was necessary to 
remove from 6 to 8 in. of wire from the measuring 
junction end of the tungsten elements since this 
portion of the wire had become brittle upon heating 
to 2000 °C. This removal of the embrittled portion 
of each tungsten element was necessary in order to 
reweld the elements for the succeeding parts of the 
test series. 

Since short portions of the tungsten elements were 
removed after the first test series, it is apparent that 


6In a few cases four thermocouple emf readings and four optical pyrometer 
observations were made in alternating sequence. 
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the emf’s that were measured in the second, third, 
and fourth test series would in some cases be incon- 
sistent with those measured in the first series. If 
the emf’s of the eleven thermocouples tested in the 
first and second series are compared with the calcu- 
lated? emf’s of the same thermocouples as derived 
from the differences measured between elements in 
the third and fourth series, it can be seen that the 
respective values are not consistent in all cases. 
These inconsistencies can be attributed to (a) the 
lack of emf reproducibility of the tungsten elements 
before and after the brittle portions ef the elements 
were removed, (b) the lack of emf reproducibility of 
the rhenium elements, and (c) experimental errors 
encountered in each of the test series. 

The long-term stability of the tungsten-rhenium 
thermocouples was not determined in this study. 
However, four test thermocouples were made up 
from the tungsten and rhenium wire lots pertinent 
to this study and were thermally cycled between 
1100 and 2000 °C. The elements of these four 
thermocouples were none of the elements listed in 
table 3. The thermocouples were initially heated to 
1100 °C (approx.) and their emf’s were measured at 
that temperature. Then the thermocouples were 
heated to 2000 °C (approx.), the emf’s again being 
measured, and then the temperature was lowered to 
1100 °C. This heating cycle was repeated twice. 
The third and final measurement at 1100 °C showed 
that the emf’s of the four thermocouples had in- 
creased (as compared to the initial measurement at 
1100 °C) between 60 and 122 uv. This is equivalent 
to 3.5 and 7.2 °C, respectively. The total heating 
time during this cycling test was approximately 
5 hr. 

None of the tungsten and rhenium elements listed 
in table 3 were heat treated prior to testing except in 
the test series involving emf measurements between 
1000 and 2000 °C. In this case, of course, the ele- 
ments were heated to 1000 °C before emf measure- 
ments were begun. 


5. Computations 


In order to represent the raw data that were 


| obtained in this study in a useful and concise form, 
several electronic computer programs were utilized. 


This data handling was performed with an IBM 
7090 computer. 

One of the first computation steps was to derive 
one or more polynomial equations that would repre- 
sent the thermocouple data over the entire tempera- 
turerange. Forinvestigative purposes, the computer 
was programmed to derive by the least squares 
method (a) a single polynomial equation to represent 
the 0 to 2000 °C range and (b) one polynomial 
equation to represent the 0 to 1000 °C range and 
another polynomial equation to represent the 1000 
to 2000 °C range. The function fitting was applied 

7 The emf was measured between the elements N-2 (CL) A and 228B as thermo- 
couple No. 11 and also between each of these elements and the other 


and 7 rhenium elements. 
emf’s can be determined. 


‘ 7 tungsten 
From this relationship the other 10 thermocouple 





to all of the thermocouple data in each of the tem- 
perature ranges, i.e., 121 data points in the low 
temperature range and 101 data points in the high 
temperature range. ‘The results of this investigation 
showed that a sixth order equation would be needed 
to adequately represent the thermocouple data over 
the 0 to 2000 °C range. The same conclusion was 
reached by Sims, Gaines, et al. [5] over the tempera- 
ture range of investigation. In the (b) part of the 
investigation it was determined that a fourth order 
equation and a third order equation would ade- 
quately represent the 0 to 1000 °C and the 1000 to 
2000 °C ranges respectively. In practical applica- 
tions, considerable time and effort is saved by using 
two equations. Hence, this method was selected to 
represent the temperature-emf relationship of the 
tungsten-rhenium thermocouple. 

In examining these two equations it was found 
that the derivative dE/dt was not the same for both 
equations at 1000 °C and that the emfs at 1000 °C 
differed by about 19 uv. A computer operation was 
used to adjust the coefficients of the fourth order 
equation (for the 0 to 1000 °C range) such that the 
emf and slope of the two equations at 1000 °C were 
equal. 

The equations for the average emf F in millivolts 
of the tungsten-rhenium thermocouples at a tem- 
perature ¢(°C) are 
for the range 0 to 1000 °C 


0.62893850 * 107-°4+ 0.20717363 
0.150672S80 * 107-7 +0.: 


Ts ig 
37778323 X 107"'t4 


for the range 1000 to 2000 °C 


9 9 


1() 1871033 107% 
52 < 1075#?—0.17634201 


oo. 5162 
0210675 or. 

The emf deviations of each of the eleven tungsten 
rhenium thermocouples from the above equations 


4 ——— - - 





are shown graphically in figure 5 for the range 0 to 
1000 °C and in figure 6 for the range 1000 to 2000 °C, 
Sd? 
—N 
where d is the emf deviation of each data point from 
the respective polynomial equation, n is the number 
of data points * and NV is the number of coefficients in 
the polynomial equation was determined for the two 
temperature ranges. For the 0 to 1000 °C range the 
standard deviation was 74.1 wv and for the 1000 to 
2000 °C range it was 70.1 uv. An estimation of the 
maximum uncertainties in the temperature measure- 
ments made between 0 and 1000 °C is +1.0 °C and 
between 1000 and 2000 °C is +0.35 percent. 


6. Results 


The standard deviation of the form S \ 
n 


Tables for the tungsten-rhenium thermocouple are 
presented giving emf at 5 °C and 10 °F intervals 
(tables 5 and 7). Inverse tables giving temperature 
in °C and °F at 20 uv intervals are also included 
tables 4 and 6). 

The emf differences were measured between each 
of the eight tungsten elements. N-2(CL)A was 
arbitrarily chosen as a reference element and the 
emfs developed between it and the other seven 
tungsten elements were measured between 0 and 
2000 °C (fig. 1).° 

Two of the eight tungsten lots were selected for 
studying the emf differences between elements of 
various diameters. The two lots selected were 
N-2(CL) and M-3(BL). From each lot three 
elements were obtained and designated with the 
suffixes B, C, and D. These were 0.020, 0.024, and 
0.028 in. in diameter (nominal) — respectively. 

2(CL)A™ was again chosen as a reference element 


* In the low temperature range 
ture range 


, n Was 121 data points an 
n Was the 101 averaged data points 
The reference element, N-2(CL)A was connected to the 
potentiometer 
N-2(CL)A and N-2(CL)B were adjacent 
» for the M-3(BL)A and 


lin the high tempera- 
negative post of the 


ts cut from one spool of 
.) B elements 


Ee 


w- 2 A 


ty 














i300 1400 


1500 


"| 


1600 1800 1900 2000 


TEMPERATURE, °C 


FIGCRE 3. 


With N-2(CL 


Emf differences between tungsten elements of various diameters from Lot N 


2(CL) 


A as the reference element] 


342 














l 


M-3 (BL)O 
M-3 GLC 0 .028"diam 
0.024"diam 


M-3 (BL)A 
0.020"diam 


M-3 (6L)B 
0.020"diam 








1300 1400 


1 
1000 iitete) 


1200 


1500 


it J 
1600 1700 1800 


TEMPERATURE, °C 


FIGURE 4 
{With M-3(BL 


and the emf’s developed between it and elements 
N-2(CL)B, C, and D in the temperature range from 
1000 to 2000 °C are shown graphically in figure 3. 
These measurements give an indication of the emf 
differences that can occur between two adjacent 
elements and also the differences that can occur 
between elements from the same lot but of varying 
diameters. 

Similar measurements were conducted with the 
elements M-—3(BL)A, B, C, and D, using M-3(BL)A 
as a reference element. Emf differences =m 
these elements are shown in figure 4. All of the 
measurements related to the tungsten © of 

varying diameters were made after the four series 
of tests previously mentioned. 

The emf differences between the eight rhenium 
elements were measured in the range 0 to 2000 °C 
by using 228B as a reference element." 
differences are shown in figure 2. 

Since rhenium elements of diameters other than 
0.020 in. were not available from the manufacturer 
for the four lots 219, 221, 226, and 228, no tests 
were conducted to denote emf differences between 
rhenium elements of various diameters. 

The thermoelectric power of the tungsten-rhenium 
thermocouples included in this study average about 
17.6, 17.6, 13.1, and 6.0 uv per deg C at 500, 1000, 
1500, and 2000 °C, respectively. A maximum 
thermoelectric power of 18.33 uv per deg C occurs at 
715 °C 


These 


7. Discussion 


In general it can be concluded that the total emf 
spread of the eleven tungsten-rhenium thermo- 
couples included in this study is not exceedingly 


1 The reference 
potentiometer 


element 228B was connected to the negative post of the 
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Emf differences between tungsten elements of various diameters from Lot M-3 (BL 
A as the reference element] 


large if one considers that the thermocouple elements 
represent four wire manufacturers and twelve lots of 
wire. The greatest emf spread of the eleven thermo- 
couples in the range 0 to 1000 °C was 246 uv at 
70 °C (fig. 5). This corresponds to approximately 
15.2 °C. Likewise, in the range 1000 to 2000 °C the 
greatest emf spread was 263 uv at 2000 °C corre- 
sponding to approximately 44.0 °C (fig. 6). The 
large emf spread of some of the thermocouples at 
370 °C was rather surprising. It is interesting to 
note that all of the thermocouples (except No. 10) 
deviate a considerable amount from the table values 
in the 0 to 700 °C range. In the range 1000 to 
2000 °C the deviations of all of the thermocouples 
from the table values are quite random. 

Figures 1 and 2 definitely show that the large emf 
deviations of some of the thermocouples from the 
table values in the 0 to 700 °C range are, for the 
most part, due to the tungsten elements. The 
maximum emf spread of the eight tungsten elements 
at 400 °C is 237 wv and the maximum spread of the 
eight rhenium elements at that same temperature is 
22 uv. The greatest emf spread between any two 
tungsten elements and any two rhenium elements 
occurred at 2000 °C in both cases. These maximum 
emf spreads at 2000 °C are 400 uv for the tungsten 
elements and 205 uv for the rhenium elements. It 
should be pointed out that the maximum emf 
spread of all of the tungsten elements should not 
be directly compared to that of all of the rhenium 
elements since the former represents three different 
manufacturers and the latter represents only one 
manufacturer. If the maximum emf spread of the 
elements produced by each of the four manufacturers 
are intercompared, it can be seen that these maxi- 
mums are not significantly different i.e., the maxi- 
mum spreads are all between 105 and 205 uv. 

The tests that were conducted to determine the 
emf differences between tungsten elements of various 
diameters were limited in scope but gave a genere! 
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Ficure 6. Deviations of thermocouples No. 1 through 
indication of the magnitude of the differences that 
canoccur. These tests show that the emf differences 
between tungsten elements of various diameters 
from one particular lot can be as large as the emf 
differences between elements from different lots and 
different manufacturers. At temperatures above 
1800 °C the two elements N-2(CL) A and N—-2(CL)B, 
which were adjacent samples on a spool of wire, 
developed a fairly large emf and at 2000 °C the 
emf was larger than that developed between the 
reference element and elements “C” and “D” 
(fig. 3). This same thing occurred between the two 
adjacent elements M-3(BL)A and M-3(BL)B at 


No. 11 from the table values in the range 1000 to 2000 °C. 


temperatures near 2000 °C (fig. 4). At temperatures 
between 1000 and 1800 °C the emf developed 
between elements various diameters remained 
relatively small. 

In general, the results of the spectrochemical 
analysis of the thermocouple elements show that the 
tungsten elements were of a higher purity than the 
rhenium elements. Of the eight tungsten elements 
that were analyzed, the element M-3(BL)A showed 
the greatest amount of impurities and M-2(BL)A 
showed the least impurities. If an average value is 
taken for the percentage range associated with each 
impurity listed (table 1), the total percent of im- 


of 
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purities for M-3(BL)A would be 0.0665 percent. 
Likewise, the total impurities for M-2(BL)A would 
be 0.0060 percent. For comparison purposes, the 
impurities in reference grade thermocouple platinum 
reportedly {14, 15] are of the order of 0.001 percent. 
Thus, the purity of M-2(BL)A is approaching that 
of reference platinum and M-—3(BL)A con- 
tains over 60 times the total impurities in reference 
platinum. 

Of the eight rhenium elements that were analyzed, 
221A showed the largest percentage of impurities 
with 1.012 percent and element 228B showed the 
least impurities with 0.0215 percent (table 2). Thus, 
element 221A contains more than a thousand times 
the impurities of reference grade platinum. 

In view of the above analysis, it can be concluded 
that although termocouple elements of the same type 
are relatively pure, they may be appreciably different 
thermoelectrically. The opposite behavior, thermo- 
electric similarity in spite of differences of composi- 
tion is shown by the two rhenium elements 221A and 
228B. The emf developed between these two ele- 
ments from 0 to 1900 °C was less than 20 uv (fig. 2) 
and vet element 221A contained nearly 50 times 
more impurities than element 228B. These results 
strongly indicate that rhenium possesses a thermo- 
electric uniqueness which manifests itself by the 
relatively small emfs produced between rhenium 
elements that contain significantly different quanti- 
ties of impurities. 

The emf differences that were recorded between 
the reference tungsten element N-2(CL)A and the 
other seven tungsten elements (fig. show no un- 
chemically cleaned”’ 
black as drawn” type of wire. 
a graphic presentation of the emf 
the for the tungsten- 
rhenium thermocouple as given in this paper and the 
values reported by et al., [5] and 
Lachman [11]. The values from the three sources 
agree reasonably well between 1000 and 2000 °C but 
in the range 0 to 1000 °C the differences between the 
NBS values and the values from Sims, Gaines, et al., 
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usual emf characteristics of the “ 
wire as opposed to the * 

Figure 7 is 
differences 
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[5] are quite large. The largest emf difference over 
the entire range occurs at 340 °C with a 415 uv 
difference between the NBS values and the values 
from Sims, Gaines, et al., [5]. This corresponds to 
approximately 26.4 °C. The greater part of the 
differences that occur between the values from the 
three investigators is probably due to chemical and/or 
metallurgical differences between the thermocouple 
elements rather than errors arising from temperature 
measurement techniques. 

The results of the various tests that were conducted 
in this study indicate that the thermoelectric quality 
of commercially available tungsten and rhenium 
wire could be improved. Causes of the large emf 
differences that may exist between tungsten elements 
from different lots and from different manufacturers 
might be identified by studying the individual effects 
of chemical and metallurgical variables on the 
thermal emf. Identification of the important vari- 
ables could then lead to their control during the 
manufacturing process. Although some of the 
rhenium elements contained relatively large per- 
centages of impurities, it was apparent that these 
impurities did not cause as large emf differences as 
in the tungsten elements. If the impurities in future 
lots of rhenium wire can be considerably reduced, 
perhaps the emf spread between elements would 
correspondingly be reduced. 

It should be emphasized that the various tests 
indicate that the emf produced by a particular 
tungsten-rhenium thermocouple at a specific tem- 
perature is highly dependent on the lot of wire from 
which the thermocouple was fabricated and the 
degree of heat treatment (annealing) the thermo- 
couple has received. In view of these factors, cau- 
tion should be exercised in using the table values cited 
herein as the correct values for a specific tungsten- 
rhenium thermocouple. 

It is anticipated that reference tables for the 
tungsten 3 percent rhenium versus tungsten 25 per- 
cent rhenium thermocouple and/or the tungsten 5 
percent rhenium versus tungsten 25 percent rhenium 
thermocouple will be prepared at NBS in the near 
future. 
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TABLE 6. Tungsten versus rhenium thermocouples 
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Standard X-ray diffraction powder patterns, H. E. Swanson, 
M. C. Morris, R. P. Stinchfield, and E. H. Evans, NBS Mono. 
25, Section 2 (May 8, 1968), 34 cents. 

Standard X-ray diffraction powder: patterns are presented 
for the following thirty-seven substances: Al (PO ;)3, SbF3*, 
Ba;(AsQ,).*, Ba(ClO,4)..3H,O, Cd(CN).*, CAWO,, Cs.0sBro*, 
Cs. sCly*, a-CrPQ,, ColHg(CNS$),]*, B-Cost 4, Dy;( ig 
(GaQO,)3*, ErMnO,;*, EusGa, (GaQ,);*, Gd;Gao(GaQ,)5*, 
Li;AsO,*, LisPs0,.3H,0*, Li,WO,.4H,0*, Lu;Ga.(Ga0,)s*, 
LuMnQ;*, MnWO, (huebnerite), HgF.*, NiSO,, NiW0O,*, 
K,ReCk, K,RuCl,NO*, RbCIO,*, RbIO *, AgsSeO,*, 
NaCNO*, Na,WO,.2H,0, 8-NayPyO».4H,0,  Sr3(AsO,)o*, 
ThAsOy*, TICLO,*, YAsO,, ZnWO,*. Eleven are to replace 
patterns already given in the X-ray Powder Data File issued 
by the American Society for Testing and Materials, and 
twenty-six patterns indicated by asterisks are for substances 
not previously included. The X-ray Powder Data File is a 
compilation of diffraction patterns from many sources and is 
used for the identification of unknown crystalline materials 
by matching spacing and intensity measurements. The 
patterns were made with a Geiger counter X-ray diffractom- 
eter, using samples of high purity. When possible, the 
d-values were assigned Miller indices determined by compari- 
son with calculated interplanar spacings and from space 
group extinctions. The densities and lattice constants 
were calculated, and the refractive indices were measured 
whenever possible. 


Mechanical behavior of crystalline solids (Proceedings of an 
American Ceramic Society Symposium, New York City, 
April 1962), NBS Mono. 59 (Mar. 25, 1963), $1.75. 

This Monograph represents the Proceedings of a Symposium 
on The Mechanical Behavior of Crystalline Solids, held under 
the auspices of the Ceramic Educational Council of the 
American Ceramic Society, with the cooperation of the 
National Bureau of Standards, and under the sponsorship 
of the Edward Orton Junior Ceramic Foundation, and the 
Office of Naval Research. The Symposium took place at 
the 64th Annual Meeting of the American Ceramic Society, 
in New York, on April 28 and 29, 1962. 


Testing of metal volumetric standards, J. C. Hughes and B. C. 
Keysar, NBS Mono. 62 (Apr. 1, 1963), 15 cents. 

The National Bureau of Standards has for many years ¢ali- 
brated and certified metal measures which are used as stand- 
ards by weights and measures officials and others in the cali- 
bration of instruments for measuring the volumes of fluids. 
No complete specifications or tolerances for these standards 
have ever been published, however, nor have standardized 
procedures for the calibration and use of the liquid measures 
been available. 

The information contained in this Monograph should assist 
in the purchase of quality instruments and the proper use 
of the standards in calibrating other measures for liquids and 
gases. 


Reduction of data for piston gage pressure measurements, 
J. L. Cross, NBS Mono. 65 (June 17, 1963), 15 cents. 
Pressure measurements made with piston gages are affected 
by gravity, temperature, pressure, and several other variables. 
For accurate determinations of pressure the calculations must 
take these variables into account. A general equation is 
developed and simplified procedures for calculating pressure 
are illustrated. 


Tabulation of data on receiving tubes, C. P. Marsden 
J. K. Moffitt, NBS Handb. 83 (May 28, 1968), $1.25. 
sedes Handb. 68.) 


and 
(Super- 


A tabulation of Receiving-Type Electron Tubes with some 
characteristics of each type has been prepared in the form of 
two major listings, a Numerical Listing in which the. tubes 
are arranged by type number, and a Characteristic Listing 
in which the tubes are arranged by tube type and, further 
ordered on the basis of one or two important parameters. 
The tabulation is accompanied by a listing of similar tube 
types and basing connections for the listed tubes. 


Transistorized building blocks for data instrumentation, 
R. L. Hill, NBS Tech, Note 168 (Apr. 1, 1963), 55 cents. 

The National Bureau of Standards has developed a rumber 
of modular transistorized digital circuits that have been 
used in automatizing many data recording and preliminary 
processing tasks encountered in its scientific operations. 
These versitile building blocks can be connected together 
systematically to form digital circuits that accept raw data 
from experimental equipment and transpose these data into 
a form suitable for input to a high-speed electronic computer 
Kach assembly of packages can be tailored to fit the special 
requirements of the project and can be used at the site of the 
experiment. The output from the system can be: 1) fed 
directly to a computer, 2) recorded on a medium (paper 
tape, magnetic tape, etc.) suitable for computer input at a 
later date, or 3) used to drive display equipment that keeps 
the scientist informed of the progress of his experiment. 

As a result of experience in the application of these units, 
some of the original packages have been modified and addi- 
tional types developed. In addition to describing the modi- 
fied and new package types, this report also includes a de- 
scription of a new series of packages consisting of identical 
circuitry, but utilizing a different type of mating connector 
and a smaller circuit-board. 


Phototypesetting of computer output, an example using 
tabular data, W. R. Bozman, NBS Tech, Note 170 (June 25, 
1963), 10 cents. 

A photocomposition machine controlled by the magnetic 
tape output from a computer was used to prepare a 559-page 
table of atomic transition probabilities at the National Bureau 
of Standards. This method makes possible the publication 
of computed data in high quality typography in a reasonable 
time and at a reasonable cost. Many styles of type are 
readily available to the programmer including Greek, italic, 
mathematical symbols, upper and lower case alphabets, etc. 


Practical methods for calibration of potentiometers, D. 
Ramaley, NBS Tech. Note 172 (Mar. 25, 1963), 30 cents. 
Potentiometer circuitry, particularly as related to calibration, 
is discussed with the primary consideration given to the 
required circuit measurements. The more feasible means 
of calibrating potentiometers are described in considerable 
detail. Emphasis is placed upon the use of the Universal 
Ratio Set as the basic implement for accomplishing the major 
portion of potentiometer calibrations. 


Table of attenuation error as a function of vane-angle error 
for rotary vane attenuators, W. Larson, NBS Tech. Note 177 
(May 20, 1963), 75 cents. 

The table of attenuation error as a function of vane-angle 
error gives the error in decibels caused by vane misalinement 
which is common in the rotary vane attenuator. The attenu- 
ation errors corresponding to vane-angle errors ranging from 
zero to 0.499° (in inerements of 0.001°) are presented for 
selected angles over the range of attenuation values from 
0.01 to 70 db. The table is divided into the following inter- 
vals of attenuation value increments: 0.01—-0.1 db in 0.01-db 
increments, 0.1—1.0 db in 0.1—db increments, 1-20 db in 1—db 
increments, and 20-70 db in 5—-db increments. 

With the aid of this table, the calibration data of a rotary vane 
attenuator can be analyzed for numerous characteristics, in- 
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cluding the following: misalinement between rotor and stator 
sections, realinement techniques, resettability, and backlash. 


Tabulation of published data on Soviet electron devices, C. P. 
Marsden, NBS Tech. Note 186 (June 3, 1968), 45 cents. 

This tabulation includes published data on Soviet electron 
devices as collected from various publications, mostly hand- 
books published by the various ministries and institutes of 
the USSR. Information is given on all active devices ranging 
from receiving to microwave devices, semiconductor devices, 
and various miscellaneous devices such as, for example, photo- 
graphic flash tubes and thermistors 


Calibration of volt-ampere converters, L. S. Williams, 
Tech. Note 188 (Apr. 25, 1963), 20 cents. 

These notes have been prepared to describe the National 
Bureau of Standards calibration services for volt-ampere con- 
verters (or transfer volt-ammeters), to suggest procedures for 
d-e standardization in the user’s laboratory, and to describe 
a voltage comparator which can be used to make such eali- 
brations quickly and easily. 


NBS 


Tables describing small-sample properties of the mean, 
median, standard deviation, and other statistics in sampling 
from verious distributions, C. Eisenhart, L. 8. Deming, and 
C.S. Martin, NBS Tech. Note 191 (June 14, 1963), 20 cents. 
This note includes a collection of tables useful for study of the 
sampling distributions of some frequently-used statistics, with 
brief discussions of their construction and use. (1) The prob- 
ability level P(e,n) of any continuous parent distribution 
corresponding to level « of the distribution of the median. 
2) Probability points of certain sample statistics for samples 
from six distributions: normal and double-exponential (mean, 
median), rectangular (mean, median, midrange), Cauchy, 
Sech, Sech? (median). In all the above tables, the sample size 
n=3(2)15(10)95 and the probability levels are e=.001, .005, 
Ol, .025, .05, .10, .20, .25. Together with the tables listed 
under (2) are given the values of certain ratios useful for 
comparing the various statistics. (3) Probability that the 
standard deviation of a normal distribution will be under- 
estimated by the sample standard deviation s and by unbiased 
estimators of ¢ based on s, on the mean deviation, and on the 
sample range. Divisors are given for obtaining the corre- 
sponding “‘median unbiased”’ estimators. 


National standard reference data program, background in- 
formation, NBS Tech. Note 194 (June 1963), 25 cents. 

Plans are proposed for a National Standard Reference Data 
System that will provide critically evaluated data in the 
physical sciences on a national basis. It will be conducted 
as a decentralized operation across the country, with central 
coordination and administration by NBS. Data will be 
centrally stored at NBS and disseminated through a series 
of services tailored to user needs in science and industry. 


New absolute null method for the measurement of magnetic 
susceptibilities in weak low-frequency fields, C. T. Zahn, 
Rev. Sci. Instr. 343 No. 3, 285-291 (Mar. 1963). 

Use is made of the magnetic equivalence of a uniformly 
polarized volume of paramagnetic material and a solenoid 
carrying electric current, to design a permanent variable 
standard of magnetic susceptibility. Such a standard is 
incorporated into a magnetic susceptibility bridge in a simple 
manner, surrounding the specimen; and the bridge is thereby 
transformed into an absolute null instrument of high ac- 
curacy and sensitivity, and of great ease and low cost of 
construction and operation. By this method numerous 
particular advantages of other methods are combined; and 
some of their notable limitations are overcome. A_ pre- 
liminary application was made showing that the bridge per- 
forms as expected. Important features in the 


; ; design of 
this bridge are discussed. A _ detailed 


consideration of 


sources of error suggests that it may eventually be possible 
by this method to obtain greater absolute : 


accuracy than by 
other known methods. 


A method for measuring the instability of resistance strain 
gages at elevated temperatures, R. L. Bloss and J. T. Trumbo, 
ISA Trans. 2; No. 2, 112-116 (Apr. 1963). 


a 
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The usefulness of resistance strain gages at elevated tempera- 
tures is frequently limited by the instability of gage resistance 
with time. Methods and equipment that 


been de- 
veloped to measure this effect are described. 


have 


of a 
Tran S 


Maximum efficiency 
R. W. Beatty, JEEE 
MTT-11, 94 (Jan. 1963). 

Given the scattering coefficients of a 2-arm waveguide 
junction, an equation is presented to calculate I'y, the 
reflection coefficient of the load for which the efficiency of 
a 2-arm waveguide junction is ny, the maximum efficiency. 
Once Ty has been calculated, an equation is given to deter- 
mine ny and the relationship between ny and A;, the in- 
trinsic attenuation of the waveguide junction, is given. 


two-arm waveguide junction, 
Vicrowave Theory and Tech. 


Applications of a semiconductor-surface-state charge-storage 
device, L. J. Swartzendruber, Solid-State Elec. 6, 59-61 
(Pergamon Press, Inc., New York, N.Y. 1968). 

Several possible applications of a new two terminal semicon- 
ductor device which utilizes surface phenomena to produce a 
charge storage effect are described. The major advantage of 
the device lies in the magnitude of the charges which can be 
stored and the ease with which it can be controlled by small 
bias currents. 


Audio-frequency compliances of prestressed quartz, fused 
silica, and aluminum, M. Greenspan and C. Tschiegg, Proc. 
Fourth Intern. Congr. on Acoustics, Part I, Paper P12 (Copen- 
hagen, Denmark, Auq. 21-28, 1962). 

An attempt was made to find the compliances 
associated with dispersions found by Fitzgerald. Com- 
pliances were obtained from resonant frequencies of fixed- 
free composite reeds. Prestress was either piezoelectrically 
or thermally induced. No excess compliances were observed. 


excess 


Ten-kilocycle pound-type kylstron stabilizer, H. Ek. Radford, 
Rev. Sci. Instr. 34, No. 3, 8304-805 (Mar. 1963). 

Through a simple circuit modification, commercial klystron 
frequency stabilizers of the FM type can be made to function 
alternatively as CW Pound-type stabilizers, with greater 
spectral purity of the klystron output. The performance of 
such a stabilizer is discussed. 


Kihara parameters and second viral coefficients for cryogenic 
fluids and their mixtures, J. M. Prausnitz and A. L. Myers, 
A. I. Ch. E. Journal 9, No. 1, 5-11 (Jan. 1963). 

The volumetric properties of sixteen fluids of interest in eryo- 
genic engineering have been used to calculate second viral 
coefficients over as large a temperature range as possible. 
These coefficients were then fitted to theoretical expressions 
based on the Kihara potential function. For helium, 
hydrogen, and neon quantum corrections were applied. For 
nitrogen, carbon dioxide, and acetylene corrections for 
quadrupole interactions were made. It was found that the 
theoretical expressions give an extremely good fit of all 
reliable experimental data. The theoretical expressions may 
therefore be used with confidence to predict volumetric 
behavior at very low temperatures where data are frequently 
unavailable. 

With the aid of semiempirical mixing rules the theoretical 
expressions may be used to predict second viral coefficients 
for mixtures. Agreement with the very limited amount of 
experimental mixture data is satisfactory. Finally it is 
shown that calculations based on the Kihara potential may 
be employed to make useful predictions of phase equilibria 
such as the solubility of a solid in a compressed gas. 


Transparent rigid mount for vacuum stopcock, M. M. Ander- 
son, Rev. Sci. Instr. 34, No. 2, 178 (Feb. 1963). 

A transparent block of plastic, moulded around a glass stop- 
cock, then drilled and tapped for mounting, reduces vacuum 
system breakage and allows visual inspection of the grease 
seal. 


Design of low voltage electron guns, J. A. Simpson and C. E. 
Kuyatt, Rev. Sci. Instr. 34, No. 3, 265-268 (Mar. 1963). 

It is shown that by use of a multistage technique in which 
electrons are drawn from a cathode by a high potential and 





decelerated to the required final energy, guns can be designed 
capable of forming beams in which the current is limited only 
by space charge in the beam itself. The design principles 
and procedures are given and illustrated by two examples of 
electron guns giving highly collimated beams and operating 
at energies of 30 and 500 eV. The measured currents 
obtained are somewhat greater than the space charge limited 
beam maximum because of ion neutralization. 


Thermometry, 'ow temperature, R. P. Hudson, Encyclopaedic 
Dictionary of Physics 7, 323-325 (1962). 

A discussion is given of apparatus and methods for thermom- 
etry in the range 1°-100°K in a style and brevity suited to 
an entry in a scientific encyclopaedia. 


Effect of outdoor exposure on some properties of chrome- 
retanned leather, T. J. Carter, J. Am. Leather Chemists 
Assoc. LVIII, No. 3, 155-160 (Mar. 1963). 

Two groups of specimens of chrome-retanned leather were 
subjected to outdoor exposure and the changes in physical and 
chemical properties were determined. One group was sub- 
jected to all weather conditions while the other group was 
shielded from the sun. The effect of the exposure was 
determined by measurements of physical properties, such as 
flex tension and impact resistance, elongation, stitch tearing 
strength, tensile strength, change in area, relative stiffness, 
and shrinkage temperature, and of chemical properties, such 
as grease content and pH. 

Results show that specimens shielded from the sun changed 
little in impact resistance, stitch tear strength, elongation of 
tensile strength specimens, and stiffness, but changed sig- 
nificantly in flex tension resistance, elongation (due to 
flexing), area, and shrinkage temperature. Specimens ex- 
posed to the sun showed deleterious changes in all the physical 
properties studied. Shrinkage temperature declined  sig- 
nificantly under both conditions, the decrease being somewhat 
greater for specimens exposed to direct sunlight. The 
chemical properties, grease content and pH values, showed 


only slight changes under either condition of exposure. 


APPA-TAPPI reference material program. II. Effectiveness 
of a reference material in reducing the between-laboratory 
variability of TAPPI standard T 414 m-49 for — 
tearing resistance of paper, T. W. Lashof, Tappi 46, No. 3, 
145-150 (Mar. 1968). 
The effectiveness of a reference material was predicted on 
the basis of the results of the first round robin which was 
reported in Part I. The analysis was modified for the 
second round robin so as to provide corrections in terms of 
measurements on the reference material. If the current 
TAPPI procedure is followed (5 replications, no standard 
reference material), the total coefficient of variation, in- 
cluding both within- and between-laboratory variability, 
is about 4% to 5%, as shown in both round robins. As 
shown in this second part, a standard reference material 
and increased replication may be used to reduce the total 
coefficient of variation to about 3% It is also shown that 
the correction curve or nomograph based on the measure- 
ments on the reference material may be used for at least 
four to five months, provided that there are no changes in 
observer, instrument, or conditions. 

The limitation to further reduction in the total coefficient 
of variation is V(A), the random interaction between inter- 
fering properties of the materials being tested and laboratory 
conditions or instrument peculiarities. It is shown that a 
portion of this is probably due to insufficient control of rela- 
tive humidity. Sinee V(A) is higher for laboratories using 
new-type instruments than for those using old-type instru- 
ments, further work must be done to determine whether this 
is due to the instruments or to laboratory conditions. 


Mullen 
Vacuum Symp., 


Some 4 eg? of a simple cryopump, L. O. 
and R. B. Jacobs, 1962 Trans. Ninth Natl. 
Am. Vacuum Soc., pp. 220-226 (1962). 

A simple and easily definable cryopump was constructed 
as a stage in a pumping system, and data were obtained to 
permit the computation of pumping speeds, performance 
decay and capture coefficients. Information on the pumping 
of CO, and Ny, as well as outgassing vapors, by surfaces at 
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77°K and 20°K is presented. 
tests is 5(10)~! torr to 5(10) and the gas flow 
3.6 10 torr liters cm.~? to 3 6(10) : 
cm.~? see.~!, Pumping speeds higher than those 
reported and higher than theory predicts were 
the results are discussed in detail. 


The 
> torr 
sec.~! 


of the 
range 18 
torr liters 
previously 
obtained, 


pressure range 


Experimental investigation of Fabry-Perot interferometer, 
R. W. Zimmerer, Proc. [EEE 51, 475-476 (Mar. 1963). 
Preliminary measurements of the microwave performance of 
Fabry-Perot interferometers with spherical mirrors is pre- 
sented and compared with theory. Of particular interest is 
the evidence of the stop band recently predicted by Boyd and 
Kogelnik. 


The measurement of moisture boundary layers and leaf tran- 
spiration with a microwave refractometer, D. M. Gates, M. J. 
Vetter, and M. C. Thompson, Jr., Nature 197, 1070-1072 
(Mar. 16, 1963). 

A microwave refractometer has been used as a hygrometer to 
measure the moisture gradient found near a free water surface 
and near the surface of aleaf. Interesting transpiration effects 
were observed for begonia and bean leaves when the leaves 
were stimulated with light. The instrument samples the air 
through a small orifice and thereby produces very little dis- 
turbance to the moisture boundary layer under investigation. 


New scale of nuclidic masses and atomic weights, 
Nature 194, No. 4829, 621-624 (May 10, 1962). 
This is an article written at the request of the 
“NATURE.” 

It reviews the considerations that led to the adoption by the 
International Unions of Physics and Chemistry of a new scale 
of nuclidic masses and atomic weights based on C!?=12. 


m W ichers, 


Editor of 


Calibration of photogrammetric lenses and cameras at the 
National Bureau of Standards, F. E. Washer, 
Eng. XXIX, No. 1, 113-119 (Jan. 1963). 

A summary of calibrations performed at the National Bureau 
of Standards on lenses and cameras that are used in precise 
photogrammetric work is given. Brief description of the 
photographic and visual calibrations most frequently re- 
quired are given. This paper includes a list of publications by 
members of the NBS staff that pertain to problems of lens 
and camera calibration. 


Photogrammetric 


Building a transistor tester, G. 
Electronics 36, 16, 56 (Apr. 19, 1963). 
A simple Baseked is de alt ha P for 
transistor parameters: leakage current 
current amplification. 


Montgomery, 


measuring two de 
and common-emitter 


A magnetic amplifier for use with diode logic, E. 
Proc. IEEE 1963 er Conf. Nonlinear 
T-149, 8.6-1 to 8.6-6 (Apr. 1963). 

A digital amplifier of pai noncritical design incorporating 
an emitter-follower and a small magnetic amplifier is de- 
scribed. Timing and some of the operating power are provided 
by a 300-ke 2-phase 7-volt sine-wave source. In structure 
and mode of operation, the amplifier is particularly suited 
for use with two-level diode gating to provide the AND and 
OR logical operations. A NOT-amplifier provides negation 
with amplification. The volt-second transfer characteristic 
of the stage critically determines the stability of propag zation 
of binary signals. Factors governing the required shi ape of 
this transfer characteristic are discussed. 


W. Hogue, 
Magnetics No. 


The speed . light, G. MeNish, 
Nos. 3 and 4, 138- 143 (Dec. 1962). 
Numerous me easure ments of the speed of light published ‘during 
the last 30 vears lead to widely divergent results as compared 
with the assigned experimental uncertainties. Because of 
wide diversity in the methods employed in the measure- 
ments, all of the data may not be combined effectively in a 
grand average. Sufficient data had been obtained by the 
geodimeter method to group them and derive a statistical 
estimate of the uncertainty in the speed of light by this 
method. This result, and conclusions reached from careful 
examination of several experiments, leads to the conclusion 


IRE. Trans. Instr: 


I-11, 





that the value 299,792.5 km per sec which has been inter- 


nationally adopted for use in radio propagation and geodetic 
is very close to the best value and not likely to be in 
error by as much as one part in one million. 


work 


Performance characteristics of split-type residential 
to-air heat pumps, J. C. Davis and P. R. Achenbach, S 
B Inst. Intern. Refrige 1-7 (1961-1962). 
This paper presents test results obtained during a laboratory 
study of six split-tvpe residential air-to-air heat pumps, a 
tvpe more widely used than others for residential application 
in the United States in the last few vears. For residential 
use, these currently selected to satisfy the cooling 

quirements of the house under design summer conditions, 
and, if necessary, the heating capacity of the compression- 
evcle is supplemented with electric resistance heaters. The 
investigation showed that the compression-cycle heating 
capacity of the heat pumps equipped with expansion valves 
nereased linearly with increasing outdoor temperature at 
constant indoor temperature and humidity, whereas the 
capacity of tne heat pumps equipped with capillary tubes also 
increased but tended toward constant values at higher 
temperatures. It was shown that the cooling capacity of the 
heat pumps decreased linearly with outdoor temperature, 
when indoor temperature and relative humidity were held 
\t constant outdoor temperatures, cooling capac- 
either with increasing indoor tem- 
relative humidity. The changes in the 
latent and sensible fractions of the total cooling capacity 
1 by change in indoor r¢ humidity, and the effect 
temperature on coefficient of performance under 
and cooling conditions, are also reported. 


air- 
ippl. 


ration, p 


Svstems are 


constant 
itv increased linearly 
perature or indoor 
cause lative 
of outdoor 


leating 


Millimeter wavelength resonant structures, R. W. Zimmerer, 
M V Anderson, G L. Strine, . \ TY; IEEE Trans. 
Wi oware Theo / Tech. MMT 11, L2 4 + (Mar. 1963). 
This paper discusses the nadia ‘a millimeter wave 
Fabry-Perot using both planar and spherical re- 
flectors. It discusses the equivalent circuits of planar 
reflectors and the method of obtaining efficient power transfer 
into the resonators. 


resonators 
also 


A simple environmental chamber for rotating beam fatigue 


testing machines, J. A. 

(80-482 (June 1963 

\ gas-tight sleeve, made from transparent plastic, permits 

fatigue testing under controlled atmosphere conditions. Tests 

of magnesium and aluminum alloys have shown that changes 
change the fatigue strength by more than 


Bennett, Mater. Res. Std. 3, No. 6, 


in humidity may 


10% 


Evidence regarding the mechanism of fatigue from studies of 
environmental effects, J. A. Bennett, Acta Met. 11, No. 7, 


799-800 (Jul 1963 


Fatigue tests of aluminum alloy speci 
under controlled humidity conditions tesults of tests in 
which the environment is changed during the test show that 
there is an initial period during which the humidity has no 
effect. This is interpreted to indicate that the deformation 
during that period is not localized 


nens are being conducted 


Calorimetric calibration of the electrical energy measurement 
in an exploding wire experiment, D. Tsai and J. H. Park, 
Erploding Wires 2, 27-107 (Plenum Press, Inc., New York, 
N.Y., 1962). 

A discussion is presented on the requirements and the methods 
for measuring the current and voltage during the transient 
discharge of a capacitor bank employed in an exploding wire 
experiment \ method is described for accurately calibrating 
the measured current, voltage, and electrical energy by com- 
paring the calorimetric heating of a fixed resistance element 
with the electrical energy dissipated in the element. Results 
show that the accuracy of the energy measurement is about 


1-27, 


Oil baths for saturated standard cells, P. H. 
ISA J. 9, No. 12, 47-50 (De c. 1962). 

The increasing use of saturated standard cells in 
has growing need for information on 


Lowrie, Jr., 


industry 


caused a equipment 


associated with tneir use. This paper discusses oil baths 
suitable for the close temperature control of these cells and 
describes the oil baths in use at NBS Boulder Laboratories. 
In these baths, the temperature is controlled by a modified on- 
off system that limits cyclic variations to less than +0.001°C 
from the mean temperature. The mean does not change more 
than 0.002°C per day in an environment in which the ambient 
temperature may change by as much as 2°C during the day. 


Performance of the ——- fluoride film hygrometer ——- 
on radiosonde flights, F. E.. Jones, J Res. 68, N 

9, 27385-2751 (May 1, 1963) 

Ten flights of the barium fluoride film electric hygrometer 
element in a modified radiosonde on the same train with a 
conventional ithium chloride element in an AN/AMT-11 
radiosonde were made from the grounds of the National 
Bureau of Standards (Wash., D.C.) during the period January 
16 through July 21 of 1961. The flights were intended to 
provide information on the performance of the barium fluoride 
element under conditions encountered in routine radiosonde 
flights and to provide information to be used in assessing the 
value of the element as a research tool 

The results for the flights verified results of laboratory 
in several areas. The element responded to changes in 
humidity over a range of indicated relative humidity, RH, 
of 1.5 to 100° in the temperature range 33.1 to & ae ah 
preflight room temperature calibrations indicated that the 
ten elements flown were typical, in this respect, of elements 
tested under laboratory conditions, exposure to high humidity 
and passage through precipitation had no apparent effect on 
the functioning of the element, the rapid response of the 
element and its ability to resolve fine humidity structure 
were Gemonstrated, indications of saturation or near-satura- 
tion were correlated with U.S. Weather Bureau surface 
observations of clouds and radar weather observations. 

In several of the flights, a strong correlation exists between 
changes in indicated RH and ambient temperature lapse 
rates. A sharp drop in indicated RH within the first several 
hundred feet above the surface in at least five of the flights 
is possibly related to boundary layer phenomena at or near 
the surface. In two of the flights the element detected 
supersaturation with respect to ice 

Although the instability with time of the barium fluoride 
element, in its present state of development, precludes its 
use in routine radiosonde flights, the ten flights indicated 
the value of the element for experimental use and as a research 
tool. 

In addition to the 
flight was made of 


Geophys. 


tests 


flights of the barium fluoride element, one 
a lead iodide film hydrometer element. 


Intererfence fringes with me path difference using He-Ne 
laser, T. Morokuma, K. Nefflen, T. R. Lawrence, and T. M. 
Klucher, J. Opt. Soc. Am. 53. No. 3, 394-395 (Mar. 1963). 
Interference fringes have been obtained in a Michelson inter- 
ferometer with path length differences up to 9 meters using a 
helium neo laser as lignt source. 


Hyrdogen retention system for pressure calibration of micro- 
phones in small couplers, W. Koidan, J/. 
35, No. 4, 614 (Apr. 1963). 

The pressure calibration of microphones using small hydrogen- 
filled couplers can be facilitated by connecting relatively 
large containers of hydrogen to the capillary tubes of the 
coupler and using additional capillary tubes to vent the 
large containers to the atmosphere. 


Acoust, Soc. Am. 


Other NBS Publication 


Journal of Research 67A (Phys. and Chem.), 

Aug. 1963), 70 cents. 

Symmetry splitting of equivalent sites in oxide 
related mechanical effects. J. B. 
Peiser, and E. P. Levine. 

Relaxation modes for trapped crystal point defects. A. D. 
Franklin. 

A note on the galvanomagnetic and thermoelectric coeffi- 
cients of tetragonal crystalline materials. W. ) 
Hernandez, Jr., and A. H. Kahn. 


No. 4 (July— 


crystals and 
Wachtman, Jr., H. 8. 
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Photolytic behavior of silver 

Correlation of muscovite sheet 
apparent optic angle, and absorption spectrum. St. 
Ruthberg, M. W. Barnes, and R. H. Noyce. 

Thermodynamic properties of magnesium oxide and beryllium 
oxide from 298 to 1,200 °K. A. C. Victor and T. B. 
Douglas. 

Heat exchange in adiabatic calorimeters. E. D. West. 

Preparation of anhydrous single crystals of rare-earth halides. 
N. H. Kiess. 

A phase study of the system: 
its significance in oxalic 
burger and J. L. Torgesen. 

Wavelength calibrations in the far infrared (30 to 1000 
microns). K. N. Rao, R. V. de Vore, and E. K. Plyler. 

On the fourth order Hamiltonian of an asymmetric rotor 
molecule of orthorhombic symmetry. Wm. B. Olson and 
H.C. Allen, Jr 

Measurement of the thickness and refractive index of very 
thin films and the optical properties of surfaces by ellipsom- 
etry Fr MeCrackin, Passaglia, R. R. Stromberg, 
and H. L. Steinberg. 

Color phenomena associated with energy 
glows and atomic flames. A. M. 


iodide. G. Burley. 


mica on the basis of color, 


oxalic acid/acetic acid/water; 
acid crystal growth. J. Strass- 


transfer in after- 
Bass and H. P. Broida. 


Journal of Research 67A (Phys. and chem.), 

Oct. 1963), 70 cents. 

Reduction of space groups to subgroups by 
strain. H. S. Peiser, J. B. 
Dickson. 

High-temperature 


homogeneous 


Wachtman, Jr., and R. W. 


thermodynamic functions for zirconium 
and unsaturated zirconium hydrides. T. B. Douglas. 

Heat of oxidation of aqueous sulfur dioxide with gaseous 
chlorine. W. H. Johnson and J. R. Ambrose. 

Thickness of adsorbed polystyrene layers by ellipsometry. 
R. R. Stromberg, E. Passaglia, and D. J. Tutas. 

Melting temperature and change of lamellar thickness with 
time for bulk polyethylene. J. J. Weeks. 
Precise coulometric titrations of potassium 

G. Marinenko and J. K. Taylor. 

Resolution limits of analyzers and 
Edith L. R. Corliss. 

Synthesis, purification, and physical 
twelve-carbon hydrocarbons. T. W. 
ki. L. Compere, Jr., and F. L. Howard. 

Reactions of polyfluorobe nzenes with nucleophilic reagents. 
L. A. Wall, W. J. Pummer, J. E. Fearn, and J. M. Anto- 


nucci. 


dichromate. 


oscillatory systems, 
properting of seven 
Mears, C. L. Stanley, 


Journal of Research 67B (Math. and Math. PHys.), No. 

3 (July—Sept. 1963), 75 cents. 

Remarks on hypo-elasticity. C. Truesdell. 

Error bounds in the pointwise approximation of solutions 
of elastic plate problems. J. H. Bramble and L. E. Payne. 

Effect of error in measurement of elastic constants on the 
solutions of problems in classical elasticity. H. Bramble 
and L. EF. Payne 

Kigenfunctions of the f* configuration. J. C. 

Zeros of first derivatives of Bessel ee of the 
Ji(z), 21<n<5l, O<2r< 100. . Ww. 
H. Reismann. 


Kisenstein. 
first kind, 
Morgenthaler and 


Journal of Research 67D (Radio Prop.), 

1963), 70 cents. 

Ionospheric VHF scattering near the magnetic equator during 
the International Geophysical Year. R. Cohen and K. L. 
Bowles. 

Radio pulse propagation by 
ionosphere. J. R. Johler. 

Field of a horizontal magnetic dipole 
magnetoplasma halfspace. G. 
H. M. Swarm. 

Reflection of VLF radio waves from an inhomogeneous iono- 
sphere. aie II. Perturbed exponential model. J. R. 
Wait and L. C. Walters. 

Collisional det: hme nt and the formation of an ionospheric 
C region. E. T. Pierce. 

Magnetic aN and Coriolis effects on a 
suspended rotating sphere. J. C. Keith. 


No. 5 (Sept.—Oct. 


a reflection process at the lower 


in the 
Tyras, A. 


presence of a 
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Radiation field characteristics of lightning discharges in the 
band 1 ke/s to 100 ke/s. W. L. Taylor. 

Low-frequency radio propagation into a moderately rough 
sea. LD. F. Winter. 

VLF superdirective loop arrays. E. W. Seeley. 

Curves of ground proximity loss for dipole antennas (a digest). 
L. E. Vogler and J. L. Noble. 

Observations and results from the “ 
ment - ey Oe ly observe the 
Watts, . . Koch, and R. M. 


hiss recorder’ an instru- 
VLF emissions. J. M. 
Gallet. 


Influence of a sector ground screen on the 
antenna, J. R. Wait and L. C. 
(Apr. 15, 1963), 25 cents. 

Refractive indices and densities of aqueous solutions of invert 
sugar, C. F. Snyder and A. T. Hattenburg, NBS Mono. 
64 (June 7, 1963), 15 cents. 

Radiobiological dosimetry. Recommendations of the 
national Commission on radiological units and measure- 
ments, NBS Handb. 88 (Apr. 30, 1963), Supersedes parts 
of Handb. 78. Handbooks 84 through 89 will completely 
replace H78, 25 cents. 

Quarterly radio noise data, 
Crichlow, R. T. Disney, and M. A. 
Note 18-15 (Mar. 1, 1963), 45 cents. 

Quarterly radio noise data, September, 
1962, W. Q. Crichlow, R. T. 
NBS Tech. Note 18-16 (June 10, 1963), 60 cents. 

Mean electron density variations of the quiet ionosphere, 
November 1959, J. W. Wright, L. R. Wescott, and D. J. 
Brown, NBS Tech. Note 40-9 (Apr. 22, 1963), 35 cents. 

Mean electron density variations of the quiet ionosphere, 
December 1959, J. W. Wright, L. R. Wescott, and D. J. 
Brown, NBS Tech. Note 40-10 (Mar. 24, 
cents. 

The error rates in multiple 
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Walters, NBS Mono. 60 
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Jenkins, NBS Tech. 
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Disney and M. A. Jenkins, 


1963), 35 

FSK systems and the signal-to- 
noise characteristics of FM and PCM-FS systems, H. 
Akima, NBS Tech. Note 167 (Mar. 25, 1963), 40 cents. 

Bibliography on atmospheric aspects of radio astronomy, 
inciuding selected references to related fields, W. Nupen, 
NBS Tech. Note 171 (May 1, 1963), $2.00. 

Tables to facilitate the determination of the ferrimagnetic 
resonance linewidth of non-metallic magnetic materials, C. C. 
Preston and W. E. Case, NBS Tech. Note 173 (Apr. 15, 
1963), 25 cents. 

Curves of ground proximity loss for dipole 
Vogler and J. L. Noble, NBS Tech. 
1963), 30 cents. 

An interpolation procedure for calculating atmospheric band 
absorptions from laboratory data, L. Droppleman, L. R. 
Megill, and R. F. Calfee, NBS Tech. Note 178 (June 3, 
1963), 20 cents. 

Relative power transmission characteristics of the ear and 
skull from hearing threshold data, E. L. Smith, Proc. 
Fourth Intern. Congr. on Acoustics, Part I, Paper H48 
(Copenhagen, Denmark, Aug. 21-28, 1962). 

Polymer ~8r at the U.S. National Bureau of Standards, 
Part I, Kline, SPE J. 19, No. 3, 278-283 (Mar. 1963). 

Polymer rese Past at the U.S. National Bureau of Standards, 
Part 2, G. M. Kline, SPE J. 19, No. 4, 403-408 (Apr. 1963). 

Radiation and the world we live in, L. 8S. Taylor, Radiology 
80, No. 3, 358-368 (Mar. 1963). 

Natural and synthetic rubbers, E. J. Parks and F. J 
Anal. Chem. 35, No. 5, 160R-178R (Apr. 1963). 

Degradation of polymers, L. A. Wall and J. H. Flynn, Rubber 
Chem. Technol. XXXV, No. 5, 1157-1221 (Dee. 1962) 

Chemical and magnetic enhancement of perturbed lines in the 
violet spectrum of CN, H. E. Radford and H. P. Broida, 
J. Chem. Phys. 38,, No. 3, 644-657 (Feb. 1, 1963). 

Profiles of Stark-broadened Blamer lines in a hydrogen plasma, 
W. L. Wiese, “i R. Paquette, ae J. E. Solarski, Phys. 
Rev. 129, No. 1225-1232 (Feb. 1, 1963). 

Low-temperature ‘thermonntry, K. D. Timmerhaus, Book, 
Applied Cryogenic Engineering, ed. R. W. Vance, ch. 4> 
60-103 (John Wiley & Sons, Inc., New York, N.Y., 1962) 

Magnetic susceptibilities and dilution effects in low-spin dt 
complexes: Osmium (IV), R. B. Johannesen and G. A. 


antennas, L. E. 
Note 175 (May 20, 
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Candela, Inorg. Chem. 2, 67—72 (1963). 

Standards and the microwave profession, J. M. Richardson, 
IRE Trans. Microwave Theory and Tech. MTT-10, No. 6, 
413-415 (Nov. 1962). 

Delay time of polar-cap blackout and its relation to decay 
time of geomagnetic disturbance, C. 8S. Warwick, J. 
Geophys. Res. 68, No. 5, 1561-1562 (Mar. 1, 1963). 

Electrolytic conductance of ammonium dihydrogen phosphate 
solutions in the saturation region, J. L. Torgesen and A. T. 
Horton, J. Phys. Chem. 67, 376-381 (1963). 

Determination of source self-absorption in the standardization 
of electron.capturing radionuclides, 8. B. Garfinkel and 
J. M. R. Hutchinson, Intern. J. Appl. Radiation and 
isotopes 13, 629-639 (1962). 

Standards for the 70s, W. A. Wildhack, Ind. Res. §, No. 3, 
15-20 (Mar. 1963). 
Nuclear orientation, E. Ambler, Book, Methods of Experi- 
mental Physies 5, sect. 2.4.2.3, 196-214 (Academic Press 

Inc., New York, N.Y., 1963). 

The ionosphere over Antarctica, W. R. Piggott and A. H. 
Shapley, Antarctic Research, Geophysical Mono. 7, p. 111 
126 (1962 

Mean first-passage 
molecules, K. E. Shuler and G. H 
38, No. 2, 505-509 (Jan. 15, 1963). 

The personal side of a research project, A. T. McPherson, J 
Wash. Acad. Sci. 53, No. 3, 63-66 (Mar. 1963). 

Wavelengths, energy levels, and pressure shifts in mercury 
198, VY. Kaufman, J. Opt. Soc. Am. 52, No. 8, 866-870 
(Aug. 1962). 

On the isomerization of isobutyl radicals, J. 


dissociation of diatomic 
Weiss, J. Chem. Phys. 


times and the 


R. MeNesby and 


W. M. Jackson, J. Chem. Phys 38, No. 3, 692-693 (Feb. 1, 
1963). 

The electrophoretic mobility of asphaltenes in nitromethane, 
J, es 


99: 
lg 


Wright and R. R. 

Mar. 1963). 

A suggestion for improving forecasts of geomagnetic storms, 
Y. Hakura and J. V. Lincoln, J. Geophys. Res. 68, No. 5, 
1563-1564 (Mar. 1, 1963). 

Observation of a 6300 A are in France, America, and Austra- 
lia, F. E. Roach, D. Barbier, and R. A. Dunean, Ann. 
Geophys. 18, 390-391 (Oct.—Dec. 1962). 

Intercomparison of national roentgen and gamma ray ex- 
pos ire-dose standards, H Rt 4 W vckoff, A. Allisy, G. H. Aston, 
G. P. Barnard, W. Hiibner, T. Loftus, and G. Taupin, 
Acata Radiol. 1, No. 1, 57-58 (Feb. 1963). 

X-ray microscopy of polymers by point projection, S. B. 
Newman, Mod. Plastics 40, No. 7, 165-179 (Mar. 1963). 
Electric fields in the ionosphere and the excitation of the red 
line s of atomic oxvgen, L. R. Megill, M. H. Rees, and L. K. 
Droppleman, Planetary Space Sci. 11, 45-56 (Jan. 1963). 

Microwave spectrum and structure of difluoramine, D. R. 
Lide, Jr., J. Chem. Phys. 38, No. 2, 456-460 (Jan. 15, 1963). 

International geophysical calendar for 1963, A. H. Shapley 
and J. V. Lincoln, J. Geophys. Res. 68, No. 4, 1157-1159 
(Feb. 15, 1963). 

Infrared absorption spectra of carbon suboxide and malonoti- 
trile in solid argon matrices, L. L. Ames, D. White, and 
D. E. Mann, J. Chem. Phys. 38, No. 4, 910-917 (Feb. 15, 
1963). 

Interactions matrix element in a shell model, U. Fano, F. 
Prats, and Z. Goldsmith, Phys. Rev. 129, No. 9, 2643-2652 

Mar. 16, 1963). 

Electron impact ionization of atomic hydrogen, S. Geltman, 
M. R. H. Rudge, and M. J. Seaton, Proce. Phys. Soc. 81, 
Pt. 2, Ne. 520, 375-378 (1963). ; 

U.S. participation in international 
McPherson, ASTM Mater. Res. 

Apr. 1963). 

Effective diffusion constant in a polvelectrolyte solution, J. L. 
Jackson and 8. R. Coriell, J. Chem. Phys. 38, No. 4, 959 
968 (Feb. 15, 1962). 

Electron attachment coefficients of some hydrocarbon flame 
inhibitors, T. G. Lee, J. Phys. Chem. 67, 360-366 (1963). 

The measurement of moisture boundary layers and leaf 
transpiration with a microwave refractometer, D. M. 
Gates, M. J. Vetter, and M. C. Thompson, Jr., Nature 
197, 1070-1072 (Mar. 16, 1963 

Nuclear resonance and the hyperfine field in dilute alloys of 


Minesinger, J. Colloid Sci. 18, 


23 236 
r 


standardization, A. T. 
Stds. 3, No. 4, 310-311 


Streever, L. H. Bennett, R. C. 
Appl. Phys. 34, No. 4, pt. 


nickel in iron, R. L. 
Force, and G. F. Day, J. 
1050-1051 (Apr. 1963). 

The history of Pt 27, E. Wichers, Book, Temperature, Its 
Measurement and Control in Science and Industry 3, pt. 1, 
259-262 (Reinhold Publ. Corp., New York, N.Y., 1962). 

Some causes of resonant frequency shifts in atomic beam 
machines. I. Shifts due to other frequencies of excita- 
tion, J. H. Shirley, J. Appl. Phys. 34, 783-788 (Apr. 1963). 

Some causes of resonant frequency shifts in atomic beam 
machines. II. The effect of slow frequency modulation 
on the Ramsey line shape, J. H. Shirley, J. Appl. Phys. 34, 
789-791 (Apr. 1963). 

On the dependence of absorption coefficients upon the area 
of the absorbent material, E. D. Daniel, J. Acoust. Soc. 
Am, 35, No. 4, 571-573 (Apr. 1963). 

The role of the International Union of pure and applied 
chemistry, E. Wichers, J. Chem. Doe. 3, No. 7, 7-11 (1963). 

Note on a subgroup of the modular group, M. Newman and 
J. R. Smart, Proc. Am. Math. Soc. 14, No. 1, 102-104 
(Feb. 1963). 

The orthobaric densities of parahydrogen, derived heats of 
vaporization and critical constants, H. M. Roder, D. E, 
Diller, L. A. Weber, and R. D. Goodwin, Cryogenics 3, 
16-22 (Mar. 1963). 

Melting pressure equation for the hydrogens, R. D. Goodwin, 
Cryogenics 2, No. 6, 1-3 (Dee. 1962). 

Isotopic fractionation of uranium in sandstone, J. N. Rosholt, 
W. R. Shields, and EK. Garner, Science 139, 224 226 
(Jan. 18, 1963). 

The total electron content of the ionosphere at middle lati- 
tudes near the peak of the solar cycle, R. 5. Lawrence, 
D. J. Posakony, 0. K. Garriott, and 8. C. Hall, J. Geophys. 
Res. 68, 1889-1898 (Apr. 1, 1963). 

Present status of our knowledge of atomic transition proba- 
bilities, W. L. Wiese, Proc. Tenth Colloquium Spectro- 
scopic Intern., pp. 37-56 (Univ. of Maryland, College 
Park, Md., 1962). 

Rubber and rubber products, W. P. Tyler and M. Tryon, 
Book, Industrial and Natural Products and Noninstru- 
mental Methods, 6th ed., Standard Methods of Chemical 
Analysis IIB, ch. 43, 2146-2226 (D. Van Nostrand Co., 
Inc., New York, N. Y., 1963). 

The information and oxidation of 
V. R. Deitz and E. F. 
Conference II, 219-232 
England, 1963). 

Pure substance and measurement, E. 
Stds. 1, No. 4, 314-315 (Apr. 1961 

The specific heat at constant 
temperatures from 15 to 90 
B. A. Younglove and D. EK. 
1-5 (Dec. 1962). 

Radiation detectors, L. 
(Mar. 8, 1963). 

Electron spin resonance of gamma-irradiated cellulose, R. E. 
Florin and L. A. Wall, J. Polymer Sci. 1, Pt. A, 1163-1173 
(1963). 

Fundamentals of measurement, A. G. 
Technol. 53, 113-128 (May 1963). 

Optimum estimators of the parameters of negative exponen- 
tial distributions from one or two order statistics, M. M. 
Siddiqui, Ann. Math. Stat. 34, 117-121 (Mar. 1963). 

Electron microscopy studies of the surfaces of magnetic 
recording media, F. Nesh and D. B. Ballard, IEEE Trans. 
Audio AU-11, No. 1, 15-18 (Jan. 2, 1963). 

Kinetics of the acid-catalyzed hydrolysis of acetal in water- 
acetone solvents at 15, 25, and 35°, R. K. Wolford, J. Phys. 
Chem. 67, 632-636 (1963). 

Crystallographic changes with the substitution of aluminum 
for iron in dicaleium ferrite, D. K. Smith, Acta Cryst. 15, 
1146-1152 (Jan. 1963). 

Methods for the analysis of rubber and related products, 
M. Tryon and E. Horowitz, Handb. Analytical Chemistry, 
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sect. 13, pp. 233-256 (McGraw-Hill Book Co., Inc., New 
York, N.Y., 1963). 

Pressure-densit y-temperature relations of 
parahydrogen to 350 atmospheres, R. D. 
genics 3, 12—15 (Mar. 1963). 

Nuclear magnetic resonance in metal powders at low 
peratures, R. J. Snodgrass and L. H. Bennett, J. 
Spectry. 17, No. 2, 53-54 (1963). 

Lunar point-to-point communication, L. E. 
Technology of Lunar Exploration, Progress in 
and Aeronautics 9, 559 (Academic 
York, N.Y., 1963). 

Comment on ‘Parametric behavior of an ideal two-frequency 
varactor,”’ G. F. Montgomery, Proc. IEEE 61, No. 3, 491 
(Mar. 1963 

Distribution of latitude of red ares, E 
Roach, J. Geophys. Res. 68, No. 7, 
1963) 

The following papers were 
Conf. lonosphere, London, 
and The Physical Society, London, England) 

A model of the atmosphere and i ionosphere in the E and 
F1 regions, R. B. Norton, T. E. VanZandt, and J.S. Denison, 
pp. 26-34 

Doppler studies of the ionospheric effects of solar flares, K. 
Davies, pp. 76-83. 

Ionospheriec variations during geomagnetic 
sushita, pp. 120-127 

The location of the irregularities responsible for ionospheric 
scintillation of a radio source, H. J. A. Chivers, pp. 258-266. 

Equatorial spread-F motions, W. Calvert, K. Davies, E. 
Stiltner, and J. T. Brown, pp. 316-322. 

Evidence for field-aligned ionization irregularities between 
100 and 1000 km above the earth’s surface, W Calvert, 
T. E. VanZandt, R. W. Knecht, and G. B. Goe, pp. 324-329 

Experimental observations and _ theoretical calculations 
leading to a model for the lower ionosphere, R. H. Doherty, 
pp. 428-434. 

Radio wave reflections at a continuously stratified plasma 
with electron collision frequency proportional to energy 
and arbitrary magnetic induction, J. R. Johler, pp. 436-445. 

Very low frequency propagation , the earth-ionosphere 
waveguide of non-uniform width, J. R. Wait, pp. 446-451 

United States participation in inte raetion al % andardization, 
A. T. McPherson, Mater. Res. Std. 3, No 310-311 (Apr. 
1963). 

Instabilitv of the equi atorial F la r after sunset, W. 
J. Geophys. Res. 68, 2591-25 13 May 1, 1963). 
Standard potential of the silver-silver chloride electrode and 
activity coefficients of hydrochloric acid in aqueous meth- 
anol (33.4 Wt ©,) with and without added sodium chloride 
at 25°, R. G. Bates and D. Rosenthal, J. Phys. Chem. 67, 

LOSS—LO90 (1963). 

International standardization. A new responsibility of the 
engineer, A. T. McPherson, Natl. Acad. Sci.-Natl. Res. 
Council Div. of Eng. and Ind. Newsletter, No. 26, 
2—4 (June 1, 1963). 

Path antenna gain and comments on 
long-distance tropospheric circuits,’”’ W. J. 
IEEE 51, 847-848 (May 1963). 

Absorption and scattering of photons by deformed nuclei, 
E. G. Fuller, Proce. Second All-Union Conf. Nuclear Reac- 
tions at Low and Intermediate Energies, p. 419, 1960 
(Russian Academy of Science, U.S.S.R., 1962). 

Microwave spectrum of tertiary butyl chloride, A. Comparison 
of tertiary butyl structures, D. R. Lide, Jr., and M. Jen, J. 
Chem. Phys. 38, No. 7, 1504-1507 (Apr. 1, 1963). 

Search for a slow component in alpha ionozation, Z. Bay and 
R. M. Pearlstein, Phys. Rev. 130, No. 1, 223-227 (Apr. 
1963). 

Dielectric friction on a rotating dipole, R. Zwanzig, 
Phys. 38, No. 7, 1605-1606 (Apr. 1, 1963). 

Field-aligned E-region irregularities identified with acoustic 
plasma waves, K. L. Bowles, B. B. Balsley, and R. Cohen, 
J. Geophys. Res. 68, 2485-2501 (May 1, 1963). 


liquid 
Cryo- 


freezing 
Goodwin, 


tem- 


Appl. 


Vogler, Book, 
Astronautics 


533 Press Inc., New 


Marovich 
LSS5 


and F. E. 
I8SSS (Apr. 1, 


published in the Proc. Intern. 
1962 (The Institute of Physics 


Mat- 


storms, 8. 


Calvert, 


Res. 


“properties of 400 Me/s 


Hartman, Proc. 
, 


J. Chem. 





The association of plane-wave electron-density irregularities 
with the equatorial electrojet, R. Cohen and K. L. 
J. Geophys, Res. 68, 2603-2611 (May 1, 1963). 

Quasi-equilibrium theory of mass spectra, H. M. 
and M. 


sowles, 


tosenstock 

Krauss, Book, Mass Spectrometry of Organic Ions, 
po. 1-64 (Academic Press, Inc., New York, N.Y., 1963). 

Determination of differential X-ray photon flux and total 
beam energy, J. 5. Pruitt and H. W. Koch, Book, Nuclear 
Physies 5, ch. 2.8.2, pt. B, 508-553, ed. L. C. L. Yuan and 
C. Wu (Academic Press, Inc., New York, N.Y., 1963). 

Index to the communication sof the ACM, 1958 1962 , W. W 
Youden, Comm. ACM 6, No. 3, 1-32 (Mar. 1963). 

Electrode potentials in fused systems. VI. Membrane 
tials, K. H. Stern, Phys. Chem. 67, 893-895 (1963 

Electron spin resonance spectra of aged, y-irradiated poly- 
styrenes, R. E. Florin, L. A. Wall, and D. W. Brown, J. 
Polymer Sci. pt. A. 1, 1521-1529 (1963). 

Magnetic properties of acetinide element alloys and 
pounds, J. C. Eisenstein, Book, Landolt-Bornstein, 6 Ed. 
2 pt. 9, p. 3.236 (Springer-Verlag, Berlin, Germany 1962). 

Planar twin boundary in aluminum, T. H. Orem, Trans Met. 
Soc. AIME 227, 786-788 (June 1963). 

fadio noise anomalies in August 1958, C. A. 
Geophys. Res. 68, 2719-2726 (May 1, 1963). 

Optical studies of the formation and breakdown of passive 
films formed on iron single crystal surfaces in inorganic 
inhibitor solutions, J. Kruger, J. Electrochem. Soc. 110, No. 
6, 654-663 (June 1963). 

The importance -of environment in fatigue failure of n 
J. A. Bennett, W. L. Holshouser and H. P. Utech, 
Fatigue of Aircraft Structures, pp. 1-18 (Pergamon 
Inc., London, England, 1963). 

Production of Lyman alpha radiation in ion-atom collisions, 
G. H. Dunn, R. Geballe, and D. Pretzer, Phys. Rev. 128, 
No. 5, 2200-2206 (Dec. 1, 1962 

Microwave spectrum of aluminum monofluoride, D. R. 
Jr., J. Chem. Phys. 38, No. 8, 2027 (Apr. 15, 1963). 

Some limits on the effect of coronal self-emission upon the 
ya agg ag coronal ions, R. N. Thomas and C. Pecker, 
Astrophys. J. 137, No. 3, 967-980 (1962 

Two-stream plasma inst ability as a source 
the ionosphere, D. Farley, Jr., Phys. 
No. 7, 279-282 (Apr. 1, 1963). 

tadiation induced polymerization of propylene 
sure, D. W. Brown and L. A. Wall, J. 
1016-1019 (1963). 

Dielectric friction on a moving ion, R. Zwanzig, J. Chem. 
Phys. 38, No. 7, 1603-1605 (Apr. 1, 1963). 
Survey of Rb/Rb* ratios in minerals, W. R. 
Garner, C. E. Hedge, and 8. 8. Goldich, J. 

68, No. 8, 2331-2334 (Apr. 15, 1963). 

Dissociation constant of pyrrolidinium ion and related ther- 
modynamic quantities from 0 to 50°, H. B. Hetzer, R. G. 
Bates, and R. A. Robinson, J. Phys. Chem. 67, 1124-1127 
(1963). 

Introduction to the 
Aeronomy, R. 


poten- 


com- 


Samson, J. 


é tals, 
300k, 
Press, 


Lide, 
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Phys. Chem. 67, 


Res. 
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International Symposium on E« 

Cohen and K. L: Bowles, J. 
68, 2359-2361 (May 1, 1963). 

Numerical computation of time-dependent 
isotopically disordered one-dimensional 
lattices, R. J. Rubin, J. Phys. Soc. 
63-69 (1963). 

The National Bureau of Standards: 
the central theme, C. 8. MeCamy, 
28 (Aug. 1962 

Statistical computation of configuration and free volume of a 
polymer molecule with solvent interaction, J. Mazur and 
L. Joseph, J. Chem. Phys. 38, No. 6, 1292-1300 (Mar. 15, 
1963). 

*Publications for which a price is indicated are available by 
purchase from the Superintendent of Documents, U.S. 
ment Printing Office, Washington, D.C., 20402 (foreign post- 
age, one-fourth additional). Reprints from outside journals 
and the NBS Journal of Research may often be obtained directly 
from the authors. 
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Apparatus for the detection of piezoelectric coupling 

Application to air bearings to an electrodynamic vibration standard 

Aqueous solution, equipment for single crystal growth 

Arrestment error in one-pan two-knife balance systems, minimization 

Automatic multichannel] correlator 


Application to rectanguler and 
* 


precise standardization of noise 


B 


Balance systems, one-pan two-knife, minimization of the 

Blaine, L. R., A far-infrared vacuum grating spectrometer 

Bowman, H. A., H. E. Almer, Minimization of the arrestment error in one- 
pan two-knife balance systems 

Brown, R. F., Jr., Automatic multichanne] correlator 


Buildup source, a power-series buildup factor formulation, 


arre stme nt error 


application. 
C 


Calibrating method, 
Calibration systems, 
accuracy 
Cassidy, E. C., D.H 
speed photography 
Cellulosic materials exposed to thermal radiation, surface flame propagation. 
Colorimeter, thermoelectric, development of filters 
Compensation and measurement, drag with manned 
study 
Conduction, steady 
heat sources 
Correlator, multichannel, automatic 
Crystal growth, single, from aqueous solution, equipment 
Cylinders with multiple continuous line beat sources, 
conduction 


1 standard volt box 
instrument, realistic evaluation of the precision and 


Tsai, New fast-opening, large-aperture shutter for high- 


satellites: feasibility 


state heat, in cylinders with multiple continuous line 


steady state heat 


D 


Darling, W. P., F. E. Washer, 
photographic objective 

Determination of optical p ith difference for a photographic objective 

Determination of residual thiosulfate in processed film 

Development of filters for a thermoelectric colorimeter 

Dimoff, T., B Pa Application of air bearings to an ee a 
vibration standard 

Disk source, a power-series buildup factor formulation, 

Drag - -eecpamaatie and measurement with manned 
stud 

Dun a4 


Determination of optical path difference for ¢ 


ne, 
application 
satellites: fe asibility 


B. L., Method for calibrating a standard 


E 


volt box 


Eisenhart, C., Realistic evaluation of the precision and accuracy of instrument 
calibration system 

Elastic constant—porosity relations for polycrystalline 

Elastic constants of thoria 
dependence 

Electrolytic hydrcgen and deuterium through iron, permeation rates 

Electron beam density of unipotential electron guns at low voltages 

Electron guns, unipotential, at low veltages, 
density 

Emara, S. H., F 
imeter 

Emittance, relation of, to other optical properties 

Equations for the radiofrequency magnetic permeameter 


Equipment for a single crystal growth from aqueous solution 


Far-infrared vacuum grating spectrometer 

Flame propagation, surface, on cellulosic 
radiation 

Film, processed, determination of residual thiosulfate 

Filters for a thermoelectric colorimeter, development. _. 

Flynn, D. R., Radial-flow apparatus for determining the the ‘rmal conduc- 
tivity cf loose-fill insulaticns to high temperatures 

Frenkel, L., Apparatus for the detection of viezoelectric coupling 


G 


Greenspan, L., Oxygen partial pressure warning instrument 

Gross, D., J. J. Loftus, Surface flame propag: 
exposed to thermal radiation 

Guns, unipotential electron, at 
density 


thoria.. 


specimens of varying porosity, temperature 


limitations on electron beam 


P. Teele, Development of filters for a thermoelectric Color- 


materials exposed to thermal 


ition on cellulosic materials 


low 


voltages, limitations on electron beam 





361 


H 
Page 
Harrington, R. D., C. A. Hoer, Parallel reversible permeability measurement 
techniques from 50 ke/s to 3 Ge/s 
High-speed one: iphy, new fast-opening, large-aperture shutter 
Hoer, C. A., D. Harrington, Parallel reversible permeability measurement 
tec hnique $ Siac 50 ke/s to 3 Ge 
, A. L. Rasmussen, Equations for the 
ameter 
Horton, A. T., J. L. Torgesen, 
growth from aqueous solution. 
Howard, F, L., Leak-resistart rotation seal for vacuum applications 
Hubbell, J. H., A power-series buildup factcr formulation. Application 
rectangular and off-axis disk source problems 
Hydrogen and deuterium, electrolytic through iron permeation rates 


I 


s 


radiofrequency magnetic perme 


C. P. Saylor, Equipment for single crystal 


Infrared, far, 
Instrument 
accuracy 
Insulations, thermal conductivity cf loose-fill, high temperature 
spparatus for determining. 
Interferometer, alinement 
Interferometers, large aperture 


vacuum grating 


calibration systems, realistic evaluation of the precision and 


, racial-flow 


with small beam dividers 


K 
Elastic 


Knudsen, F. P., S. Spinner, 
for polycrystalline thoria 
Temperature dependence of the elastic constants of thoria speci 

mens of vi irying porosity 
Kuyatt, C. E., J. A. Simpson, Limitations on electron 
unipotential electron guns at low voltages 


L 


, J. P. Vinti, Drag compensation 
manned satellites: feasibility study- 

Large aperture interferometers with small beam dividers 

Leak-resistant rotation seal for vacuum applications 

Limitatiors on electron beam density of unipotential electron guns 
voltages 

Loftus, J. J., D. Gross, Surface flame propagation on cellulosic m 
posed to thermal radiation. 

M 


Magnetic permeameter, equations for the radiofrequency 
Manned satellites, 
study. 
Measurement 
study 
Method for calibrating a standard bolt box 
Microfilm, processed, stability of residual thiosulfat 
Microwave noise sources, radiometer for precise standardization 
Mills, R. M., New radiofrequency mass spectrometer h 
Cyc le_. 
Minimization of the 
Multichanne!] corre! 


L. Stone, constant—porosity relations 


beam density of 


Langer, R. M and measurement with 


at low 


iterials ex 


drag compensation and measurement with: feasibility 


and compensation, drag, with manned satellites: feasibility 


ving 


irrestment error in one-pan two-knife balance systems 
ator, automatic. 


N 


New fast-opening, large-aperture shutter for high-speed photography 

New radiofrequency mass spectrometer having high duty cycle 

Newton, C. J., Residual stresses and their relaxation cn the surfaces of sections 
cut from plastically deformed steel specimens 

Noise, electrical, analysis of a microwave radiometer for precise stand 
tion 


ardiza 


O 


One-pan two-knife balance systems, minimization of the 
Optical path difference, for a photographic objective 
Optical properties, relation of emittance 

Oxygen partial pressure warning instrument. 


P 


Parallel reversible permeability measurement techniques from 
3 Ge/s. 

Partial pre -ssure warning instrument, Oxygen___. 

Payne, B. T. Dimoff, Application of air bearings to an electrodynami 
vibration - andard __ 

Peary, B. A., Steady state heat conduction in cylinders with m vitiple con- 
tinuous ine heat sources 

Permeability measuremert techniques from 50 ke/s to 3 Ge/s 

Permeameter, magnetic, equations for radiofrequency 

Permeation rates of electrolytic hydrogen and deuterium through iron 

Photography, high-speed, new fast-opening, large-aperture shutter 

Piezoelectric coupling, apparatus for the detection of 

Pitts, J. W., Permeation rates of electrolytic hydrogen 
through iron_. 

Poirting interferometer. ‘on 

Polycrystalline thoria, e lastic constant -porosity relations 

Pope, C. 1., Determination of residual thiosulfate in processed film. 

, Stability of residual thiosulfate in processed microfilm. 

Porosity relations for polycrystalline thoria 


irrestment error 


50 k« 


s to 


and deuterium 





trument, partia Steel specimens, plastically deformed 
lation of emittanc stresses and their relaxation 
Stone, L., S. Spinner, F. P 


for polycrystalline thoria 
onductivity of loose-f 


Page 
surfaces of sections cut from, residual 

101 
Knudse 


Elastic constant—porosity relations 


on cellulosic 


3Y 
, lemperature dependence of the elastic constants of thoria specimens 
of varying porosity 
materials Studies on the tungsten 
5 Surtace flame propa 


henium thermocouple to 2000 
tion 


( : 
ition on cellulosic materials exposed to thermal radia- 
i 


Teele, P, P., S. H, Emara, Developme 1 thermoelectric 
imeter 
Temperature ¢ 


color- 
ndence of the 


varyin® poros 
rhermal conductivity 
flow apparatus fo 
sections cut fron radi 


ela thoria specimens of 
of loose-fill 
r determinin 


surface flame pre 


in h temperatures, radial- 


cellulosic materials 
, tun’sten-rhenium, to 2¢ 
microwave radiometer for 


i¢ colorimeter 


tudies on 


, development 
idual, in processed film, de 

optic il properti hilitv ofr 
udies 


termination 
esidual, in processed micro‘ilm 
son the tungsten-rhenium thermocouplk 
, elastic constant-—porosity 
ng porosity, temypx« 
P : , A. T. Herton, € ylor, Equipment for sinvle crystal 
' Vv eous solution 
Equipment nz Cassidy, New fast-opening, large-aperture shutter 
h-speed p 


C, studi 
electron in 


relations 

iture dependence of the elastic 
qu 
( for 
ions 


on 


, applic 
Dra ¢ ensation ane surement with 
Vv study 


for 
for ¢ 


247 
ilibrat 1 
W 

i ‘ hardson, Analy r 1 microwave radiometer for 
precise t } rdi m of noise source 
Washer, } ; iling, Detert 
, photographic ob 


tical path difference for 


GOVERNMENT PRINTING OFFICE: 1963 








UNITE 
GOVERNMENT FRINTING OFFICE 
CIVISION OF PUBLIC DOCUMENTS 
WASHINGTON, D.C, 20402 


OFFICIAL BUSINESS 


JOURNAL OF RESEARCH >i the National Bureau of Sondre 
C. Engineering and Instrumentation 


Limitations on electron beam density of » ; 
low erie: J. Arol Simpson and C. a“ 


M. Mills. 


A power-series buildup factor formulation. ag pr 
and off-axis disk source problems. J. H. seat 


An alinement interferometer. J. B. ‘Saunders. _ 
Determination of optical path difference 


Francis E. Washer and Walter R. Pee : re 


Emara and Bar P. P. Teele; 


Application of air bearings to an 
T. Dimoff and B. F. Payne, 


Leak-resistant rotation seal for vacuum ¢ 


Studies on the tungsten- shenium : 
Thomas. 


Publications of the National Bu: sau of ‘ a 
Index to volume 67C, January to Decen 


For sale by the Superintendent of SOMERS 
(single evpy). Subseription price: 
UNITED STATES GOV 








